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An extensive theoretical and experimental study is presented of the magnetic susceptibility versus 
temperature xiT) of spin 5 = 1/2 two- and three-leg Heisenberg ladders and ladder oxide com- 
pounds. The xi^r) of isolated two-leg ladders with spatially anisotropic antiferromagnetic (AF) 
Heisenberg exchange was calculated by quantum Monte Carlo (QMC) simulations, with and with- 
out ferromagnetic (FM) second-neighbor diagonal intraladder coupling, and for two-leg ladders 
coupled into a two-dimensional (2D) stacked ladder configuration or a 3D LaCu02.5-type inter- 
ladder coupling configuration. We present accurate analytical fits and interpolations of these data 
and of previously reported related QMC x{T) simulation data for the isolated ladder with spatially 
isotropic exchange, for the 2D trellis layer configuration and for isotropic and anisotropic three-leg 
ladders. We have also calculated the one- and two-magnon dispersion relations for the isolated 
2 X 12 ladder with 0.5 < J' /J < 1, where J is the AF coupling constant along the legs and J' is 
that along the rungs. The exchange constants in the two-leg ladder compound SrCu203 are esti- 
mated from LDA-I-U calculations. We report the detailed crystal structure of SrCu2 03 and of the 
three-leg ladder compound Sr2Cu305. New experimental x(^) data for the two- leg ladder cuprates 
SrCu203 and LaCu02.5, and for the nominally two-leg ladder vanadates CaV205 and MgV205 
which are structurally similar to SrCu203, are presented. These and literature xiT) data for these 
compounds and Sr2Cu305 are modeled using the QMC simulation fits. For SrCu2 03, we find that 
J'/J « 0.5 and J/ks ~ 1900 K, assuming a spectroscopic splitting factor g ~ 2.1, confirming 
the previous modeling results of D. C. Johnston, Phys. Rev. B 54, 13 009 (1996). The interlad- 
der coupling J'" / J = 0.01(1) perpendicular to the ladder layers is found to be very weak and on 
the spin-gapped side of the quantum critical point (QCP) at Jq'cp/J ~ 0.048(2) for J'/J = 0.5. 
Sr2Cu305 is also found to exhibit strong intraladder exchange anisotropy, with J'/J = 0.66(5) and 
J/fcB « 1810(150) K for 5 = 2.1(1). The x{T) data for LaCu02.5 are consistent with J'/J f» 0.5 with 
J/ks ~ 1700 K, again assuming g « 2.1, and with a 3D FM interladder coupling J^^/J ~ —0.05 
which is close to and on the AF ordered side of the QCP at JQCp/i = -0.036(1) for J'/J = 0.5, 
consistent with the observed AF-ordered ground state. The surpisingly strong spatial anisotropy of 
the bilinear exchange constants within the cuprate ladders is discussed together with the results of 
other experiments sensitive to this anisotropy. Recent theoretical predictions are discussed including 
those which indicate that a four-spin cyclic exchange interaction within a Cu4 plaquette is important 
to determining the magnetic properties. On the other hand, CaV205 is found to be essentially a 
dimer compound with AF intradimer coupling 669(3) K {g = 1.96), in agreement with the results 
of M. Onoda and N. Nishiguchi, J. Solid State Chem. 127, 359 (1996). The leg and rung exchange 
constants found for isostructural MgV2 05 are very different from those in CaV2 05, as predicted 
previously from LDA+U calculations. 

PACS numbers: 71.27.-Ha, 71. 70. Cm, 75.30.Et, 75.50.Ee 



I. INTRODUCTION 

Low-dimensional quantum spin systems have attracted 
much attention over the past decade mainly due to 
their possible relevance to the mechanism for the high 
superconducting transition temperatures in the layered 
cuprate superconductors, which contain Cu+^ spin S = 
1/2 antiferromagnetic (AF) square lattice layers. One 
avenue to approach the physics of the two-dimensional 
(2D) AF square lattice Heisenberg antiferromagnet is 



to study how the magnetic properties of AF Heisen- 
berg spin ladders evolve with increasing number of legs 
and/or with coupling between them.EJ The study of such 
systems is also interesting in its own right. The spin- 
ladder field has been motivated and guided by theory. 
Odd-leg ladders with AF leg and rung couplings were 
predicted to have no energy gap ("spin gap") from the 
spin singlet S = ground state to the lowest magnetic 
triplet S ^ 1 excited states (as in the "one-leg" isolated 
chain), whereas, surprisingly, even-leg ladders were pre- 
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dieted ta.laa.ve a spin gap for any finite AF rung cou- 
pling J'.Bl3 For even- leg ladders in which the ratio of 
the rung to leg exchange constants is J' / J < 1, the spin 
gap [dccspaans-. exponentially with increasing number of 
legs.Enj'EJ til'll3 A close relationship of these generic spin 
gap behaviors of 5* = 1/2 even- and odd-leg AF Heisen- 
berg spin ladders was established with AF integer-spin 
and half-integer-spin Heisenheijg .phains, which are gap- 
ful and gapless, respectively MMO A spin gap also oc- 
curs for AF leg coupling if J' is any finite ferromag- 
netic (FM) value, although the dependence of the gap 
on the magnitude of J' is different from the dependence 
when J' is AF; a second-order transition between the two 
spin-gapped ground states occurs when the spin gap is 
zero as the rung corusJing passes from AF values through 
zero to FM values .E3 Subsequent developments in the 
field involved a close interaction of theory and experi- 
ment, mainly on oxide spin ladder compounds. Inter- 
est in the properties of oxide spin ladder materials was 
stimulated by the stripe picture for the high-Tc layered 
cuprate superconductors, in which the doped Cu02 lay- 
ers may be viewed as containing undoped S — 1/2 AF 
n-leg spin ladders sepjarated by domain walls contain- 
ing the doped charges. E5~E£l Undoped non-oxide 5* = 1/2 
two-leg AF spin ladders also exist in nature. The best 
studied of these is Cu2(C5Hi2N2)2Cl4, for which the Cu- 
Cu exchange interactions are much weaker than in the 
layered and spin ladder cuprate compounds, which in 
turn has allowed extensive studies to be done of the low- 
temperature magnetic field-temperatui£ phase diagram 
and associated critical spin dynamics. Ell 

To provide a basis for further understanding of the 
properties of both undoped and doped spin ladders, it is 
important to determine the values of the superexchange 
interactions present in undoped spin ladder compounds. 
The work reported here was motivated by the surpris- 
ing inference by one of us in 1996,E3 based on several 
theoretical analyses of the magnetic susceptibility ver- 
sus temperature x(r) of the two- leg ladder cuprate com- 
pound SrCu2 03, that the exchange interaction J be- 
tween nearest-neighbor Cu spins-1/2 along the legs of 
the ladder is about a factor of two stronger than the ex- 
change interaction J' across the rungs. This conclusion 
strongly disagrees with J' /J > 1 as expected from well- 
established "empirical rules" for superexchange in oxides, 
and has wide ramifications for the understanding of su- 
perexchange interactions not only in cuprate spin ladders 
but also, e.g., in the high-Tc layered cuprate supercon- 
ductors. 

We therefore considered it important to conclusively 
test the modeling results of Ref. To do so, we car- 
ried out extensive high-accuracy quantum Monte Carlo 
(QMC) simulations of x{T) not only for isolated two- 
leg spin S = 1/2 Heisenberg ladders with spatially 
anisotropic exchange as suggested in Ref. |2^, but also 
for both isotropic and anisotropic two-leg ladders cou- 
pled together in a two-dimensional (2D) stacked-ladder 
configuration and in a 3D "LaCu02.5-type" configura- 



tion. These QMC simulations of x(^) ^ire in addition 
to QMC simulations that have already been reported in 
the literature for isolated and spatially isotropic two- and. 
three-leg ladders and, as several of us reported recently,L£l 
for both isotropic and anisotropic intraladder exchange in 
the 2D trellis layer coupled-ladder configuration present 
in SrCu2 03. In order to reliably and precisely model 
experimental x{T) data, we obtained high-accuracy an- 
alytic fits to all of these QMC simulation data, and the 
fit functions and parameters are reported. Two addi- 
tional calculations were carried out. First, we computed 
the one- and two-magnon dispersion relations for two-leg 
2 X 12 ladders in the parameter range 0.5 < J'/ J < 1 
for use in determining exchange constants in two-leg lad- 
der compounds from inelastic neutron scattering data. 
Second, for comparison with our modeling results for 
SrCu203, the exchange constants in this compound were 
estimated using LDA-I-U calculations. 

On the experimental side, we report the detailed 
crystal structure of SrCu203, required as input to our 
LDA-I-U calculations for this compound, along with 
that of the three- leg ladder compound Sr2Cu305. New 
x{T) data are reported for the two-leg ladder cuprates 
SrCu203 and LaCu02.5, and for the two-leg ladder vana- 
dates CaV205 and MgV205 which have a trellis-layer 
structure similar to that of SrCu203. The intraladder, 
and in several cases the interladder, exchange constants 
in each of these materials and in Sr2Cu305 were deter- 
mined by fitting the new as well as previously reported 
x{T) data for one to three samples of each compound us- 
ing our fits to the QMC x{T) simulation data. In each of 
the three cuprate ladder compounds, the intraladder ex- 
change constants were found to be strongly anisotropic, 
with J'/ J w 0.5-0.7, confirming the results of Ref. ^ 

In the following sections we discuss in more detail the 
important previous developments in theory and exper- 
iment on spin ladder oxide compounds relevant to the 
present work, which is necessary to better place in per- 
spective our own comprehensive theoretical and experi- 
mental study on undoped spin ladders, and then give the 
plan for the remainder of the paper. There is an exten- 
sive literature on the theory of doped spin ladders which 
we will not cite or discuss except in passing. 

A. SrCu203 and Sr2Cu305 



Rice, Gopalan and SigristoCI recognized that a class of 
undoped layered strontium cuprates discovered by Hiroi 
and Takano, withj— general formula Srm_iCu„i+i02m 
(to — 3, 5, . . .)j2Hi£3 may exhibit properties character- 
istic of nearly isolated spin 5=1/2 n-leg ladders with 
n = {m + l)/2, due to geometric frusteation between the 
ladders in the "trellis layer" structuretD which effectively 
decouples the ladders magnetically. A sketch of the Cu 
trellis layer substructure of the n — 2 two-leg ladder com- 
pound SrCu203 is shown in Fig. |l|. The above spin-gap 
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Trellis Layer Structure of SrCu203 and CaV205 
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FIG. 1. Sketch showing the basic trellis layer structure of 
the M — Cu and V sublattices in the two-leg ladder compound 
SrCu2 03 and in CaV2 05 and MgV2 05, respectively. The 
figure shows the intraladder leg (or chain) (J) and rung (J') 
exchange couplings and the interladder (J") coupling within 
the trellis layer. The intraladder diagonal coupling (J'^"'^) 
between second-nearest-neighbors within each M4 plaquette 
of each ladder and the nearest-neighbor interladder coupling 
(J'") between pairs of M atoms in adjacent stacked layers are 
not shown. 

predictions were subsequently verified experimentally 
from x(^) measurements on the two- leg ladder com- 
pound SrCu203 which exhibits a spin-gap and for 
the pthcee-leg ladder compound Sr2Cu305 which does 

Normand et aLLj have investigated the ground state 
magnetic phase diagram of the trellis layer in exchange 
parameter space for antiferromagnetic (AF) spin inter- 
actions. They find spin-gap, Neel-ordered and spiral or- 
dered phases, depending on the relative strengths of the 
interactions. Thermodynamic and other properties of the 
trellis layer with spatially isotropic coupling within each 
ladder have been calculated using the Hubbard model 
by Kontani and Ueda.m At half filling, with a ratio of 
interladder {t") to isotropic intraladder (t) hopping pa- 
rameters t" /t = 0.15 and with an on-site Coulomb repul- 
sion parameter U given by U/\t\ = 2.9, they find that a 
pseudogap opens in the electronic density of states with 
decreasing temperature T, and that d-wave superconduc- 
tivity develops in the presence of this pseudogap at lower 
T, similar to behaviors observed for "underdoped" high- 
Tc layered cuprate superconductors. 

Experimental work on spin ladder systems has been 
strongly motivated by such theoretical predictions that 
even-leg ladder systems (with spin gaps) may exhibit a d- 
wave-like superconducting ground state via an electronic 



mechanism when appropriately doped.EHai23~l23 Theoret- 
ical studies suggest that superconductivity may also oc- 
cur in doped three-leg ladderSj-which have no spin gap, 
but with interesting subtleties .EJlZI Thus far, it has not 
proved possible to dope SrCu203 or Sr2Cu305 into the 
superconducing state, although electron doping has been 
achieved bji substituting limited amounts of Sr by La in 
SrCu203.E!l 

Surprisingly, introducing small amounts of disorder in 
the Cu sublattice by substituting nonmagnetic isoelec- 
tronic Zn+^ for Cu+^ in Sr(Cui_a;Zn2;)203 was found to 
destroy the spin gap and induce loiifcpange AF ordering 
at Tn ~ 3-8 K for 0.01 < a; < 0.08.EIE3t1I Specific heat 
measurements above Tn for x = 0.02 and 0.04 indicated 
an electronic specific heat coefficient 7 ~ 3.5mJ/molK^ 
and a gapless ground state; from this 7 value, an (av- 
erage) exchasee constant in the ladders J/Zcb ~ 1600 K 
was derived. Ej Many theoretical studies have been car- 
ried out on site-depleted and otherwise-disordered two-leg 
ladders to interpret these experiments E5L3 The essential 
feature of the theoretical results is that the spin vacancy 
induces a localized magnetic moment around it as well 
as a static staggered magnetization that enhances the 
AF correlations between the spins in the vicinity of the 
vacancy. The enhanced staggered magnetization fields 
around the respective spin vacancies interfere construc- 
tively, resulting in a quasi- long range AF order along the 
ladder, so that even weak interladder couplings are pre- 
sumably sufficient to induce 3D AF long-range order at 
finite temperatures. To our knowledge, no quantitative 
calculations have yet been done of the 3D AF ordering 
temperature Tn in the Neel-ordered regime versus inter- 
ladder coupling strengths for the 3D stacked trellis layer 
lattice spin coupling configuration of the type present in 
SrCu203 either with or without Zn doping. 



B. LaCu02.5 



Another candidate for doping is the two-leg-S|pia-ladder 
compound LaCu02.5 (high-pressure form),L3L3 which 
has the oxygen-vacancy-ordered CaMn02.5 (Ref. [79| ) 
structure. The interladder exchange coupling is evi- 
dently stronger than in SrCu203, since long-range AF 
ordering was observed to occur in LaCu02.5 from ^^Cu 
NMR and muon spin rotation/relaxation (/zSR) 
surements at a Neel temperature Tn = 110-125 K J 
Metallic holepdoaed compounds Lai_2;Sra;Cu02 5 
be formed,BE!lea but superconductivity has not 
been observed at ambient pressure above 
< a; < 0.20 or at high pressures up to 8 GPa.L3lIZl 

Normand and coworkers have cpji'ied out detailed ana- 
lytical calculations for LaCuPj>^E3"E^ The 3D exchange 
coupling topology proposecEil^EZI for this compound is 
shown in Fig. ||. If the leg coupling J = 0, the spin lattice 
is a 2D spatially anisotropic honeycomb lattice, whereas 
if the rung coupling J' = 0, one has a 2D anisotropic 



3 



3D LaCuO^ g-Type Couplings 
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FIG. 2. Exchange coupling topology proposedSS for the 
two- leg ladder compound LaCu02.5. The ladder legs run per- 
pendicular to the page, along which the exchange coupling is 
J (not shown). The figure shows the intraladder rung (J') 
and interladder (J^^) exchange couplings. 

square lattice. —From a tight-binding fit to the LDA 
band structureE3 Normand and Rice inferred that the 
ratio of the rung to leg exchange coupling constants is 
J'/ J ~ 1, with an AE-|interladder superexchange inter- 
action J^^ / J Ki O.25.E3 Then using a mean- field analy- 
sis of the spin ground state, they found (for J'/ J — 1) 
that with increasing J^^ the spin gap disappears at a 
quantum critical point (QCP) at Jqcp/J — 0.121 sep- 
arating the spin-liquid from the long-range AF ordered 
ground state, and thereby inferred that the ground state 
of LaCu02.5 is AF ordered {J^^ / J is on the ordered side 
of Jqcp/<^)- The critical on-site Coulomb repulsion pa- 
rameter Uc necessary to induce AF order was found to 
be given by Uc/W « 0.2 where W is the bandwidth; the 
small value of this ratio is a reflection of the substan- 
tial ID character of the bands, the interladder coupling 
notwithstanding. |— . 

Troyer, Zhitomirsky and UedafiJ confirmed using large- 
scale QMC simulations of x(T) for J' / J ^ 1 that 
Jqcp/J ~ 0.11, and confirmed the predictedE^ behavior 
X oc (up to logarithmic corrections as three dimen- 
sions is the upper critical dimension) at they-Sp QCP. Ad- 
ditional calculations by Normand and RiceE3 in the vicin- 
ity of the QCP predicted that x = x{0)+bT^ with further 
increases in J'^^ /J, where x(0) increases with increase, 
ing (J^^ — Jqcp)/J, consistent with the simulations .eZI 
Comparison of the calculated T^{J^^ /J) with the exper- 
imental results substantiated that 'I^^-/jJ in LaCu02.5 
is only slightly larger than Jqqp/J.E3EZI Other calcula- 
tions, using as input the results of x-ray photoelectron 
spectroscopy measurements, indicate that the interlaiii 
der coupling is ferromagnetic (FM) with \ J^^ /J\ < 0.1,E2l 
rather than AF. However, the same generic behaviors 
near the QCP described above are expected regardless of 
the sign of J^^ , which can be determined from rpagnetic 
neutron Bragg diffraction intensities below Tn;" these 



measurements have not been done yet. 

The possibility of superconductivity occurring in the 
doped system Lai_a;Sr2;Cu02.5 with 0.05 < x < 0.2 
was recently investigated within spin fluctuation theory 
by Normand, Agterberg and Rice, who found that d- 
wave-like superconductivity should occur within this en- 
tire doping range.EJ They suggested that the reason that 
superconductivity has not been observed to date in this 
system may be associated with crystalline imperfections 
and/or with the disruptive influence of the intrinsic ran- 
dom disorder which occurs upon substituting La by Sr. 
They suggested that future improvements in the materi- 
als may allow superconductivity to occur. 

Two conclusions from Refs. ^ and ^ are important 
to the present experimental xC^) studies and modeling. 
First, for intraladder exchange couplings in the proxim- 
ity of the QCP, the existence of relatively weak inter- 
ladder coupling does not change xC^) significantly for 
T > J^^/kB from that of the isolated ladder s [ex cept 
for the usual mean- field shift of x as in Eq. (20a) be- 
low]. Second, the onset of long-range AF ordering has 
a nearly unobservable effect on the spherically-avepaged 
x{T) of polycrystalline samples, as was also foundEil for 
the undoped high-Tc layered cuprate parent compounds. 
Although the studies of Refs. ^ and ^ were carried out 
primarily for ladders with spatially isotropic exchange 
(J' /J = 1), these conclusions are general and do not 
depend on the precise value of J' / J within a ladder. 
They explain both why the experimental x(XXdata for 
LaCu02.5 could be fitted assuming a spin gapllj [Eq. (|l|) 
below] even though this compound does not have one, 
and why this experimental x{T) study did not detect the 
AF ordering transition at Tn. 



C. (Sr,Ca,La)i4Cu2404i 

A related class of compounds with general formula 
A14CU24O41 has been extensively investigated over the 
past several years, especially since large single crystals 
have become available. The structure consists of CU2O3 
trellis layers, as in SrCu203, alternating with Cu02 chain 
layers and A layers, where the ladders and chains are 
both oriented in the direction of the c-axis.ES^EJ The 
Cu02 chains consist of edge-sharing Cu-centered CUO4 
squares with an approximately 90° Cu-O-Cu bond angle, 
so from the Goodenough-Konamori- Anderson superex- 
change rules the Cu-Cu superexchange interaction is ex- 
pected to be -vpakly FM. 

From x(2^)iEju specific heat and polarized, and un- 
polarized neutron diffraction measurements, e3 the un- 
doped compound with An — LagCag shows long-range 
magnetic ordering of the chain Cu spins below Tn = 
12.20(5) K, but, the detailed magnetic structure could 
not be solved.L^ Similar measurements on a single crys- 
tal of the slightly doped compound with An = LasCag 
showed long-range commensurate AF ordering below 
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Tn = 10.5 K with FM alignment of the spins within 
the chains and AF ahgnment between nearest-neighbor 
chains; the ordered Cu moment is ^ 0.2/lxb/Cu with, 
an intrachain FM-ahgned O moment ~ 0.02 /iB/O.EZl 
Long-range AF ordering of the chain- Cu spins has also 
been fomid in single crystals of the system with An 
= Sri4_a:Laa; for X =1-6, 5 and 3 at Tn = 16 K, 12 K 
and 2K, resaectively,L3 and for A14 — Sr2.5Caii.5 at 
Tn ~ 2.1K.EI From Cu NMR and NQR measurements 
on a Sr2.5Caii.5Cu2404i crystal, Ohsugi et al. found that 
the Cu spins in the two-leg ladder trellis layers have 
an ordered moment of only ~ 0.02 /ib, whereas the or- 
dered morftOHt on the magnetic Cu sites in the chains is 
~ 0.56AtB.tH2 

The doped-chain compounds Cao.83Cu02, Sro.73Cu02, 
Cao.4Yo.4Cu02 and Cao.55Yo.25Cu02, containing the 
same type of edge-sharing Cu04-plaquette Cu02 chains 
as in the yli4Cu2404i materials, exhibit long-range 
AF orderinffi_atj-2i^ « 12.2 K, 10.0 K, 29 K and 23 K, 
respectively]lHi2a 

The stoichiometric compound Sri4Cu2404i is self- 
doped; the average oxidation state of the Cu is +2.25, 
corresponding |ta^ a doping level of 0.25 holes/ Cu. 



McElfresh et a/.E23 found that single crystals show highly 
resistive semiconducting behavior with an activation en- 
ergy of 0.18 eV between 125 and 300 K for conduc- 
tion in the direction of the— chains and ladd ers. Va- 
lence-bpnd sum calculations Jl2Zl x{T) (Refs. |95| , |l08| ) and 
ESR£23 measurements indicated that the localized doped 
holes reside primarily on O atoms within the Cu02 
chains, ^nith their spins forming pcnmagnetic Zhang-Rice 
singletal23 with chain Cu spins.Lj The remaining chain- 
Cu spins show a maximum in x(T) at ~ 60-80 K (after 
subtracting a Curie C/T term due to ^ 1.5% of iso- 
lated Cu defect spins) arising from short-range AF or- 
derinft Midj.thti formation of a spin-gap A//cb ~ 100- 



|io^,^|iio| ^jjg ladders do not contribute signifi- 
cantly to x(T) below 300 K due to their larger spin-gap. 
An inelastic neutron scattering investigation on single 
crystals by Regnault et alxll found that at temperatures 
below 150 K, spin correlations develop within the chain 
layers. The data were modeled as arising from chain 
dimers with an AF Heisenberg intradimer interaction 
« 116 K, a FM interdimer intrachain interaction — 12.8K 
and an AF interdimer interchain interaction 19.7 K. Sim- 
ilar measurements on a single crystal by Matsuda et al. 
yielded somewhat different values of these ptkree exchange 
constants.ti^l Specific heat measurementstil from 5.7 to 
347 K and elastic constant measurementalli on a single 
crystal from 5 to llOK showed no evidence for any phase 
transitions. An elastic constant study to 300 K indicated 
broad anamolies in cn and C33 at ~ 110 K and 230 K, pos- 
sibly aaaopiated with charge-ordering effects in the CUO2 
chains. Eij 

In the series Sri4_a;Caa;Cu2404i, substituting iso- 
electronic Ca for . Sr up to a; — 8.4 increases the 
conduct ivity.Pi'Pl At the composition a; = 8, a spin-gap 
A/^B — 140(20) K was found on the chains which was 



modeled as due to AF-coupled dimers comprising about 
29% of the Cu in the chapsJ£3 where /cb is Boltzmann's 
constant. Osafune et aLtia inferred from optical con- 
ductivity measurements that holes are transferred from 



the Cu02 chains to the CU2O3 ladders with-ii 



^asmg x; 
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high pressure enhances this redistribution.ti£ 

The spin excitations in the CU2O3 trellis layers 
in Sri4Cu2404i single crystals were studied usiag in- 
elastic neutroij—SjCattering by Eccleston et alt3l and 
Regnault et aLli£3 and in Sr2.5Caii.5Cu2404i single,crys- 
tals using the same technique by Katano et a/.llH3 The 
spin gaps for ladder spin excitations were found to be 
A/fcB = 377(1) K, 370K and 372(35) K, respectively. 
Essentially the same spin gap (380 K) was obtained by 
Azuma et alSJi from inelastic neutron scattering measure- 
ments on a poly crystalline sample of SrCu203. The good 
agreement among all fouc-spin gap values indicates that 
the hole-doping inferredHj to occur in the CU2O3 trellis 
layer ladders in Sri4Cu2404i and Sr2.5Caii.5Cu2404i has 
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with 

Dagotto et al.t^ recently 



L^oin gap, consistent 



little influence on the ladde: 
NMR results of Imai et al. 
predicted theoretically that lightly hole-doped two-leg 
ladders should exhibit two branches in the lowest-energy 
spin excitation spectra from neutron scattering experi- 
ments, with different gaps for each occurring at wavevec- 
tor (7r,7r). These two branches have evidently not (yet) 
been observed or at least distinguished experimentally. 

Many Cu NMR and NQR studies of the paramag- 
netic shifts and spin f^Y.^"!]^''";^,''^ i^ii^"2i^ii compoui|id^ 
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have been reported.LJ'liiitt£j E£j Melzi andpCarretta, 
Kishine aiid_,Fukuyama,ll23 Ivanov and Ledi23 and Naef 
and Wangi£0 have discussed the spin gaps obtained from 



these measurements and have presented analyses which 
may explain why some of these inferred spin gaps do not 
agree with each other and/or with the spin gaps derived 
independently from other measurements such as inelastic 
neutron scattering and x(T). | — , 
Superconductivity was discovered by Uehara et al.t^B 
under high pressure (3-4.5 GPa) in Sri4„a:Caa;Cu2404i 
for X = 13.6 at temperatures up to ^ 12 K, and 
subsequently confirmed. ollij According to NMR mea- 
surements by Mayaffre et al., a spin gap is absent 
at high pressure in jthp CU2O3 ladders of the super- 
conducting materialJHj a result subsequently studied 



theoretically.Ef3 Metallic interladder conduction within 
the CU2O3 trellis layers occurs at high presswe in a su- 
perconducting single crystal with x = 11.5.t3 These re- 
sults suggest a picture in which the superconductivity 
originates from 2D metallic trellis layers with no spin 
gap. However, the ladder spin gap determined from 
^■^Cu NMR measurements by Mito et al. in a crystal with 
a; = 12 at ambient pressure and at 1.7 GPa, when ex- 
trapolated into the superconducting pressure region, sug- 
gested that the spin gap may persist in the normal stafe 
at the pressures at which superconductivity is found. tl] 
Further, inelastic neutron scattering measurements of 
Sr2.5Caii.5Cu2404i single crystals under pressures up to 
2.1 GPa, which is somewhat below the pressure at which 
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superconductivity is induced, suggested that the spin gap 
does not change significantly with pressure, although ths, 
scattered intensity decreases with increasing pressure. 
Thus whether the superconductivity occurs in the pres- 
ence of a spin gap or not is currently controversial. 

The crystal structure of the ^i4Cu2404i compounds 
can be considered to be an ordered intergrowth of CU2O3 
spin ladder layers and Cu02 spin chain layers, and the 
composition can be written as [A2Cu203]7[Cu02]io-i--^ 
different configuration occurs as [A2Cu203]5[Cu02]7,l12 
corresponding to the overall composition A10CU17O29. 
Single crystals of this phase have been grown at ambient 
pressure with a deficiency (~ 10%) in Cu and in which 
AiQ is a mixture of Sr, Ca, Bi, Y and Pb, and some- 
times Al, which were found to become superconducting 

om both 



iftihe spin gap min- 
itB Substituting V 



at a temperature of 80 K at ambient pressur^ 
resistivity and magnetization measurements 
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D. CaV2 05 and MgV2 05 

The (ii—. vanadium oxide CaV2 05 has a crystal 
structurdllil containing (puckered) V2O3 trellis layers 
with one additional O above or below each V atom, where 
adjacent two-leg ladders are displaced to opposite sides 
of the trellis layer plane. CaV2 05 is a member of the 
^V205 {R = Li, Na, Cs, Ca, Mg) family of compounds, 
each of which exhibits intecesting low-dimensional quan- 
tum magnetic properties .1113 Because the structure is sim- 
ilar to that of SrCu203, and CaV205 was-fiMind to pos- 
sess a spin gap from NMR mcasurements,tl£l this vana- 
dium oxide was suggested to be a possible candidate for 
a. S = 1/2 two-legi ladder compound.lil3 However, Onoda 
and Nishiguchillla found that the x(^) could be fitted 
well by the prediction for isolated dimers, with a spin 
singlet ground state and an intradimer exchange con- 
stant and spin gap J'/^b — A/fce — 660 K. On the 
other hand, Luke et al. concluded from fiSK and mag- 
netization measurements that spiaJreezing occurs in the 
bulk of CaV205 below - 50 K.N This is most likely 
caused by impurities and/or defects in the samples. The 
spin gap is thus evidently destroyed or converted into a 
pseudogap upon even a small amount of doping and/or 
disorder, as was also seen by Azuma et al. in Zn-doped 
SrCu2-a;Zn2,03 as described above. 

For the isostructural compound MgV205 which con- 
tains thp — jifaofi type of Y^Ps trellis layers as in 



CaVsOsJiirN Millet et alSM found from analysis of 
x(T) measurements using Eq. (@) below that A/fce = 
14.8 K, a remarkable factor of 44 smaller than in CaV205. 
/iSR measurements did not show any static magnetic or- 
dering above 2.5 K, and x(r) and high-field (< 30 T) 
magiietazation measurements yielded a value A/Zcb ^ 
17KJ1£3 similar to the result by Millet et al. Inelas- 
tic neutron scattering measurements iiidicated a gap 



A/kB ~ 20 K at a wave vector of (tt, 7r),E£j which how- 



cause strong frustration effects can-s 
imum away from this wavevector.E^ 
by up to 10% Ti introduces a strong local moment Curie 
(C/r) contribution to x(^); however, no long-range AF 
ordering jsjias induced in contrast to lightly Zn-doped 
SrCu203.E£il The exchange interactions in CaV205 and 
MgV2 05 have recently beeii estimated by three of us us- 
ing LDA-t-U calculations ££3 Additional calculations were 
carried out which explain why the spin gaps and excha,iig 
interactions are so different in these two compounds 
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E. Exchange Couplings from Fits of the Uniform 
Susceptibihty by Theoretical Models 

Of particular interest in this paper are the signs, mag- 
nitudes and spatial anisotropics of the exchange inter- 
actions between the transition metal spins in undoped 
spin-ladder oxide compounds. These interactions have 
a direct bearing on the electronic (including supercon- 
ducting) properties predicted for the doped materials 
and are of intrinsic interest in their own right. The 
primary experimental tool we employ here is xi^) mea- 
surements. The_first x(T) measurements on SrCu203 
by Azuma et al£3 were modeled by the low-temperature 
approximation to x(T) of a spin S* = 1/2 tjwa-leg ladder 
derived by Troyer, Tsunetsugu, and Wiirtai^ 



X(T) = A e-A/(fcBT) 



ever need not be the magnon dispersion minimum be- 



(1) 



where a spin-gap A/fce = 420 K was found and /cb is 
Boltzmann's constant. If one assumes spatially isotropic 
exchange interactions J' — J within isolated two-leg lad- 
ders, then using A/J « 0.5 appropriate to this case 
(see Sec. HI) yields J/fce ~ 840 K, about a factor of 
two smaller thaa. in the layered high Tc cuprate par- 
ent compounds.Ell On the other hand, one of us inferred 
from analysis of the value of the prefactor A, which is 
not an adjustable parameter but instead is a function 
of A which in turn is a function of J and J', and from 
fits to the x{T) data by numerical calculations for iso- 
lated ladders, all assuming a g-factor g = 2.1, that a 
strong anisotropy exists between the rung coupling con- 
stant J' and the leg coupling constant J: J'/ J ~ 0.5, 
J/fce ~ 200OJC, for which the spin-gap is similar to that 
cited above£3 If confirmed, which we in fact do here for 
SrCu203 and LaCu02.5 as well as for the three-leg ladder 
cuprate Sr2Cu305, this suppression of J' with respect to 
J is predicted to suppress, pupeccpnducting correlations 
in the doped spin-ladders J3'e1[cjH 

The large spatial anisotropy in the exchange interac- 
tions and the large value of J inferred in Ref. |2^ for 
SrCu203 were very surprising. The Cu-Cu distance 
across a rung in this compound is 3.858 A, and that 
along a leg is 3.934 A (see the crystal structure refine- 
ment data in Sec. VI), so if the nearest-neighbor Cu-Cu 



6 



distance were the only criterion for determining the ex- 
change constants, one would have expected J'/ J > 1, 
not J'/J <C 1. This inference is strengthened when one 
notes that the Cu-O-Cu bond angle across a rung is 180°, 
whereas that along a leg is smaller (174.22°). Further, 
the Cu-Cu distance in the layered cuprates is ~ 3.80 A, 
shorter than cither the rung or leg Cu-Cu distance in 
SrCu203, with similar w 180° Cu-O-Cu bond angles, and 
it is well-established that J/ /cb « X^O K in the undoped 
layered cuprate parent compoundsJ£j so on this basis one 
would expect J and J' in SrCu203 to both be smaller 
than this value, not one of them much larger. On the 
other hand, an O ion in a rung of a ladder in SrCu203, 
with two Cu nearest neighbors, is not crystallographi- 
cally or electronically equivalent to an O ion in a leg, 
with three Cu nearest neighbors, so the respective su- 
perexchange constants J and J' involving these different 
types of O ions are not expected to be identical. Second, 
the experimentally inferred spin-gap is approximately re- 
produced assuming J'/J ~ 0.5 and J/Zcb ^ 2000 K, as 
noted above. Third, the nearly ideal linear Heisenberg 
chain compound Sr2Cu03 with 180° Cu-O-Cu bonds has 
a Cu-Cu exchange cpn.stant-eptimated from x(r) data as 
J/kB = 2150tig^K,ElM:i£3 and from optical measure- 
ments as 2850-3000 K,t£il which arc much larger than 
in the layered cuprates even though the Cu-Cu distance 
along the chain is 3.91 A, significantly larger than the 
3.80 A in the layered cuprates; this itself is surprising. 

Regarding the subject of the present paper, it is impor- 
tant to keep in mind that fitting experimental x(T) data 
by theoretical predictions for a given model Hamiltonian 
can test consistency with the assumed model, but can- 
not prove uniqueness of that model. A recent example 
in the spin-ladder area clearly illustrates this point. The 
compound vanadyl pyrophosphate, {Y0\^2f^ 
("VOPO"), has an orthorhonibic crystal structureE£2 
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which can be viewed crvstallographically as containing 
S = 1/2 two-leg ladders.ll£a However, the x{T) was ini- 
tially fitted by one of us to high precision by the predic- 
tion for the-nS =1/2 AF alternating-exchange Heisen- 
berg chainjl£j a spin-ladder model fit was not possi- 



ble at that time (1987) due to lack of theoretical pre- 
dictions for x(r) of this niodel. When such calcula- 
tions were eventually doneJl£3 it_was found that the 
same experimental x(^) data setli£3 could be fitted by 



the spin ladder model to the same high precision as fot 



the very different alternating-exchange chain model.E££ 



Inelastic neutron scattering measurements on a poly- 
crystalline sample reportedly confirmed the spin-ladder 
model.E£ZI However, subsequent inelastic neutron scat- 
tering results on single crystalSpproved that (VO)2P207 
is not a spin-ladder compound. E£3 The current evidence 
again indicates that (VO)2EsX)7 may be an alternating- 
exchange chain compoundJl£3 although the conjpauBd. 
has continued to^i^rStudied both experimentall}il£3"EI2l 



field magnetization measurements have indicated that 
there are two magnetically distinct types of alternating- 
exchange V chains in (VO)2P207, interpenetrating with 
each other, each with its own spin gap;EZj this find- 



ing has important implications for the interpretation of 
the neutron scattering data. A high-pressure phase of 
(VO)2P207 was recently discovered by Azuma et al. 
which has a simpler structure containing a single type-pf 
S" = 1/2 AF alternating-exchange Heisenberg chain.tZZi 



F. Plan of the Paper 

Herein we report a combined theoretical and experi- 
mental study of the x(T) of S' = 1/2 spin-ladders and 
spin-ladder oxides. Extensive new quantum Monte Carlo 
(QMC) simulations of x{T) are presented in Sec. || for 
isolated two-leg ladders with spatially anisotropic in- 
traladder exchange, including a FM diagonal second- 
neighbor intraladder coupling in addition to J and J', 
and of two-leg ladders interacting with each otlie|n-with 
stacked ladder (for J'/J = 0.5, 1) and proposedy'EZl 3D 
LaCu02.5-type interladder exchange configurations (for 
J'/J = 0.5). We discuss the previous QMC simulation 
data of Frischmuth, Ammon and TroyertZa for the iso- 
lated ladder with J'/J — 1 and for three-leg ladders 
with spatially isotropic and anisotropic intraladder cou- 
plings, of Miyahara et aLcJ for anisotropic two-leg ladder 
trellis layers and of Troyer, Zhitomirsky and UedaEZi for 
isotropic {J'/J = 1) two-leg ladders with 3D LaCu02.5- 
type interladder couplings. In this section we also obtain 
accurate estimates for J'/J = 0.5 and 1 of the values of 
the interladder exchange interactions at which quantum 
critical points occur for the 2D stacked ladder exchange 
coupling configuration and for ladders coupled in the 3D 
LaCu02.5-type configuration. 

A major part of the present work was obtaining a 
functional form to accurately and reliably fit, interpo- 
late and extrapolate the multidimensional (T and one 
or two types of exchange constants) QMC x(T) simula- 
tion data. Four sections are devoted to this topic, where 
we discuss and incorporate into the fit function some of 
the physics of spin ladders. The spin gap of the iso- 
lated two-l eg la dder is part of our general fit function, 
so in Sec. [II A we obtain accurate analytic fits to the 



and theoreticalUiJlI3tj and an alternative 2D model has 



been proposed .1113 Recent ^^^P and ^^V NMR and high- 



reported literature data for the spin gap versus the in- 
traladder exchange constants J and J' . High tempera- 
ture series expansions (HTSEs) of x{T) for isolated and 
coupled isotr opic and anisotropic ladders are considered 
in Sec. |lIIB| . The first few terms of the general HTSE 
for an arbitrary Heisenberg spin lattice with spatially 
anisotropic exchange, which to our knowledge have not 
been reported before, are given and are incorporated into 
the fit function so that the function can be accurately 
extrapolated to arbitrarily high temperatures, and also 
so that the function can be used to reliably fit QMC 
data sets which contain few or no data at high temper- 
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atures. In this section we also give the HTSE to lowest 
order in 1/T for the magnetic contribution to the spe- 
cific heat of = 1/2 spin ladders, and correct an error 
in the literature. The fit function itself that we use for 
most of the fits to the QMC x{T) data is then presented 
and discussed in Sec. [II C . In special cases where suffi- 



cient low-temperature QMC data are not available, a fit 
function containing a minimum number (perhaps only 
zero, o ne or t wo) of fitting parameters must be used, so 
in Sec. Ill D| we consider such functions formulated on 
the basis of the molecular field approximation. The Ap- 
pendix gives the detailed procedures we used to fit our 
new QMC xC^) simulation results and the previously re- 
ported QMC data cited above. Tables containing the fit- 
ted parameters obtained from all fifteen of the one-, two- 
and three-dimensional fits to the various sets of QMC 
data are also given in the Appendix. We hope that the 
general fit function and the extensive high-accuracy fits 
we have obtained will prove to be generally useful to both 
theorists and experimentalists working in the spin ladder 
and low-dimensional magnetism fields. 

Our calculations of the one- and two-magnon dis- 
persion relations for isolated two-leg ladders with 
0.5 < J'/ J < 1, in increments of 0.1, are presented in 
Sec. 0; 



]liese extend the earlier calculations by Barnes 
and RieraMj for J'/J = 0.5, 1 and 2. Additionally we 



discuss the influence of interladder couplings on these 
dispersions. Dynamical spin structure factor S{q, lu) cal- 
culations for J'/J = 0.5, the experimentally relevant ex- 
change constant ratio, are also presented in this section 
and compared with previous related work. These cal- 
culations show that two-magnon excitations should be 
included in the modeling of inelastic neutron scatter- 
ing data when using such data to derive the full one- 
magnon triplet dispersion relation including the higher- 
energy part. Our calculations of the intraladder and in- 
terladder exchange constants in SrCu2 03, obtained using 
the LDA-f-U method, are presented in Sec. ^ 

We begin the experimental part of the paper by pre- 
senting in Sec. VI a structure refinement of SrCu2 03, 
necessary as input to t he L DA-fU calculations, and also 
of Sr2Cu305. In Sec. VII, we present our new experi- 
mental x{T) data for SrCu203, LaCu02.5, CaV205 and 
MgV205 and estimate the exchange constants in these 
compounds and in Sr2Cu305 by modeling the respective 
x{T) data using our fits to the QMC x{T) simulation 
results. 



The paper concludes in Sec. VIIl with a summary and 
discussion of our theoretical and experimental results, 
their relationships to previous work and a discussion of 
how the presence of four-spin cyclic exchange, as pre- 
sented in the literature, can affect the magnetic proper- 
ties of spin ladders and the exchange constants derived 
assuming the presence of only bilinear exchange interac- 
tions. Further implications of the presence of this cyclic 
exchange interaction are also discussed. 



II. QUANTUM MONTE CARLO SIMULATIONS 

Throughout this paper, the Heisenberg Hamiltonian 
for bilinear exchange interactions between spins is as- 
sumed. 



H — ^ Jij Si ■ Sj 



(2) 



<ij> 



where Jij is the exchange constant linking spins Si and 
Sj, Jij is positive (negative) for AF (FM) couphng and 
the sum is over distinct exchange bonds. For notational 
convenience, we define the reduced spin susceptibility x* > 
reduced temperature t and reduced spin gap A* as 



^spin^ma 



J" 



A* = 



J" 



(3a) 



(3b) 



(3c) 



where x^^"^ is the magnetic spin susceptibility, J™^'^ is 
the largest (AF) exchange constant in the system, N is 
the number of spins, g is the spectroscopic splitting fac- 
tor, fiB is the Bohr magneton and /cb is Boltzmann's 
constant. All of the QMC simulations presented and/or 
discussed here are for spins S = 1/2. We summarize be- 
low the definitions of the exchange constants to be used 
throughout the rest of the paper: 



J 
J' 

jdiag 

J" 

J'" 
J3D 



Nearest — neighbor, leg 2 x 

Nearest — neighbor, rung 1 x 

Second — neighbor, 2 x 
diagonal intraladder 

Trellis layer, interladder 2 x 

Stacked ladder, interladder 2 x 



(4) 



3D interladder, 
LaCu02.5 — type 



2 X 



The uniform susceptibilities x(T) of S" = 1/2 Heisen- 
berg spin ladder models were simulated using the contin- 
uous time versi|Oii,of the quantum Monte Carlo (QMC) 
loop algorithm.llIj This algorithm uses no discretization 
of the imaginary time direction and the only source of 
systematic errors is thus finite size effects. The lattice 
sizes were chosen large enough so that these errors are 
much smaller than the statistical errors of the QMC sim- 
ulations. The simulations of the trellis layer suffer from 
the "negative sign problem" caused by the frustrating 



interladder interaction J" . Improved estimators^'' 
used to lessen the sign problem in this caser3 



were 
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Magnetic Susceptibility of Isolated 2-Leg Ladders 

'^■^^ ~^ Uniform Chain; 

J'/J = 

0.10 



0.05 



0.00 
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k^T/J 



Uniform Chain 
J'/J = 
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FIG. 3. (a) Quantum Monte Carlo (QMC) x(r) simulation 
data for isolated two-leg ladders with J'/J = 0.1, 0.2, 0.3, 
0.4, 0.5, 0.6, 0.7, 0.8 and 1.0 (symbols, top to bottom). Ad- 
ditional simulation data for J'/J = 0.25, 0.35, 0.45 and 0.55 
are not shown. Trbe data for J'/J = 1 are the QMC data of 
Frischmuth et ai.tlj The set of solid curves through the data 
points is a two-dimensional fit to all the QMC data. Atea 
shown are the fit for the 5* = 1/2 uniform Heisenberg chainc£j 
(J'/J = 0) and an interpolation curve for J'/J = 0.9. (b) Ex- 
panded plots at low temperatures of the data, now including 
error bars, and curves in (a). 



A. Isolated Ladders 

X* (t) was simulated for isolated two-leg ladders of size 
2 X 200 spins 1/2 with J'/J = 0.1, 0.2, 0.25, 0.3, . . ., 0.6, 
0.7 and 0.8 with maximum temperature range t = 0.01 
to 3.0, comprising 348 data points; here, J'^'^^ = J. A 
selection of the results for i < 2 in J'/J increments of 0.1 
is shown as open and filled symbols in Fig. |^(a) along4vith 
the QMC simulation results of Frischmuth et aZJlZa for 
J'/J = 1 (30 data points from t — 0.05 to 5). Expanded 
plots of the data, now including error bars, are shown in 
Fig. ||(b), where the error bars are seen to be on the order 
of or smaller than the size of the data point symbols. 



Isolated 2-Leg Ladders 
with Ferromagnetic Diagonal Couplings 
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FIC. 4. (a) Quantum Monte Carlo simulations of the mag- 
netic spin susceptibility x vs temperature T for isolated 
S — 1/2 two-leg ladders with J'/J = 0.4 and 0.6 and 

each with ferromagnetic second-neighbor diagonal couplings 
jdiag^j ^ _Q Q5 g^^^ (-j^^j^ £qj. jdiag^j ^ „o.iii are 

not shown) . The set of solid curves is a three-dimensional fit 
to all the data, including those not shown for J'/J — 0.45, 
0.5, 0.55 and 0.65, each with J'^'='s/j ^ q, -0.05, -0.1 and 
—0.111. (b) Expanded plots at low temperatures of the data 
and fit in (a), where the error bars on the data points are also 
plotted. 

We have also simulated x*it) for isolated two-leg lad- 
ders with J'/J = 0.4, 0.45, 0.65, each with fer- 
romagnetic diagonal coupling J'^^^^/J = —0.05, —0.1 
and —0.111, and with maximum temperature range t = 
0.025 to 2, comprising 457 data points. This FM sign 
of the diagonal intraladder coupling was motivated by 
the LDA+U results reported below in Sec. The 
QMC results with error bars for J'/J — 0.4 and 0.6 and 
jdiag / J _ —0.05 and —0.1 are shown as the symbols in 
Fig. H along with our results above for J'/ J = 0.4 and 0.6 
and J'^'^s — 0. The x*{t) is seen to be only weakly af- 
fected by the presence of J'^^'^^. Hence one expects that 
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Magnetic Susceptibility of Isolated 2-Leg Ladders 
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FIG. 5. Quantum Monte Carlo simulations of x{T) for iso- 
lated ladders with J/J' = 0.1, 0.2, 0.3, 0.4, 0.5, 0.7, 0.9 



and 1.0 (symbols, top to 



fee- 



are from Frischmuth et al. 



;tom). The data for J/J' = 1 
In all cases, the error bars are 
smaller than the data symbols. The set of solid curves is 
a two-dimensional fit to all the data. The solid curve for 
J/J' = is a fit to the exact theoretical x{T) for the isolated 
dimer. Interpolation curves for J/J' = 0.6 and 0.8 are also 
shown. 

fits of experimental x(T) data for spin ladder compounds 
by the simulations will not be capable of determining 
J^i'^^g quantitatively if | J'^'^^s/J] < 1. 

For stronger interchain couplings J'/J > 1, 'X*{t) was 
simulated for isolated two-leg ladders with J/J' — 0.1, 
0.2, 0.3, 0.4, 0.5, 0.7 and 0.9 with maximum temperature 
range t — 0.06 to 1.5, comprising 119 data points; here, 
in Eqs. (^) one has J'^^^^ = J'. The results with error 
bars are shown as open and filled symbols in Fig. |^, alaaj 
with the QMC simulation results of Frischmuth et al.Hl 
for J/J' = 1. The susceptibihty ;^*.dimor^^) ^^^^^^ isolated 
antiferromagnetically-coupled S = 1/2 dimer (J/J' — 0) 
is shown for comparison, where 



*, dimer 



t{3 + eyt) 



(5) 
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FIG. 6. (a) Quantum Monte Carlo x data (symbols) vs T 
for the trellis layer with J'/J = 0.5 and 1, where the set of 
solid curves for J" 7^ is the molecular-field-theory prediction 
obtained from the fit for J" = with no adjustable param- 
eters, (b) Expanded plots at low temperatures of the data 
and fit in (a). The error bars for all data points are shown in 
both (a) and (b). 

The x*{t) for the S = 1/2 uniform Heisenberg chain 



(J'/J = 0) was obtained essentia, 
Affleck and Takahashi in 1994.t8i 



ly exactly by Eggert, 
A fit ta 



recently 

refined numerical calculations of KliimpeiH^d'E^ for this 
chain in the temperature range 0.01 < t < 5 is shown in 
Fig. H for comparison with the data; this high-accuracy 
(15ppm rms fit deviation) seven-parameter analytical 
fit to these calculated data, using the fitting sche me in 
Sec. Ill below, was obtained ("Fit 1") in Ref. 182. 



B. Coupled Ladders 

1. Trellis Layer Interladder Interactions 

Miyahara et al^ carried out QMC simulations of x* (t) 
for trellis layers with J'/J = 0.5 (64 data points) and 1 
(72 data points) over a maximum temperature range 
0.1 < t < 1.5, with trellis layer interladder couplings 
J"/ J = -0.2, -0.1, 0.1 and 0.2 for each J'/J value. 
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FIG. 7. (a) Quantum Monte Carlo simulations of the mag- 
netic spin susceptibility x of the S = 1/2 trellis layer vs 
temperature T for intraladder couplings J /J' = 0.1 and 
terladder couplings J" /J' = -0.2,-0.1, 0, 0.1 and 0.2 
(b) Expanded plots of the data in (a) at low temperatures. In 
(a) and (b), the error bars on the data points are all smaller 
than the data point symbols. The set of solid curves is a 
three-dimensional fit to all these data together with those in 
Fig. |, which total 162 data points for J" /J' ^ 0. 



and for additional exchange parameters J" /J = ±0.5 
which we do not discuss here. The results are shown in 
Fig. ^, along with the above isolated ladder results for 
J'/ J = 0.5 and 1 and J" = 0. We also show their QMC 
simulations for the strong-coupling regime over the max- 
imum t range 0.08 < t < 1.5, for the same values of 
J" /J and with J/J' = 0.1 (84 data points, Fig. 0) and 
J/ J' = 0.2 (78 data points, Fig. ^). The results are seen 
to be quite insensitive to the frustrating interladder ex- 
change interaction J". In addition, for ladders with spa- 
tially isotropic exchange, Gopalan, Rice and Sigrist have 
deduced that the spin gap is nearly independent of J" 
for weak coupling.!!! Thus one expects that fits of exper- 
imental x{T) data by our fits to the QMC data will not 
be able to establish a quantitative value of J" for trellis 
layer compounds with weak interladder interactions. 
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FIG. 8. (a) Quantum Monte Carlo magnetic spin suscep- 
tibility X data (symbols) for the S = 1/2 trellis layer vs 
temperature T for intraladder couplings J/ J' — 0.2 and 
terladder couplings J" /J' = -0.2,-0.1, 0, 0.1 and 0.2 
(b) Expanded plots of the data in (a) at low temperatures, 
where the error bars are also plotted. The set of solid curves 
is a three-dimensional fit to all these data together with those 
in Fig. |. 



2. Stacked Ladder Interladder Interactions 

Another interladder coupling path, with exchange con- 
stant J'", is from each spin in a ladder to one spin in each 
of two ladders directly above and below the first ladder, 
a 2D array termed a "stacked ladder" configuration. We 
have carried out QMC simulations of x*(t) for J'/ J = 0.5 
(96 data points) and 1 (94 data points) over a maximum 
temperature range 0.02 < i < 1 with AF stacked ladder 
couplings J"7 J = 0.05, 0.1, 0.15 and 0.2 for each J' /J 
value, and also for J'" /J = 0.01, 0.02, 0.03 and 0.04 for 
J'/ J = 0.5 (106 data points). The results for J'/ J = 0.5 
are shown in Fig. |[ A log-log plot of the low-t data for 
J'" / J — 0.01 to 0.05 from F ig .. Igl lis-slnown separately in 



Fig. hrt. According to theory,t 



the quantum crit- 



ical point (QCP) separating the spin-gapped phase from 
the AF ordered phase in a 2D system is characterized by 
the behavior oc t at low t. A comparison of the data 
in Fig. O with the heavy solid line with slope 1 indicates 
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FIG. 9. (a) Quantum Monte Carlo magnetic spin suscep- 
tibility X data (symbols) for stacked 5 = 1/2 2-leg ladders 
vs temperature T for intraladder couplings J'/J = 0.5 and 
stacked interladder couplings J'" / J = 0, 0.05, 0.10, 0.15 
and 0.2. (b) Expanded plot of the data in (a), plus additional 
data for J'" /J = 0.01, 0.02, 0.03 and 0.04, at low temper- 
atures. In (a) and (b), the error bars are smaller than the 
data symbols. The solid curves are two two-dimensional fits 
to the data in the gapped regime (J"'/J = 0-0.04) and gap- 
less regime {J'" / J = 0.05-0.2), respectively. Note the unique 
crossing point of all the curves at ksT/J « 0.16, where x is 
independent of J'" / J. 

that a QCP occurs for J'/J = 0.5 at 0.04 < Jqcp/J < 
0.05. In order to obtain a more precise estimate of J'" / J 
at the QCP, in Fig. |ll| we plot x*{t = 0) vs J'"/ J, where 
the x*{'t = 0) values were determined by fits to the data 
as described later. By fitting these x*it = 0; J'" I J) data 
by various polynomials, such as the third order polyno- 
mial shown as the solid curve in the figure, and by non- 
integer power laws, we estimate Jqqp/J — 0.048(2) with 
conservative error bars. 

The x*{t) data for isotropic {J'/J — 1) stacked ladders 
are shown in Fig |l2|, along with the above isolated ladder 
results for J'/ J — 1 and J'" = 0. A QCP is seen to occur 
at Jq'cp/ J ~ 0.16. This value of Jqcp/ J is much smaller 
than the various values 0.32(2) (Ref. 0.43 (Ref. ||) 
and 0.30 (Ref. 187) inferred for spatially isotropic ladders 
arranged in a nonfrustrated flat 2D array, because the 
interladder spin coordination number for the 2D stacked 
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FIG. 10. Log-log plot of the low-temperature QMC mag- 
netic spin susceptibility x data (symbols) for stacked S = 1/2 
2-leg ladders vs temperature T for intraladder couplings 
J'/J = 0.5 and interladder couplings J'" /J = 0.01 to 0.05 
(bottom to top data sets) from Fig. |^. The error bars (not 
shown) are smaller than the data point symbols. The lines 
connecting the data points are guides to the eye. Compari- 
son of the data with the heavy line with a slope of 1 indicates 
that the quantum critical point occurs for J'" /J between 0.04 
and 0.05. 
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FIG. 11. Magnetic spin susceptibility at zero temperature 
x{T = 0) data (•) for gapless stacked 5=1/2 2-leg ladders vs 
stacked interladder coupling J'" / J for intraladder couplings 
J'/J — 0.5. The solid curve is an exact third-order poly- 
nomial fit to the data. The quantum critical point (QCP) 
occurs when x{T = 0,J"'/J) 0, from which we find 
Jqcp /J = 0.048(2). 

ladder configuration is two whereas for ladders arranged 
in 2D flat layers it is only one. On the other hand, 
our Jqcp/ J for J' / J = 1 is somewhat larger than that 
(w 0.11, Refs. ||,|^ for the 3D LaCuOs.s-type coupHng 
configuration (see below), even though the interladder 
spin coordination number is the same; here the effects of 
dimensionality evidently come into play, with fluctuation 
effects generally being stronger in lower dimensional sys- 
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FIG. 12. Quantum Monte Carlo magnetic spin susceptibil- 
ity X data (symbols) for stacked S = 1/2 2- leg ladders vs 
temperature T for intraladder couplings J' / J — 1 and in- 
terladder couplings J'" /J = 0, 0.05, 0.1, 0.15 and 0.2. The 
system has a spin gap for J'" /J = 0, 0.05, 0.1 and 0.15 but is 
gapless for J'" /J = 0.2. The set of solid curves through the 
data for J'" /J — 0, 0.05, 0.1 and 0.15 is a two-dimensional fit 
to these data. A separate fit was obtained for J'" /J = 0.2, 
shown as the solid curve through those data. Note the unique 
crossing point of all the curves at k^T / J « 0.31, where x is 
independent of J'" / J . 

terns (other factors being equal). 

Another interesting aspect of the QMC data for the 
stacked ladders is the striking well-defined crossing points 
of x*(t) versus J'" / J at t « 0.16 for fJJ = 0.5 in Fig. | 
and at t k. 0.31 for J' / J = 1 in Fig. M. Similar crossing 
points in the specific heats of strongly correlated electron 
systems versus some thermodynamic variable (e.g., pres- 
sure, magnetic field, interaction parametej) have been 
pointed out by VoUhardt and coworkers .t23ll2j 

We have also obtained x* (t) data for strong interchain 
intraladder couplings J/J' = (42 data points), 0.1 (39 
data points) and 0.2 (38 data points) over a maximum 
temperature range 0.06 < t < 1.5 with stacked ladder 
coupfings J'"/ J' = 0.1 and 0.2 for each J/ J' value. The 
results are shown in Fig. along with x*(0 f^^' the 
isolated dimer from Eq. (H) and the above isolated ladder 
results for J/ J' = 0.1 and 0.2 and J'" = 0. Over these 
exchange interaction ranges, the spin gap persists and a 
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FIG. 13. (a) Quantum Monte Carlo magnetic spin suscep- 
tibility X data (symbols) for stacked S = 1/2 2-leg ladders 
vs temperature T for intraladder couplings J/ J' =0, 0.1 
and 0.2 and interladder couplings J'" / J' = 0, 0.1 and 0.2 for 
each J/ J' value, (b) Expanded plot of the data in (a) at low 
temperatures. The error bars for each data point are shown in 
both (a) and (b). Note that the pairs of data sets with a fixed 
value (0.1, 0.2 or 0.3) of J'" / J' + J / J' are nearly coincident, 
and thus all these data closely follow the prediction of molecu- 
lar field theory. The set of solid curves is a three-dimensional 
fit to all the data. 



QCP is not traversed. Note that the pairs of data sets 
with a fixed value (0.1, 0.2 or 0.3) of J'"/ J' -I- J/ J' are 
nearly coincident, and thus closely follow molecular field 
theory which predicts (see Sec. IIIDD that x'{i) of cou- 
pled dimers only depends on the sum of the exchange in- 
teractions between a spin in a dimer and all other spins 
outside the dimer. 
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FIG. 14. (a) Quantum Monte Carlo magnetic spin 
susceptibility x data (symbols) for isolated {J'^^/J 
= 0) and three-dimensionally ferromagnetically-coupled 
(J™/J = -0.02, -0.04, -0.06, -0.08 and -0.1) two-leg 
ladders vs reduced temperature t = ksT/J for intraladder 
couplings J'/J = 0.5. The 3D interladder coupling topol- 
ogy is that proposed for LaCu02.5. (b) Expanded plots at 
low temperatures of the data in (a) plus additional data for 
J^°/J = -0.025, -0.03 and -0.035. In (a) and (b), the error 
bars on the data are smaller than the data symbols. The two 
sets of solid curves are two fits to the data with and without 
spin gaps, respectively. A spin gap occurs for J^^ / J — 0, 
-0.02, -0.025, -0.03 and -0.035, and the set of solid curves 
through these data is a two-dimensional [t, J^^ /J) fit at 
fixed J'/J = 0.5. For J^^/J = -0.04, -0.06, -0.08 and 
—0.1, the magnetic excitation spectra are gapless, and the 
set of solid curves through these data is a three-dimensional 
{t, J'/J, J'^^ / J) fit to these data together with the data for 
J'/J = 0.6, 0.7, 0.8, 0.9 and 1.0 in Figs, [g] and [l| with J™/J 
values in the gapless regime given in Sec. |A4 . Note that x 
is nearly independent of J^^ / J at t ~ 0.21, even though the 
different data sets do not cross. 



3. 3D LaCu02.5-Type Interladder Interactions 



Additional simulations of were performed which 
incorporated the nonfrustrated 3D interladder coupling 




k^T/J 

FIG. 15. Quantum Monte Carlo magnetic spin suscepti- 
bility X data (symbols) from Fig. |l^(b), replotted on log-log 
scales. The light lines are guides to the eye. At the quan- 
tum critical point (QCP), x oc T'^ at low T as indicated 
by the heavy line with slope 2. A comparison of this line 
with the data indicates that the QCP occurs in the range 
0.035 < t J^°/J[ < 0.040. 

configuration proposecS0 for the two-leg ladder com- 
pound LaCu02.5, in which each spin is coupled by ex- 
change constant J^^ to one nearest-neighbor spin in each 
of two adjacent ladders diagonally above and below the 
first ladder in adjacent layers, respectively. Simulations 
are reported here for J'/J = 0.5 and J^^/J — 0.2, 
0.1, 0.05, -0.02, -0.025, -0.03, -0.035, -0.04, -0.06, 
—0.08 and —0.1 over the maximum temperature range 
0.021 < t < 3 (204 data points), for a 3D lattice of 
size 6x6 ladder^ x 40 spins. The simulations for ferro- 
magnetic (FM, negative) interladder couplings were mo- 
tivated by the recent findings of Mizokawa et al^B men- 
tioned in the Introduction. The data for FM couplings 
at t < 2, plotted in Fig. |l^, indicate a loss of the spin 
gap for 0.035 < | J^°/J| < 0.040. At the QCP in a 3D 
systepii^* is predicted theoretically to be proportional 
to f^.E30 We have plotted our y*(i) simulation data on 
double logarithmic axes in Fig. ^ where by comparison 
of the data with the heavy line with slope 2, the quantum 
critical point for J'/J = 0.5 is indeed seen to be in this 
range. j_. 

The previously reporteoEl x*{t) simulation data for 
J'/J = 1 and = 0.05, 0.1, 0.11, 0.12, 0.15 and 0.2 

(a total of 169 data points), isihich show a quantum crit- 
ical point at Jqcp/ J « 0.11, Ell are shown in Fig. 16 



Since one expects that upon approaching the QCP 
from the AF-ordered side that x*(i = 0) — > 0, to deter- 
mine more precisely the QCPs we have plotted x*(t = 0) 
vs for each of J'/J = 0.5 and 1 in Fig. |l]. The 

extrapolated x*{t = 0) values were determined by fits to 
the x*{t) data described in Sec. Ill below and in the Ap- 
pendix. From exact polynomial fits to the data in Fig. |l^, 
we find J^p/J = -0.036(1) for J'/J = 0.5 and J^p/J 
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FIG. 16. (a) Quantum Monte Carlo magnetic spin sus- 
ceptibility X data (symbols with error bars) for isolated 
{J"^^ / J = 0) and three-dimensionally antiferromagneti- 
cally-coupled (J^°/J = 0.05, 0.1, 0.15, 0.2) two-leg ladders 
vs reduced-temperature t = k^T/ J for intraladder couplings 
J'/J — I.lZI The interladder coupling topology is that pro- 
posed for LaCu02.5. (b) Expanded plots of the the data 
in (a) and additional data for J'^^/J = 0.11 and 0.12. In 
(a) and (b), the data for J^°/J = 0, 0.05, 0.1 and 0.11 are 
in the gapped regime, whereas those for J^°/J = 0.12, 0.15 
and 0.2 are in the gapless regime. The set of solid curves 
for the gapped regime is a two-dimensional {t, J"^^ / J) fit to 
these data, whereas the set of solid curves for the gapless 
regime is a three-dimensional (t, J'/J, J^^ / J) fit to these 
data together with the data for J'/J — 0.5, 0.6, 0.7, 0.8 and 
0.9 in Figs, and |l^ in the gapless regime with J^^ /J values 
given in Sec7 ]A 4 Note the crossing point of all the curves at 
t ~ 0.32, where x is nearly independent of J^^ /J. 



= 0.115(1) for J'/J = 1. There also exist QCPs for the 
opposite sign of J^^ / J in each case, respectively. 

We report here additional x*(i) simulation data for 
intermediate values of J'/J — 0.6, 0.7, 0.8 and 0.9, each 
with = 0.1, 0.15 and 0.2 at temperatures 0.3 < 

t < 3 {a total of 325 data points). A selection of these 
data for J^^/J = 0.1 and 0.2 are plotted in Figs. |l|(a) 
and [l^(b), respectively. 
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FIG. 17. Quantum Monte Carlo magnetic spin suscepti- 
bility data extrapolated to zero temperature x{T = 0) (sym- 
bols) for gapless two- leg ladders vs interladder coupling J"^^ / J 
and intraladder couplings J'/J = 0.5 and 1. The interlad- 
der coupling topology is that proposed for LaCu02.5. The 
solid curves are exact polynomial fits to the respective sets 
of data. The quantum critical point (QCP) is determined by 
the condition x{T = 0) = 0, yielding Jqcp/J = -0.036(1) 
for J'/J = 0.5 and J^p/J = 0.115(1) for J'/J = 1. There 
is a QCP for J^^ /J of the opposite respective sign for each 
value of J'/J. 



C. Three-Leg S = 1/2 Ladders 

The x*{t) for S — 1/2 n-leg ladders with n = 1, 
2, . . . , 6 and isotropic exchange {J'/J — 1) and for 
n — 3 with spatially anisotropic exchange was com- 
puted using QMCL simulations for 0.02 ^ t < 5 by 
Frischmuth et a/JlZj As noted in the Introduction, the 
even-leg ladders exhibit a spin-gap but the odd-leg lad- 
ders do not. In this paper we will be fitting experimental 
x{T) data for the three-leg ladder compound Sr2Cu305. 
The QMC x* (t) data for three-leg laddeps-with both spa- 



tially isotropic and anistropic exchangetlS are shown in 
Fig. |l^. The x*(0 is seen to be sensitive to the value 
of J'/J. As a consequence, we expect to be able to ob- 
tain an accurate value of J'/J for Sr2Cu305 from fits to 
the experimental xiJ") data. In contrast, Kim et al. have 
found from QMC simulations that x*(^) is very insensi- 
tive to the interladder coupling-tietween isotropic three- 
leg ladders arranged in a layer ,E2ZI although it should be 



noted that they calculated x* {t) for nonfrustrated inter- 
ladder couplings and not for the frustrated trellis layer 
interladder coupling configuration present in Sr2Cu305. 



III. FITS TO THE QMC x* (t) SIMULATION 
DATA 

In order to precisely fit experimental x(^) data by the 
QMC x*{t) simulations, it is essential to first obtain accu- 
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FIG. 18. Quantum Monte Carlo magnetic spin susceptibil- 
ity X vs temperature T data (symbols) for three-dimensionally 
antiferromagnetically-coupled gapless two-leg ladders with in- 
traladder couplings J'/J = 0.5 to 1 and interladder couplings 
J^°/J = 0.1 (a) and 0.2 (b). The interladder coupling topol- 
ogy is that proposed for LaCu02.5. The set of solid curves 
in (a) and (b) is a single three-dimensional fit to these data 
together with additional data (not shown) for J^^ /J = 0.05 
and J'/J = 0.6, J^°/J = 0.15 and J'/J = 0.6-0.9, and the 
additional data for J'/J — 0.5 and 1.0 in Figs. ^ and ^ in 
the gapless regime with J^^ / J values given in the respective 
captions. Extrapolations of the fit to low temperatures are 
also shown. Note that the extrapolations for J^^ / J = 0.1 
and J'/J = 0.8 and 0.9 are negative (unphysical) at the low- 
est temperatures. 
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FIG. 19. (a) Quantum Monte Carlo (QMC) magnetic spin 
susceptibility x temperature T data (symbols) for three-ilea 
ladders with spatially anisotropic intraladder exchange. tl3 
The set of solid curves through the data points is a global 
two-dimensional fit. Also shown are extrapolations of the 
QMC data fits for 0.4 < J'/J < 0.7 by a factor of about 
three to higher T, and two interpolation curves for J' /J = 0.8 
and 0.9. The data for the lipear chain were calculated by Eg- 
gert, Affleck and TakahashiJifj (b) Expanded plots of the data 
and fit at low temperatures. 

necessary to use a fit function, containing a minimum 
number (even zero) of fitting parameters, derived from 
the molecular field theory for coupled subsystems. We 
therefore also discuss such fit functions. 



rate analytical fits to the simulation data. As part of 
our QMC data fit function for two-leg ladders, we first 
obtain fits to the known dependence of the spin gap A 
on J and J' for isolated ladders. We then discuss the 
first few terms of the high-temperature series expansion 
(HTSE) for the magnetic susceptibility, which will also 
be incorporated into the fit function so that the function 
can be accurately extrapolated to arbitrarily high tem- 
peratures. The fit function itself is then presented and 
discussed. For some sets of exchange parameters (for the 
frustrated trellis layer) it was not possible to obtain ex- 
tensive QMC X* (t) data at low temperatures due to the 
"negative sign problem" . For such cases, it is sometimes 



A. Spin Gap for Isolated Two-Leg Ladders 



For any finite J'/J > an energy gap ("spin gap") A 
exists in the magnetic excitation spectrum between the 
singlet ground state and the lowest triplet excited states 
of 5 = 1/2 two-leg Heiaenberg ladders. □ The numerical 
A values of Barnes et aln in the range < J' /J < 1, ob- 
tained by extrapolating exact diagonalization results for 
short two-leg ladders to the bulk limit— were previously 
fitted by one of us with the expressions 
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S = 1/2 Two-Leg Ladder Spin Gap A 
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Fit: Greven A/J' Data 
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FIG. 20. Spin gap A of isolated two-leg ladders versus (a) 
J'/JJor < ,/'/J < 1 fronjWjiumerical calculations of Barnes 
et alB and of Greven et alH and (b) J/ J' % < J/ J' < 1 
from pumerical calculations of Dagotto et aln and of Greven 
et alH The dashed curve in (b) is the exact dimer series 
expansion to fourth order in J/ J' in Eq. (^. In (a) and (b), 
the solid curves are weighted fits to the data of Greven et al. 
by Eqs. (^) and (|), respectively. 



A* 



A 

7 



0.4 



0.1 



(6) 



Higher accuracy A values were subsequently obtained 
numerically from Moate Carlo simulations by Greven, 
Birgeneau and WiesetHi as shown in Fig. ^(a) where the 
previous data of Barnes et al. are shown for comparison. 
We find that the data of Greven et al. for J' /J < 1 in 
Fig. |2^(a) are fitted better by 



^ = 0.4030 



0.0989 ( y 



(7) 



with a statistical x^/DOF = 0.18 (for the definition, see 
the Appendix), as shown by the solid curve in the figure. 
For isotropic ladders (J'/J = 1), the fitted A/J = 0.5019 
is in good agreement with the values 0.50(1) of Barnes 
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FIG. 21. Dimer series expansions from 9th to 13th order 
in J/ J' of the spin gap A/J' of isolated two-leg ladders ver- 
sus J/J'i-feolid curves) as computed by Weihong, Kotov and 
Oitmaa.E23 The symbols identify the curves. 

et al.^ 0.504(7) of Oitmaa, Singh, and Weihong,EI 
0.504 of White, Noack and Scalapinol and 0.5028(8) of 
Weihong, Kotov and Oitmaa.E£il Using exact diagonal- 
izations for ladders of size up to 2 x 15 spins, Flocke 
obtained an extrapolated value 0.49 996 for |-tiM5 bulk 
limit, and conjectured that the exact value is i.ll23 Equa- 
tion (0) predicts values of Aq/ J which are systematically 
larger than the extrapolated bulk-limit values of Flocke 
at smaller J'/J, e.g. by 0.003 at J'/J = 0.8 and by 0.02 at 
J'/J = 0.2. The fitted initial slope iHpEq. (^ agrees with 
the estimates 0.41(1) of Grevea-|et aM and 0.405(15) of 
Weihong, Kotov and Oitmaa.P^ 



For the strong interchain coupling regime < J/ J' < 1 
of the two-leg ladder, the exact dimer secies-expaiision for 
the spin gap to seventh order in J/ J' isEHl22l'E2jil2a 
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The fourth-order series is plotted as the dashed curve 
in Fig. ^(b). Comparison of this prediction wth the 
numerical daia of Dagotto, Riera and ScalapinoO and of 
Greven et al£3 in the figure shows that the fourth-order 
series is a poor description of the data for J'/J > 0.5. 
The dimer series expansion has been qsHiputed to 13th 



Plots of their 
nl. The series 



order by Weihong, Kotov and Oitmaa.t22 
9th to 13th order series are shown in Fig 
is seen to converge very slowly with increasing order for 
J/ J' > 0.5. In fact, Piekarewicz and Shepard concluded, 
on the basis of dimer series expansions of the ground 
state energy per site at 50th order in perturbation theory 
for 4-, 6- and 8-rung ladders, that the radius of conver- 
gence of the dimer series is only « 0.7-0. 8. E£3 Therefore, 
for both of these reasons, to obtain an expression for 
A(J/J') to use in our QMC x*(t) data fit function for 
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the entire strong-coupUng range < J/ J' < 1, we-rar- 
ried out a weighted fit of the data of Greven et al£3 for 
0.1 < J/ J' < 1 in Fig. ^(b) by the simple two-parameter 
third-order polynomial 



Ao 
J' 



yielding the parameters 



a = 0.6878, 6 =-0.1861 



(9a) 



(9b) 



The first two terms in Eq. ( pa| ) were set to be the same 
as the corresponding exact terms in Eq. (||). The fit 
in Eqs. (^ is shown as the solid curve in Fig. ^(b). 
The high precision of the fit is characterized by the small 
xVDOF = 0.16. The value of A/J' at J/ J' = 1 is 
0.5017, which matches very well the value 0.5019 of the 
fit for < J'/J < 1 in Eq. (|^) at this isotropic-ladder 
crossover point between the two fits. 



where z is the coordination number of a spin by other 
spins and J is the unique exchange constant in the sys- 
tem. The same form of the HTSE of xiT)T is valid 
for any spin S, but where of course the numerical coef- 
ficients in Eq. (|l^) depend on S. Each term listed on 
the right-hand-side (but not the higher-order terms) de- 
pends only on z (and S) and not on any other feature 
of the spin lattice or magnetic behavior; as noted above, 
however, additional term(s) are added to the numera- 
tor of the last term if geometric frustration is present or 
if second-neighbor interactions are present (see below). 
Hence, one can generalize Eq. ( |To| ) to systems containing 
equivalent spins but unequal exchange constants Jij by 
the replacement zJ" J2j ^ij^ yielding using Eqs. (||) 



24i3 



8^2 

n -1 



(11) 



B. High- Temperature Series Expansions for 
and the Magnetic Specific Heat Cmag(T) 

As the second component of our fit function for the 
QMC x*{t) simulation data, we next consider the high 
temperature series expansion (HTSE) of x*(^) for a 
general Heisenberg spin lattice containing magnetically 
equivalent spins. Spins are magnetically equivalent if 
they have identical magnetic coordination spheres. Note 
that the HTSEs we discuss here, and HTSEs in general, 
are not restricted to AF couplings (with a positive sign as 
defined in this paper); the expansions are equally valid 
if the couplings are all FM (negative) or if they are a 
mixture of AF and FM couplings. 

HTSEs for x*(^) £^re calculated as, and the results are 
normally expressed directly as, a power series jin. 1/t. 



which we write as 



However, as mentioned by Rushbrooke and WoodE2ZI and 
discussed in Ref. the expressions for the expansion 
coefficients for a general spin lattice containing magneti- 
cally equivalent spins considerably simplify if the HTSE 
for 3x*t/[S'(S' -I- 1)] in powers of 1/t is inverted (the un- 
derlying pbjisics of this is unclear). Indeed, Rushbrooke 
and WoodE23 presented their calculated expansion coef- 
ficients in precisely this form. In fact for any Heisen- 
berg spin lattice (in any dimension) containing magneti- 
cally equivalent spins interacting with spatially isotropic 
nearest-neighbor AF or FM Heisenberg exchange, a sim- 
ple universal HTSE of x* {t)t exists up to second order in 



1/t, an 
order ,E 



for geometrically nonfrustrated lattices to third 
which we write for S = 1/2 as 



4fcBTx(r) 



zJ' 



zJ 
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(10) 



di do d-i 
1 + — + — + — 



with 



i j 



24J" 



(12a) 



(12b) 



(12c) 



Including only the first term on the right-hand-side of 
Eq. (12a) gives the Curie law x*{t) — C* /t with re- 
duced Curie constant C* = C/(iVg2^2 ) ^ g^g -t- 1)/3 = 
1/4, whereas the first and second terms together yield 
the Curie- Weiss law x*(*) = C*/{t - 9*) with re- 
duced Weiss temperature 6* = k ^O / J^^"^^ = — rfi = 
-(1/4) J-iilJ'^'""' (see Sec. PTd| below). 

A geometrically frustrated spin lattice is one in which 
there exist closed exchange path loops containing an odd 
number of bonds. Usually, the exchange path loops are 
triangles containing three bonds (such as in the 2D trian- 
gular lattice), where at least two nearest neighbors of a 
given spin are nearest neighbors of each other, although 
e.g. spin rings with any odd number of spins (and there- 
fore an odd number of bonds) are also geometrically frus- 
trated. Another example of a system containing triangu- 
lar exchange path loops is the dimer system SrCu2 (603)2 
(intradimer interaction J\ = J™^'^) with a partially frus- 
trating interdimer interaction J2 . One can show— tbat 
Eq. (y) agrees exactly to C'(l/t2) with the HTSeN-ES 
for x{T)T of this system. The frustration first becomes 
apparent in the HTSE as an additional additive term 
[(— 15/4)( J2/ Ji)2 in this case] in the numerator of the 
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coefBcient in the square brackets in Eq. (|Tl|). In the 
context of the present discussion, a second- or further- 
nearest-neighbor interaction is equivalent to a nearest- 
neighbor one in a system with geometric frustration, and 
hence the general expansion (O) for x*t is still exact to 
0{\/t^) for such systems, provided again that all spins 
are magnetically equivalent. 

For our isolated and coupled ladder QMC simulation 
fits, the three c?„ HTSE coefficients in Eqs. ( p^ ) are 



4xr 



4J" 



8J" 
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2jdiag 
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2 J'" 4 
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2jdiag 


%2J 




2J'"' 


+ 2 J3D'" 


2 J"^ - 


f 2 J'"' 




2J3D' 


- 9 JJ"' 



(13) 



where the last term in c?3, given in the HTSE in Ref. |2^ 
for x*{i) of the trelhs layer, arises due to the geomet- 
ric frustration in the trellis layer interladder coupling. 
The dn in Eqs. (|l3|) are the correct HTSE coefficients in 
Eq. (12a), except for da in the case of diagonal second- 
neighbor intraladder couplings J'^^'^s- this latter case 
we will not use da in the fit function. 

Weihong, Singh and Oitmaa have computed the HTSE 
for the product x(T)r of the isolated S* = 1/2 two- 
leg Heisenberg ladder with spatially anisotropic exchange 
in the leg and rung to 9th order in 1/T, which con- 
tains a total of 54 nonzero coefficients in powers of J 
and/or J'.t2a As anticipated above, the series simplifies 
if it is inverted. In addition, this inversion allowed us 
to easily estimate the rational fractions approximated 
by the ten-significant-figure decimal coefficients given by 
Weihong et al. Our result for the inverted ninth-order 
series, containing 42 nonzero terms, is 



1 + (2 J + J') I + (2 J2 + J' 



2\ X_ 
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116J^ -I- 99JJ'* - 32J"'^) ^ 
317J^ - lllJ'^J'^ - 96J^ J'^ 



^ 6427^'^^ + 297 J J'^ - 32 J'^^ ^ 



+ (792 + 3444JV^ -I- 5068 

- 22932 J'" - 10332JV^ 
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10395 JJ'^ + 2h%.f^ 



1440 

4 T/3 



40320 



6165J* - 411JV^ -f 604JV^ 



264777"* J'" - 8220 JV^ 
9580 JV^ - 9702 J J''' + 64J'^^ ^ 



40320 

(^23674J^ - 29916 - 46269J^ J'^ 

228168 jV' - 65340 J" J'^ 
78516JV^ - 51840JV'^ 



where 



- 68607 J J' + 128 J'' 



1 



9\ x-- 



362880 



(14) 
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Upon inverting the series in Eq. (14) and then con- 
verting each resulting rational fraction coefficient to a 
ten-significant-figure decimal value to compare with the 
HTSE of Weihong et ai, each of the 54 coefficients is 
found to be identical to the corresponding ten-significant- 
figure coefficient given by Weihong et alxEB The HTSE 
in Eq. ([l|) is identical to order l/T^ with the HTSE for 
(xT)~^ in Eq. (12a) in which the di, c?2 and ds coeffi- 
cients are given for the general S" = 1/2 two- leg ladder 
by Eq. (|3|), but where in the present case only J and 
J' are nonzero. An interesting aspect of the HTSE in 
Eq. (|lj) is that in the expression for the coefficient of 
each term shown, the coefficient of the J"~* J' term 
vanishes. 

Gu, Yu and Shcn have derived an analytic expression 
for the magnetic field- and temperature-dependent free 
energy of the two-leg ladder for strong interchain cou- 
plings J/ J' sS^l using perturbation theory to third or- 
der in J/J'!Ei Our HTSE of l/[4:X*{t)t] obtained from 
their free energy expression is identical to order l/T^ 
with Eq. ([1^). As expected, the coefficients of the fourth 
order and higher order terms of the HTSE do not agree 
with the corresponding correct coefficients in Eq. (H). 

Just as there is a universal expression for the first three 
to four HTSE terms for x{T)T of a Heisenberg spin lat- 
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tice containing magnetically equivalent spins as discussed 
above, a universal HTSE for the magnetic specific heat 
Cmag{T) of such a spin lattice exists to order to 
and is given for 5 = 1/2 by 



Nkn 



3 

32 



(15) 



The sums are over all exchange bonds from any given spin 
Si to magnetic nearest- neighbor spins Sj . The first term 
holds for any spin lattice containing magnetically equiv- 
alent spins, but the second term holds only for geomet- 
rically nonfrustrated spin lattices in which the crystallo- 
graphic and magnetic nearest-neighbors of any given spin 
are the same. Higher order terms all depend on the struc- 
ture and dimensionality of the spin lattice. The HTSE 
for Cinag(T') to (lowest) order 1/T^ is the specific heat 
analogue of the Curie- Weiss law for the magnetic sus- 
ceptibility, i.e., they can both be derived from the same 
lowest (first) order term in l/T of the magnetic-nearest-i 
neighbor instantaneous two-spin correlation function.^ 
Equation (^5|) is therefore accurate in the same high- 
temperature region in which the Curie- Weiss law for the 
magnetic susceptibility is accurate. Physically, the rea- 
son that the lowest-order HTSE terms of Cniag{T) are 
of the form Jfj/T"^ is that Cmag(r) cannot be negative, 
regardless of the sign(s) of the Jij. 

For the two-leg spin ladder couplings considered in this 
paper, to lowest order in l/T Eq. ( p^ ) yields 

CraasiT) _ 3 ^2 + j' ^ 2 J^^'^S^ 3 j"' 



iVfcB 32(fcBr)2 



2 J'"' + 2 J^D'") . (16) 



From comparison with Eq. (|T|), the HTSE for Cmag(r) 
of isolated S* = 1/2 spatially anistropic two- leg Heisen- 
berg ladders (J, J' 7^ 0) given to lowest order (1/T^) in 
Ref. 157 is found to be incorrect. 



C. General x*(*) Fit Function 

The following fit function incorporating the above 
considerations, and containing the Fade approximant 
V^^^it), was found capable of fitting the QMC x*(i) data 
for a given exchange parameter set to within the accuracy 
of those data (i.e., to within a x^/DOF ^ 1) 



At 



1 + 
1 + 



(17a) 



(17b) 



which satisfies the Curie law at high temperatures, where 
A^j. is not necessarily the same as the true spin gap. In 
order to further constrain the fit and also to produce a 
fit which can be accurately extrapolated to high temper- 
atures, we require that a HTSE of Eqs. ( [l^ ) reproduce 
the HTSE in Eqs. (|ll])-(|3|), which in turn yields the 
constraints 
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D3 ^ {ds + d2Ni + diN2 + N3) - Afit(d2 -I- diNi + N2) 

(18c) 



Al^ AL^ 
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In general, one has 



n—p 



m— 



(19) 



Unless otherwise explicitly noted for a specific fit, 
Di, D2 and D3 are not independent fitting parameters 
but are rather determined from the fitting parameters 
iVi, N2, N3 and A* J in Eqs. (|l|), where A^^ can also 
be a fitting parameter. To obtain a fit to a QMC x*{t) 
data set for a specific set of exchange constants to within 
the accuracy of the data, i.e. which yielded x^/DOF ~ 1, 
typically required a total of 6-9 independent fitting pa- 
rameters, which was essentially independent of the num- 
ber of data points in the data set. 

Finally, we reformulated the fit function into a two- 
or three-dimensional one so that it could not only in- 
terpolate and extrapolate x* versus t for a given set of 
exchange constants but could also interpolate x*(0 for 
a range of exchange constants. To do this, we expressed 
the parameters iV„, D„ and sometimes A^^ in Eqs. (|T^ ) 
as power series in the exchange constants; this also con- 
siderably reduced the total number of fitting parameters 
required to obtain a global fit to x*(0 data for a given 
range of exchange constants. This scheme was success- 
fully used except for exchange constant ranges traversing 
a QCP, for which two piecewise continuous interpolation 
fits were required for the two exchange constant ranges 
on opposite sides of the QCP, respectively. The result- 
ing fits to the QMC x*(0 simulation data and several 
exchange parameter interpolations are shown as the sets 
of solid curves in the above QMC data figures, as de- 
scribed in the captions. For most of the QMC simula- 
tions, a x^/DOF ^ 1 was obtained. The high quality of 
the fits may therefore perhaps be appreciated from the 
small errors estimated for the QMC data, especially at 
the higher temperatures, which varied from 1-10% for 
0.01 <t<0.1 to 0.03-0.1% for t > 0.5. The details of 
the fitting procedures and tables of fitted parameters are 
given in the Appendix. 
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D. Fit Functions for X*(*) Derived from 
Molecular Field Theory 

For Heisenberg spin lattices consisting of identical spin 
subsystems which are weakly coupled to each other, it is 
sometimes necessary to use a fit function for theoretical 
X*{t) data in the paramagnetic phase which contains a 
minimum number (perhaps only zero, one or two) of fit- 
ting parameters, and which still provides a reasonably 
good fit to the data. Such fit functions can be provided 
by molecular field theory (MFT) and its extensions as 
described in this section. Each isolated spin subsystem 
is assumed to have a known susceptibility Xo(^)- It can 
be easily shown that if each spin in the entire system is 
magnetically equivalent to every other spin, with spins 
in adjacent subsystems coupled by Heisenberg exchange, 
then the reduced susceptibility x* (t) in the paramagnetic 
state of the system is given by MFT as 



or equivalently 



1 1 



+ A , 



x*it) xm 

where the MFT coupling constant A is given by 



A 



/ J Jn\a.yi 
j 



(20a) 



(20b) 



(20c) 



the prime on the sum over j signifies that the sum is only 
taken over exchange bonds from a given spin Si to 
spins Sj not in the same spin subsystem, and J™^'^ is the 
exchange constant in the system with the largest magni- 
tude. By definition, the expression for Xoi^) does not 
contain any of these Jij interactions which are external 
to a subsystem. Within MFT, Eqs. (^ ) are correct at 
each temperature in the paramagnetic state not only for 
bipartite AF spin systems, but also for any system con- 
taining subsystems coupled together by any set of FM 
and/or AF Heisenberg exchange interactions. The only 
restriction, as noted above, is that each spin in the system 
is magnetically equivalent to every other spin in the sys- 
tem. Thus, our fit functions derived in this section could 
have been used to fit x*{t) data for any of the coupled 
two-leg ladder spin lattices discussed in this paper, al- 
though in general to much lower accuracy than obtained 
in the above section and the Appendix. However, they 
do not apply, e.g., to trellis layers containing three- leg 
ladders, because in this case the spins are not all mag- 
netically equivalent since the magnetic environment of a 
spin in the central leg of such a ladder is different from 
that of a spin in the outer two legs of the ladder. 

Before proceeding further, we first make contact with 
the familiar case in which a subsystem consists of a single 
spin. Then XoC?") is the Curie law, 



Xa{T) = ^ , 

where the Curie constant is 

Ng^^JilS{S+l) 
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In reduced units, the Curie law reads 
_ Xo(r)J™^'' _ C* 



xm 



Ng^ii t 



(21a) 



(21b) 
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where the reduced Curie constant C* and reduced tem- 
perature t are defined as 



C 5(5 + 1) 



Ng^ixi 3 
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Then our general expression ( 20a ) incorporating inter- 
actions between the spins yields the Curie-Weiss law 
X{T) = C/{T - 9) in reduced form as 



x*it) 



C* 
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with reduced Weiss temperature 6* given by 
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where we have removed the prime from the sum because 
in this case all exchange interactions in the system are 
external to a subsystem which consists only of a single 
spin Si. 



Equation (20a) can be used as a fit function containing 
no adjustable parameters to parametrize numerical x*{i) 
data for spin systems with weak intersubsystem coupling. 
In the following, we extend the MFT framework to pro- 
vide latitude for including adjustable fitting parameters 
to improve the quality of the fit. To obtain a general form 
for the fit function for 5 = 1/2 Heisenberg spin systems 
we first rewrite the HTSE for x*(i) in Eqs. (p^), absorb- 
ing the HTSE terms for a subsystem back into the exact 
Xo {t) for the subsystem which already implicitly contains 
the correct HTSE for the subsystem, leaving only the ex- 
ternal interactions explicit, yielding the modified HTSE 



1 



1 



r{t) xm t ^ 



(23a) 



with 
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44 = - — y , 44 = — 2 y Ji, , (23b) 
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(23c) 



where the prime on the sums has the same meaning as 
in Eq. ( ^Oc| ). Again, dg has additional terms if geomet- 
ric frustration and/or second-neighbor interactions are 
present in the intcrsubsystem couphngs. Comparison of 
Eqs. (11) and (|20|) shows exphcitly that MFT exactly 
yields the first expansion term (A — 4:di) of the quan- 
tum mechanical HTSE for x*(t) in terms of the exchange 
constants external to a subsystem. Note that the Weiss 
temp erat ure in the Curie- Weiss law is always given by 
Eq. ( ^2b| ), where the sum over j is over all magnetic 
nearest neighbors of a given spin and not just over those 
external to a subsystem. 

We now rewrite the HTSE in Eqs. (ESl) in the form 



with 



Ad' 



ArV 

f{J,^,t)^4d[ + -^^^ 



(24a) 



(24b) 



where the Ad'^ parameters for n — 1-3 are the same as 
given in Eqs. (23b) and (23c). Equation (|24a|), which 



is an extension of the MFT prediction in Eq. (20a), can 
be used as a function to fit x*(t) data for coupled two- 
leg ladders, where Xo(^) then the susceptibility of iso- 
lated ladders. The function /(Jy,t) in Eq. ( |24bD , which 
contains (apart from J™^'^) only the intcrsubsystem ex- 
change constants coupling the ladders to each other 
and is expected to be most accurate at high tempera- 
tures, can be modified to provide for the introduction of 
adjustable fitting parameters as will be further discussed 
in the Appendix where we obtain fit functions for our 



QMC x*{t) data for the two-leg ladder trellis layer. In 
addition, especially when the intcrsubsystem interactions 
significantly change the spin gap, terms which include the 
Jij interactions external to a subsystem and additional 
fitting parameters could be included in the Xo(^) function 
itself. 



IV. DISPERSION RELATIONS 

The one- and two-magnon dispersion relations E{ky) 
were computed for S" = ^ 2 x 12 ladders by exact di- 
agonalization using the Lanczos algorithm for J' / J = 
0.5, 0.6, 1.0 and in each case for wavevectors 

kya/iT = 0, 1/6, 2/6, . . . , 1, where ky is the wavevector 
in the ladder leg direction and a is the nearest-neighbor 
spin-spin distance; in the discussion below we set a — \. 
The results are given in Table |[ Our value of the spin- 
gap for J'/ J = 1, A/ J = 0.514 784, is identical to the 
six significant figures with that calculated for the same 



spin lattice by Yang and HaxtonHHj This value is about 
1.8% higher than for the 2 x 16 ladder (0.505460 384, 
Ref. 12£) and for the bulk hmit discussed in Sec. III. 
The one-magnon dispersion relation data are shown as 
the symbols in Fig. |2^. In the limit J' / J —>■ 0, the exact 
dispersion relation calculated for the 5=1/2 AF uniform 
Heisenberg chain by des Cloizeaux and Pearson in 1962-is. 
E{ky) = (^J/2)|sin(^)| = (7rJ/2)[i - i cos(2fc^lli/2>02 
also shown in Fig 



Eiky 



J' 



2 2 

2L wjlp^cas for J' /J ^ 1 one has 
Jcos(fcy).Elll23 Since our dispersion rela- 



tions are for exchange constant ratios 0.5 < J'/ J < 1 
closer to the former limit, we obtained exact fits to the 
data bM,4h|e square root of an even seven-term Fourier 
series jHa't£a shown as the solid curves in Fig. |2^. The 
spin-gap A for each J'/ J occurs at wavevector ky — n/a. 



TABLE I. One-magnon dispersion relations {k^a/Ti = 1) and the lower boundary of the two-magnon continua {kxa/-n — 0) for 
S = 1/2 2 X 12 two-leg Heisenberg ladders with J'/ J = 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0, calculated by exact diagonalization using 
the Lanczos algorithm. The spin-gap between the singlet ground state and the lowest triplet excited state is A = E{kya/n = 1) 
for the one-magnon dispersion relation. 

kya/n 








1/6 


1/3 


1/2 


2/3 


5/6 


1 


J'/J 








One-Magnon E/J 








0.5 


1.13926 


1.25058 


1.62977 


1.79736 


1.52857 


0.908547 


0.289208 


0.6 


1.28127 


1.38810 


1.70428 


1.82438 


1.52561 


0.902179 


0.319648 


0.7 


1.42954 


1.52908 


1.78656 


1.85470 


1.52828 


0.904685 


0.358123 


0.8 


1.58158 


1.67096 


1.87366 


1.88871 


1.53729 


0.915965 


0.403962 


0.9 


1.73510 


1.81177 


1.96301 


1.92667 


1.55290 


0.935748 


0.456448 


1.0 


1.88782 


1.94993 


2.05281 


1.96863 


1.57507 


0.963567 


0.514784 


J'/J 








Two-Magnon E/J 








0.5 




0.960355 


1.56230 


1.86742 


1.78996 


1.43543 


1.15334 


0.6 




1.01850 


1.61171 


1.92615 


1.87887 


1.57866 


1.35375 


0.7 




1.08763 


1.67135 


1.99573 


1.98118 


1.73555 


1.55798 


0.8 




1.16790 


1.74149 


2.07589 


2.09499 


1.90220 


1.76404 


0.9 




1.25911 


1.82206 


2.16602 


2.21834 


2.07507 


1.96985 


1.0 




1.36070 


1.91271 


2.26527 


2.34932 


2.25116 


2.17341 
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Spin S = 1/2 2-Leg Ladder (2 x 
One-Magnon Dispersion Relati 
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y 

FIG. 22. One-magnon dispersion relations for the isolated 
two-leg 2 X 12 ladder calculated for J'/J = 0.5, 0.6, . . . , 1.0 
using the Lanczos algorithm. The wavevector ky is parallel to 
the legs of the ladder and a is the nearest-neighbor spin-spin 
distance {kxa/n — 1). The data for J'/J = 2 calculated 
for the-same lattice using the same algorithm by Barnes asd 
Rierat£j and the dispersion relation for the isolated chainEi!3 
are shown for comparison. The solid curves for J'/J > are 
exact fits to the data by the square root of an even seven-term 
Fourier series with period Akya/n — 2. 

Also include/d in Fig. |22| are the earlier results of Barnes 
and Rierali£3 for J'/J = 2 computed using the same algo- 
rithm on the same spin lattice. Their data for J'/ J = 1 
and 0.5 (not shown) are in good agreement with our data 
for these J'/J values. 



A notable feature of the data in Fig. 22 is that the ra- 
tio E™-'^^/ A of the maximum to the minimum energy of 
each dispersion relation estimated using the above exact 
fits to the data is a strong function of J'/J, as shown 
in Fig. 23, which facilitates obtaining highly precise esti- 
mates of J'/J from neutron scattering data. This obser- 
vation, previously made by us on the basis-jof the earlier 
dispersion relations of Barnes and Rieraji£3 was used by 
Eccleston et al. to estimate J'/J for the two-leg ladders 
in Sri4Cu2404i from their neutxan scattering data on 



single crystals of this compound.Ei_l Their data, in turn, 
motivated us in the present work to compute the disper- 
sion relations on a finer grid of J'/J values than existed 
previously. We will discuss t hese calculations and exper- 
imental data further in Sec. 



VIII 



The lower boundary of the two-magnon continuum 
{kx = 0) is shown in Fig. |2^ for J'/J < 1 over much 
of the ky range. From our data on the finite-size ladder, 
we cannot clearly distinguish between the two-magnon 
scattering states and bound states lying near the lower 




0.4 0.5 0.6 0.7 0.8 0.9 1.0 
J'/J 

FIG. 23. Ratio of the maximum to the minimum 
one-magnon excitation energy, i?™™/A (•), versus J'/J for 
the 2 X 12 two-leg ladder, for magnons with wavevector par- 
allel to the legs of the ladder, from the fits to the data in 
Fig. p3. The solid curve is a polynomial fit to the data. 
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FIG. 24. The lower boundary of the two-magnon contin- 
uum {kxa/-n = 0, open symbols) for the isolated two-leg 
2x12 ladder calculated for J' / J = 0.5, 0.6, . . . , 1.0 using the 
Lanczos algorithm. The wavevector ky is parallel to the legs 
of the ladder and a is the nearest-neighbor spin-spin distance. 
For comparison, the dispersion relations for one-magnon ex- 
citations {kxa/-K = 0, filled symbols, from Fig. ^2|) are in- 
cluded. The dashed curves for J'/J > are exact fits to 
the two-magnon data by the square root of an even six-term 
Fourier series with period Akyu/ir — 2. 



boi 



203 



ary of the two-magnon continuum for ky 



Each dashed curve is an exact fit to the respec- 
tive two-magnon data by the square root (see above) of 
a six-term Fourier series. The two-magnon excitations 
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are degenerate with the one-magnon spectra over much 
of the Brillouin zone. 

Interladder coupHng within the trelhs-layer (J") has 
virtuaUy no influence on the one-magnon dispersion and 
on the spin structure factor close to the dispersion mini- 
mum, as the contributions due to this couphng interfere 
destructively at kya = tt. A coupling J'" in the third (z) 
dimension, perpendicular to the trellis layer, will how- 
ever give an additional dispersion J'" cos{kz) that has 
to be taken into account. Away from the minimum the 
magnon band disperses also due to the trellis layer inter- 
ladder coupling J", most strongly close to the dispersion 
maximum apojand fcj,a ~ 7r/2. As was shown by Lidsky 



and Troyer,li^the band center is not moved substantially 



by J", and averaging over all momenta perpendicular to 



the ladders, as done by Eccleston et al. 
recovers the uncoupled ladder dispersion 



12: 



essentially 



The neutron 

scattering function depends of course on the relative con- 
tributions of all the magnetic excitations. Our calcula- 
tions of the dynamical spin structure factor S{q^ lo) for 



a 2 X 12 ladder and for the experimentally relevant (see 
Sec. Vn) intraladder coupling ratio J'/J — 0.5, shown 
in Fig. ^5|, demonstrate that around the top of the one- 
magnon band the two-magnon states have energy and 
weight comparable to those of the single magnon band. 
Therefore the two-magnon states should be taken into 
account when fitting inelastic neutron scattering data to 
obtain the part of the one-magnon dispersion relation at 
the higher energies. 

Related results for S{q,uj) oi S — 1/2 two-leg Heisen- 
berg ladders have been obtained previously. A study of 
the 2 X 12 ladder with spatially anisotropic exchange by 
Yang and Haxton indicated that the contribution of the 
lowest one-magnon triplet states to the response func- 
tion at wavevector {tt,tt) increased from 91.7% to 96.7% 
of the total response at that wavevector as J'/J increased 
from 0.4 to 1; the totaLjjesponse itself at this wavevector 
peaked at J'/J ~ 0.5.c£d Calculations of the odd number 
of magnons sector of S{q, uj) for the 2x16 ladder with 
J'/J = 1 were reported by Dagotto et al&3B Their 




1 2 3 1 2 3 



co/J ro/J 

FIG. 25. Dynamical spin structure factor S{qy,ij) (in arbitrary units) for the isolated S — 1/2 AF two-leg Heisenberg ladder 
with L = 12 rungs for the physically relevant coupling J'/J =1/2 and for both the odd number of magnons (left panel) and 
even number of magnons (right panel) sectors. As in the text, the y direction is defined to be along the ladder legs and the x 
direction along the ladder rungs. In the abscissa label we have set h = 1. Note that the ordinate scale is different for the odd 
number of magnons sector at qy = 57r/6 and qy = n than for the other plots. It is apparent that for momenta along the ladder 
qy near the top of the one-magnon band, the two-magnon continuum (lowest peaks in the even sector) has energy and spectral 
weight comparable to those of the single magnon collective mode. 
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results showed that the one-magnon contribution to 
S'(q, w) continuously decreases as the wavevector de- 
creases from (tt, tt) to (7r,0),llHa which is qualitatively the 
same as we have found for J' / J — 0.5. 



LDA+U CALCULATIONS OF EXCHANGE 
CONSTANTS IN SrCuaOs 



An ab-initio calculation using the LDA+U methodcB: 
was enlisted to compute the electronic structure of 
SrCu203 and to extract from it the exchange cou- 
plings. The atomic coordinates used are those given in 
Sec. VI below. The LDA-I-U method has been shown 



to give good results for insulating—tcansition metal ox- 

The exchange in- 



20£ 



ides with a partially filled d-shell. 
teraction parameters can be calculated using a proce- 
dure based on the Greens func.tiO|n rnethod which was 
developed by A. I. Lichtenstein.c23E23 This method was 
successfully applied for calcu lation of the exchange cou- 
plings in KGwFs (Ref. 206) and in layered vanadates 
CaV„02„+i.N 

The LDA-fU methodEj'M is essentially the Local 
Density Approximation (LDA) modified by a potential 
correction restoring a proper description of the Coulomb 
interaction between localized d-electrons of transition 
metal ions. This is written in the form of a projection 
operator 



LSDA 



E 



inlma)V^^. 



I {inlr 



(25) 



' = yZ{{m,m" 



mm' / J 
{m} 



Vee I m\m"')n- 



(26) 



+ {{m,m" I Vee I m',m"') 
-{m,m" I Vee I m"',m')X, 

-L/(A.-i)+,/(A.^-l) 



where | inlma) {i denotes the site, n the main quantum 
number, I- orbital quantum number, m- magnetic num- 
ber and a- spin index) are c?-orbitals of transition metal 
ions. The density matrix is defined by 



1 



ImG: 



inim.tnim' 



(E) dE 



(27) 



where G" , . , ,(E) 



[inlma \ {E — H) ^ \ inlm a) 
are the elements of the Green function matrix, N°' = 
Tr{n^^,), and A'^ = A^^ + A^^. U and J are screened 
Coulomb and exchange paraipeters calculated via the so- 
called "supercell" procedurec23 and found to he U — 
7.79 eV and J = 0.92 eV, respectively. The calcula- 
tion scheme was realized in the framework of the Lin- 
ear Muffin-Tin Orbitals (LMTO) methodEB based on the 
Stuttgart TBLMTO-47 computer code. 



The inter-site exchange couplings were calculated with 
a formula which was derived using the Green function 
method as the second derivative of the ground state 
enerey— spt|i respect to the magnetic moment rotation 
anglec2Z 



/ J mm -^mm' m" m'" m"m"' 
{m} 



(28) 



where the spin-dependent potentials / are expressed in 
terms of the potentials of Eq. (Effl as 



■n 

mm' 



min' 



(29) 



The effective inter-sublattice susceptibilities are defined 
in terms of the LDA-I-U eigenfunctions -0 as 

?^nkT - '^n'kj l ilm* i jlm" ,ilm' J.jl"^"" 



mm m m 



E 



(30) 

Equation (^) was derived as a second derivative of 
the total energy with respect to the angle between spin 
directions of the LDA+U solution. The LDA+U method 
is the analogue of the Hartree-Fock (HF, mcaripEeld) ap- 
proximation for a degenerate Hubbard model.t£3 While 
in the multi-orbital case a mean-field approximation gives 
reasonably good estimates for the total energy, for the 
non-degenerate Hubbard model it is known to underes- 
timate the triplet-singlet energy difference (and thus the 
value of effective exchange parameter Jj^) by a factor of 

two for a two-site problem (i?HF = ^ and inexact = , 
where t <^ U is the inter-site hopping parameter) . Thus 
the J value calculated by expression ( |28| ) was multiplied 
by a factor of two to correct the Hartree-Fock value. The 
calculated results are presented in Table |l[ 



TABLE II. Exchange constants in SrCu2 03, calculated 
using the LDA+U method, along wit h j t^oe e obtained for 
CaV2 05 and MgV2 05 by Korotin et alJ^it^ using the same 
method. For the definitions of the exchange constants, see 
Eq. (^. Antiferromagnetic (AF) and ferromagnetic (FM) ex- 
change constants are positive and negative, respectively. J™^" 
is the largest (AF) exchange constant in the system. 





SrCuaOs 


CaVaOs 


MgVaOs 


J/fcB (K) 


1795 


122 


144 


J'/kB (K) 


809 


608 


92 


J''"'<'/kB (K) 


-200 


20 


19 


J"/kB (K) 
J"'/kB (K) 


4 
3 


-28 


60 


jrmax 


J 


J' 


J 




1 


0.201 


1 


j/yjmax 
^diag 1 jmax. 


0.451 


1 


0.64 


-0.111 


0.033 


0.13 




0.002 


-0.046 


0.42 




0.002 
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VI. STRUCTURE REFINEMENTS OF SrCuaOg 
AND Sr2Cu3 05 



VII. EXPERIMENTAL MAGNETIC 
SUSCEPTIBILIES AND MODELING 



Up to now the detailed crystal structure of the two- 
leg ladder compound SrCu203 has not been reported. 
This is necessary as input to our LDA+U calculations 
in the previous section. Here we give our refinement re- 
sults for this compound as well as those for the three-leg 
ladder compound Sr2Cu305. Powder X-ray diffraction 
data were collected with a Rigaku RAD-C diffractometer 
equipped with a graphite crystal monochrometer (CuKa 
radiation, 30 kV, 100 mA). Data were collected from 20 to 
120 ° with a step width of 0.02 °. Lattice parameters and 
the atomic positions wece refined by a Rietveld technique 
using Rietan software.EH 



Our crystal data for SrCu203 and Sr2Cu305 are given 
in Tabic III. The atomic positions for Sr2Cu3CXa are 
in good agreement with those of Kazakov et aZH de- 
termined from Rietveld refinement of powder neutron 
diffraction data, but our lattice parameters are signifi- 
cantly smaller, indicating a possible homogeneity range 
and/or variable density of defects in this compound. Se- 
lected interatomic distances and bond angles are listed 



for both compounds in Table IV 



A. Introduction 

For the three cuprate spin ladder compounds for which 
experimental magnetic susceptibility x{T) = M/H data 
are presented below, where M is the magnetization and 
H = 0.1 or IT is the applied magnetic field, we fitted 
the x{T) data per mole of Cu by the expression 



X{T) = Xo + 



-'imp 

T -e 



where 



Xo = X 



(31a) 



(31b) 



X^P'"(T) 



^max 

= (0.3751 



x*{t) 

cm^ K 
mol I J'^^'^/ks 



xi^) (31c) 



and Na is Avogadro's num ber. The first term xo hi 
Eq. (31a), according to Eq. (31b), is the sum of the or- 
bital diamagnetic core contribution x'^"'^'' ^^^d paramag- 
netic Van Vleck contribution x^^ which are normally 
essentially independent of T; these contributions, espe- 
cially x^^j are very difficult t o est imate accurately a pri- 
ori. The second term in Eq. ( |31a| ) is the extrinsic impu- 
rity Curie- Weiss term with impurity Curie constant Cimp 
and Weiss temperature 9 which gives a low-temperature 



TABLE III. Crystal structure refinement data for SrCu203 
and Sr2Cu305. Sample color: black; radiation: CuKq; tem- 
perature: 290 K; monochromator: graphite; 29 range: 20 
to 120°; step width: 0.02°; crystal system: orthorhombic; 
space group: Cmmm. SrCu2 03: lattice parameters (A): a — 
3.9300(1), h = 11.561(4), c = 3.4925(1); density: 5.497 g/cm^ 
= 4.51%; Rp = 2.78%; i?e = 2.76%; Ri = 2.76%; 
Rf = 3.97%. Sr2Cu305: lattice parameters (A): a — 
3.9292(1), b = 19.396(5), c = 3.4600(1); density: 5.614g/cm^ 



Rwp 


= 5.03%; Rp 


= 3.35%; i?. 


= 3.00%; Ri = 


6.92%; 


Rf = 


6.28%. 








atom 


site 


X 


y 


2 


SrCu2 03 


Sr 


2d 








1/2 


Cu 


4i 





0.3348(5) 





0(1) 
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1/2 








0(2) 


4i 





0.1728(24) 





Sr2Cu305 


Sr 


4i 





0.3987(5) 


1/2 


Cu(l) 


2a 











Cu(2) 


4:i 





0.1996(5) 





0(1) 


2b 


1/2 








0(2) 


4i 





0.1028(24) 





0(3) 


4i 





0.3001(24) 






TABLE IV. Selected bond lengths (A) and angles (°) for 
SrCu203 and Sr2Cu305. 



SrCu2 03 



Sr-O(l) 




2.629(0) 


Sr-0(2) 




2.654(20) 


Cu-O(l) 


rung 


1.910(6) 


Cu-0(2) 


leg 


1.967(13) 


Cu-0(1)-Cu 


rung 


180.00 


Cu-0(2)-Cu 


leg 


174.22 


Cu-0(2)-Cu 


interladder 


92.56 


Sr2Cu305 


Sr-O(l) 




2.617(7) 


Sr-0(2) 




2.618(1) 


Sr-0(3) 




2.580(36) 


Cu(l)-0(1) 


center leg 


1.965(1) 


Cu(l)-0(2) 


rung 


1.995(55) 


Cu(2)-0(2) 


rung 


1.877(53) 


Cu(2)-0(3) 


side leg 


1.965(0) 


Cu(2)-0(3) 


interladder 


1.949(2) 


Cu(l)-0(1)-Cu(l) 


center leg 


180.00 


Cu(l)-0(2)-Cu(2) 


rung 


180.00 


Cu(2)-0(3)-Cu(2) 


side leg 


180.00 


Cu(2)-0(3)-Cu(2) 


interladder 


90.17 



26 



upturn in x(T) not predicted by theory for the intrinsic 
spin susceptibiHty x^'""(^) and is assumed to arise from 
paramagnetic impurities and/or defects. The fitted Cimp 
values are typically ^ 10~^ cm^ K/mol Cu, corresponding 
to a few tenths of an atomic percent with respect to Cu of 
paramagnetic species with 5 = 1/2 and g = 2. The 9 is 
typically ~ —2 K which may indicate AF interactions be- 
tween the impurity magnetic moments, the occurrence of 
single-impurity-ion crystal field effects, and/or compen- 
sate for paramagnetic saturation of the magnetic impu- 
rities at low temperatures in the fixed field of the mea- 
surements. For a given sample Cimp and 9 are nearly 
independent of the model and parameters for x''P'"(T), 
for which the QMC data were presented and fitted in 
previous sections. 

For a given experimental x{T) data set, there are po- 
tentially at least six fitting parameters: xqj Cmp, 9, g, 
jm&x least one add itional exchange parameter. 

In addition, from Eq. (31c), the fitted g and J™^'^ are 
intrinsically strongly correlated and thus even minor in- 
accuracies in the experimental data can cause the fitted 
g and J™'*'^ parameters for a given type of fit to vary 
significantly (±20% or more) from sample to sample of 
a given compound. Therefore it is important to con- 
strain j-the g-value to lie within a physically reasonable 
range .c3 Unfortunately, ESR measurements of the in- 
trinsic (bulk) Cu+^ g-values in SrCu203, Sr2Cu305 and 
LaCu02.5 are not available, although g — 2.14 was 
ported for Cu defects in both SrCu2 03 and Sr2Cu305.Ei2 
Fortunately, in all Cu+'^-containing oxide compounds for 
which ESR data are available of which we are aware, the 
powder- aver age (7- value is always in the narrow approxi- 
mate range 1 n(.'V^ as JIln^sitratpH for representative com- 



pounds in Table ^.L 



We will also be inves- 
tigating the xiTfof CaV2 05 and MgV2 05 below. The 
^-values for 5 = 1/2 compounds are observed to be 
in a narrow range about g = 1.96, the sign and magni- 
tude of the deviation from 2 being respectively opposite 
and smaller than for Cu+^ due to the opposite (positive) 
sign and smaller magnitude of the spin-orbit coupling 
constant for V compared to that of Cu. Representative 
5'- values nVigpr-irprl frir V+4 ty^pplpq in sBveral materials are 
given in Table 0|llj>E3EI3&Ea In our fits, we will use 
the fixed value g — 1.96 determined for powder samples 
of CaiiaQ^and MgV205 by Onoda and coworkers using 
ESRi 



144,150 



Before proceeding to presentation and modeling of the 
experimental x(T) data, we comment briefly on units. 
The popular commercial Quantum Design SQUID mag- 
netometer, used also here, reads out the magnetic mo- 
ment of a sample in units of "emu" . Unfortunately 
this "unit" is useless for unit conversions, and authors 
use this same "unit" variously for magnetic moment 
and magnetic susceptibility, which of course are not the 
same quantities. Here we use cgs units with one excep- 
tion (T). The unit for magnetic moment (an "emu") is 
Gcm3 = erg/G (e.g., 1 = 9.274 x IQ-^iGcm^), for 
molar magnetization Gcm'^/mol, for magnetic field H 



and magnetic induction B G = Oe, for susceptibility 
of a sample cm'^ and for molar susceptibility cm'^/mol. 
For convenience, we occasionally quote applied magnetic 
fields using the SI magnetic field unit T= 10'' G. 



B. SrCuaOs 

x{T) data for three polycrystalline samples of SrCu203 
were fitted by the above QMC x*(^) simulations. Data 
for sample 1 in the temperature range from 4 tip-650K, 
shown in Fig. have been previously reported£3 Sam- 
ples 2 and 3 are new; data were obtained for these sam- 
ples from 4 to 400 K, as shown in Fig. 30 below. x{T) for 
each of the samples increases monotonically with T from 
~ 70K up to our high-T measurement limit. At lower 
T, an upturn in x{T) is seen for each sample which we 
attribute to paramagnetic impurities and/or defects. 



TABLE V. (;-factors parallel (gn) and perpendicular {g±) 
to the principal local crystal field and/or crystal structure axis 
and the powder-averaged value [(g) = -\/ (gf + ffl 4- pD/S ] 
for bulk Cu^^ and V^'* S = 1/2 species in several representa- 
tive copper and vanadium oxide compounds, and for defects 
in compounds for which a bulk ESR signal has not been ob- 
served. Samples are polycrystalline unless otherwise noted. 
The literature references are given in the last column. 



Compound 


511 




{g) Rcf. 


CuO 






2.125(5) 


213 


crystal (fe-axis) 


2.185(5) 






213 


BaCu02+i 


2.23 


2.06 


2.12 


214 




2.223 


2.041, 2.103 


2.124 


215 




2.21 


2.057, 2.12 


2.13 


215 


YzBaCuOs 


2.24 


2.06 


2.12 


216 




2.22 


2.08 


2.13 


217 


La2BaCu05 


2.229 


2.037 


2.103 


218 


CuGeOs (crystal) 


2.338 


2.064 


2.159 


219 


Sri4Cu24 04i (crystal) 


2.26 


2.05, 2.04 


2.12 


108 


powder 


2.30 


2.05 


2.14 


114 


Cao.85Cu02 






2.0796 


224 


Cu^"^ defects in: 




SrCu02 


2.25 


2.05 


2.12 


22c 


YBa2CU306.7-6.9 


2.28 


2.03 


2.12 


221 


(crystal) 










YBa2Cu30, 






2.08 


222 
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FIG. 26. Magnetic spsceptibility x versus temperature T 
for SrCu2 03 sample 1.e3 The curves are theoretical fits to 
the data by Eqs. (pH) assuming that x'''"°(T') is given by the 
theory for isolated ladders with g = 2.1. The solid curve 
shows the best fit to the data, where J'/J was allowed to 
vary, for which J'/J = 0.488(3) and J/fce = 1905(5) K were 
obtained. For comparison, the long- and short-dashed curves 
are the best fits with fixed J'/J = 0.3 and 0.7, respectively. 

We measured the low-T magnetization M versus ap- 
plied magnetic field H in detail for sample 3, as shown 
in Fig. 1^. For 25K < T < 300 K, M is found to be 
proportional to H for this sample to within our preci- 
sion. Significant negative curvature arises at 10 K and 
below, attributed to saturation of paramagnetic impuri- 
ties and/or defects. The slope of the lowest- _ff (1-3 kG) 
molar M{H) data from 2 to 10 K yielded the true x(r) 
(i.e., the low-field limit), which was fitted very well by a 
constant term plus a Curie- Weiss term 



X = Xo + 



imp 



T -B 



with parameters 



Xo = 0.52(1) X 10- 



molCu ' 



(32a) 



(32b) 



a. 



//gS'inip(*S'inip + 1) 



0.000654(9) 



3fcB 

cm^ K 
mol Cu ' 



B ^ -0.61(3) K 



(32c) 



(32d) 



where /imp is the fraction of impurities with respect to 
Cu with spin 5imp and g-factor (?imp- 
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SrCu^03 




FIG. 27. Magnetization M versus applied magnetic field H 
for SrCu2 03 sample 3 at low temperatures T . The set of solid 
curves is a global theoretical fit to all the data by Eqs. (p3[), 
containing a modified Brillouin function. 
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We have performed fits to the M{H) data in Fig. B7 



M = xo-ff + /iinpA^gimp<S'impMBSSi,„p(a;) , (33a) 



ffimp "S'impA'B 

kB{T-0) 



(33b) 



where yq a-nd 9 were fixed to the obsers^ed values in 
Eqs. (p^), i?Si„,p(a;) is the Brillouin functioned describing 
the magnetization vs field of the impurity spins, and tiiej 
parameter x has been modified from the usual formc22l 
(without 9) so that the expansion of Eqs. (^ ) for H 
gives the correct observed behavior in Eq. (32a). The 



impurity Weiss temperature 9 can arise from interactions 
between the impurity spins and/or from single-ion effects 
associated with splitting of the impurity spin energy lev- 
els. The /imp is fixed unique ly b y the observed Qmp 
and by Simp and gimp in Eq. {'i2c). Fitting the M{H) 
data at 10 K showed that gimp 



2 and Si-, 



-•imp ~ 3/2, 

lower spin values cannot give the strong negative cur- 
vature in M{H) extending up to ^ 10 K. Then fixing 
'S'imp = 3/2, the only remaining adjustable parameter is 
gimp, and from a two-dimensional global fit to all the 
M{H) data in Fig. |^ we obtained gimp = 2.093(7); the 
fit is shown by the set of solid curves in Fig. This g- 
value is within the range expected for Cu"'"^ as discussed 
above, but is slightly smaller than thaL-^2.14) found by 
ESR for magnetic defects in SrCu203.Ei3 Hence the im- 
purity spin may consist of ferromagnetically coupled Cua 
clusters. The fit is very good from 5 to 10 K, but deterio- 
rates progressively for T = 4, 3 and 2 K. Thus including 
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9 in Eq. (|33b|), which is a high-T mean- field like approx- 
imation, is not nor was expected to be accurate at the 
lowest temperatures. We have carried out fits of the ex- 
act expression (without 0) for M{H, T) of a spin S — 3/2 
impurity with the Sz levels split by a single-ion inter- 
action DS^, for which the high-T approximation gives 
9 = -4D/5. The Sz = ±1/2 levels were indeed found to 
be the ground levels (I? > 0), wit h D close to that pre- 
dicted from the observed 9 in Eq. ( |32cj ). This treatment 
improved the fit to the low-T data at large H at the ex- 
pense of a poorer fit at low H , but with an improvement 
in the overall fit. We will not present or further discuss 
such detailed fits here. 



1. Isolated Ladder Fits 

Figure |2^ shows the x(T) for sample 1 along with the 
best fit (sohd curve) by Eqs. (|3l|), where J™ax — j g ^ 
2.1 and x*(0 is our fit to the QMC data for isolated 
ladders ( J'^'^'s = 0) with spatially anisotropic exchange, 
for which the parameters are 



Isolated 2-Leg Ladder Parameters 
For Three SrCu^Og Samples 




FIG. 28. Exchange constants J'/ J (filled symbols, left 
scale) and J/Zcb (open symbols, right scale) versus fixed 
g-value assumed in the fits by Eqs. (^l|) to the data for 
SrCu2 03 samples 1 (circles), 2 (squares) and 3 (diamonds). 



Xo = -0.27(4) X 10' 



cm 
mol Cu 



0.00105(2) 



cm^ K 
mol Cu ' 



9 = -2.4(1) K, 



/ V 

— = 1905(5) K , — = 0.488(3) 
Kb J 



(34) 



We will not continue to quote xo, Cimp and 9 values 
from the fits since these were essentially the same for 
all the fits to be described. Also shown as the dashed 



curves in Fig. |26| are the best fits obtained by setting 
J'/ J at the fixed values of 0.3 and 0.7; for these very poor 
fits, J/fcB = 1870(100) K and 1827(21) K were obtained, 
respectively. 

Next, we fixed g at 2.0 to 2.2 in 0.05 increments and for 
each g-value determined the best-fit parameters for each 
of the three samples for the isolated ladder model, which 
are plotted versus g in Fig. Over the physically most 
reasonable g-value range 2.10(5) discussed above, from 
the fit parameters for all three samples taken together we 
estimate that J'/ 0.48(3) and J/Zcb ^970(150) K. 



Eccleston et 



and Azuma et oLLlI have found 
from inelastic neutron scattering measurements on the 
CU2O3 two-leg ladders in Sri4Cu2404i and SrCu203 that 
A/fcB = 377 (1) K and « 380K, respectively. As dis- 
cussed in Sec. VIIl| below, Eccleston et al. also infer that 
J'/ J = 0.55 and J/k-Q = 1510 K. Using the parameters 
of Eccleston et al., we fitted the data for our SrCu203 
samples after setting g = 2.1. The fit for sample 1 is 
shown as "Fit 1" (dashed curve) in Fig. similarly bad 
fits were obtained for samples 2 and 3. We then allowed 



Isolated Ladder Fits: g = 2.1, 
A/kg = 377 K (Eccleston et al. 
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FIG. 29. Magnetic sapceptibility x versus temperature T 
for SrCu2 03 sample 1,e3 showing theoretical fits to the data 
by Eqs. ( pl[ ) assuming g = 2.1 and a spin gap A/fce = 377 K 
which was measured using neutron scattering for the CU2O3 
two-leg ladders in Sri4Cu2404i crystals.E£ij Fit 1 is for 
fixed J and J' / J values deduced by Eccleston et al 
whereas for Fit 2 these parameters were allowed to vary, sub- 
ject to the constraint that the spin gap is given by the neutron 
scattering result A/fce = 377 K. 



29 



Isolated Ladder Fits 
12 



g = 2.1, J'^'^s/J = -0.1 



Isolated Ladder Parameters for J'^'^^/J = 0, -0.1 
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FIG. 30. Magnetic susceptibility x versus temperature T 
for SrCu203 samples IjEi 2 and 3. The curves are theo- 
retical fits to the data by Eqs. (^f|) assuming g = 2.1 and 
jdiag^j _ _Q i^ g^jjfj extrapolations up to TOOK are shown. 
The fit parameters are given in the text. 
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For Three SrCu^Og Samples 
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FIG. 31. Exchange constants J'/J (open symbols, left 
scale) and J/fce (filled symbols, right scale) versus fixed 
(/-value assumed in the fits by Eqs. (^l|) to the data for 
SrCu2 03 samples 1 (circles), 2 (squares) and 3 (diamonds); 
the fit parameters obtained for both fixed J'^'^^/J — —0.1 
and are shown. 



2. Coupled Ladder Fits 



J'/J and J/Zcb to vary, subject in Eq. (|^) to the con- 
straint that the spin gap is given by the accurate neutron 
scattering result A/fce — 377 K; the result is "Fit 2" in 
Fig. ^ (the fit parameters are given in the figure) , which 
is obviously a much better fit. 

Continuing with the isolated ladder model, we now in- 
clude the FM nonfrustrating diagonal intraladder cou- 
pling J'^'^^s in fits. As discussed in Sec. ^ above, 
our LDA-hU calculations predict that J'^'^'s/J « -0.1. 
Shown in Fig. ^is the best fit for each of the three sam- 
ples assuming fixed g = 2.1 and J'^'^'s/J = -0.1. The 
fits yielded J'/J = 0.481(4) and J/fce = 1854(6) K for 
sample 1 with a relative rms deviation (Tims = 1-98%, 
J'/J = 0.493(2), J//CB = 1883(3) K and (t„„s = 0.54% 
for sample 2, and J'/J = 0.471(2), J/k^ = 1930(4) K 
and (Trms = 0.71% for sample 3. These J'/J and J/fce 
parameters are well within the respective ranges deter- 
mined above for J'^'^s = q. To determine the sensitivity 
of these parameters to the assumed value in the pres- 
ence of a J'^^^s. = —0.1, we again fixed g at 2.0 to 2.2 
in 0.05 increments and for each (7-value determined the 
best-fit parameters for each of the three samples, which 
are plotted versus g in Fig. For the fixed parameter 
ranges g = 2.10(5) and J'^'^s = -0.05(5), from the fit pa- 
rameters for all three samples taken together we estimate 
that J'/J = 0.48(4) and J/fce = 1950(170) K, nearly the 
same as for J'^'^'s = above. Thus the fit parameters are 
not very sensitive to the precise value of J'^^^^/J, at least 
if its magnitude is much less than unity. 



It has been suggested that the frustrating trellis layer 
interladder coupling should be ferromagnetic because 
it involves Hund's rule coupling through 90° Cu-O-Cu 
bonds. eI Its magnitude was estimated to be about an or- 
der of magnitude smaller than the intraladder coupling, 
i.e. J"/ J ^ -0.1. Our LDA-I-U calculations in Sec. 
indicate an even smaller magnitude. We carried out fits 
to the x(r) data for all three SrCu203 samples by this 
model with fixed g = 2.1 but allowing all three param- 
eters J/kB, J'/J and J" / J to vary. The fitted values 
of J" /J were in the range -0.7 < J" / J < 0.9 (two of 
the three fitted values are outside the range of validity 
of the fit), showing that J" /J is too strongly correlated 
with the other two parameters to allow all three exchange 
constants to be simultaneously varied. Shown in Fig. p3 
are the best fits to the data for samples 1, 2 and 3 as- 
suming the fixed values g = 2.1 and J" / J = —0.1. The 
fit parameters are J/fcs = 1944(5) K, J'/J = 0.476(2) 
and a,ms = 1-46% for sample 1, J/kB = 2000(3) K, 
J'/J ^ 0.474(2) and a^ms = 0.65% for sample 2 and 
J/kB = 2051(3) K, J'/J = 0.455(2) and cr,ms = 0.79% 
for sample 3. 

In Fig. ^ we compare the parameters obtained for 
fixed J" / J = —0.2 and and for fixed values from 
2 to 2.2. From this figure we infer from the fit pa- 
rameters for all three samples taken together that for 
the ranges g = 2.10(5) and J" / J = -0.1(1), the in- 
traladder exchange constants are J'/J ~ 0.465(40) and 
J/kB = 2000(180) K. 
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MFT-QMC Trellis Layer Fits: g = 2.1, J7J = -0.1 
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FIG. 32. Magnetic sasceptibility x versus temperature T 
for SrCu2 03 samples 2 and 3. The curves are theoretical 
fits for the trellis layer to the data by Eqs. (E^) assuming 
g = 2.1 and trellis layer interladder coupling J /J = — 0.1. 
The fit parameters are given in the text. 
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FIG. 33. Intraladder exchange constants J'/ J (open sym- 
bols, left scale) and J/fce (filled symbols, right scale) versus 
fixed p- value assumed in the fits by Eqs. (Bll) to the data for 
SrCu203 samples 1 (circles), 2 (squares) and 3 (diamonds), 
for fixed trellis layer interladder couplings J" / J = —0.2 and 
as indicated. 

The quality of stacked-ladder fits to the data is very 
sensitive to the value of J' / J and of the interladder cou- 
pling J'" I J perpendicular to the plane of the ladders, 
because the shape of the theoretically predicted spin sus- 
ceptibility x*(0 depends strongly on these parameters 
due to the proximity to a QCP. Thus when we fit the 
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FIG. 34. Magnetic susceptibility x versus temperature T 
for SrCu2 03 sample 1.e3 The two solid curves are fits to the 
data by Eqs. ( ^l| ) assuming g = 2.1 and J' j 3 = 0.5 and 1, 
respectively. The other exchange constants determined from 
the respective fits are listed. 



Stacked Ladder Parameters for J'/J = 0.5 
For Three SrCu^O^ Samples 
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FIG. 35. (a) Exchange constants ,]"'/,] (filled symbols, 
left scale) and J/fce (open symbols, right scale) versus fixed 
g-value assumed in the stacked-ladder fits by Eqs. (^l|) to the 
data for SrCu2 03 samples 1 (circles), 2 (squares) and 3 (dia- 
monds), (b) The relative rms fit error for the fits shown in (a). 
The lines connecting the parameter data points in both (a) 
and (b) are guides to the eye. 
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experimental data by the predictions only small values 
of J'" / J can fit the data. We find that of the two possi- 
bilities J' / J — 0.5 and 1 for which we carried out QMC 
X*{t) simulations, fits with J'/ J = 1 are very poor, in 
contrast to the excellent fits obtained with J'/J = 0.5. 
Shown in Fig. 34 are the respective exchange parame- 
ters and best fits to the x{T) data for SrCu203 sample 1 
assuming g = 2.1. For either assumed J'/J value, the 
fitted J'" I J values are consistent with SrCu203 being in 
the gapped part of the phase diagram. Concentrating 
now on fits with J'/J = 0.5, the sensitivities of the fit- 
ted J'" I J and J//cb values to the assumed g- value are 
shown in Fig. ^(a) for fixed (/-values from 2.0 to 2.2. 
In contrast to other fits discussed above, the rms devi- 
ation of a fit from the data depends rather strongly on 
the assumed g-value, as illustrated in Fig. ^(b), where 
the optimum ^-values for the fits to the data for two of 
the three samples are seen to be consistent with the range 
2.10(5) we have assumed when quoting the exchange con- 
stants derived from the other fits above. From Fig. p5|(a), 
the exchange constants for the (/-value range 2.10(5) are 
J'"/ J = 0.01(1) and Jjk^ = 1920(70) K. These parame- 
ters and the fit quality are essentially identical with those 
determined for the optimum isolated ladder fit in Fig. |2^, 
so we will not plot the fit for the present case. 



C. SrzCusOs 

The x(2^) measured in H = 1 T by Azuma et a/.@ 
for the three- leg ladder trellis layer compound Sr2Cu3 05 
from 5 to 650 K is shown in Fig. |36|(a). The expanded plot 
of the low-T data in Fig. ^(b) exhibits a cusp at ~ 50 K, 
evidently associated with the short-range spin-glass-type 
ordering observed £oe this compound at « 52 K from 
/iSR measurementsEJ We fitted the data in Fig. |3^(a) 
by Eqs. (|^), where x*(0 is our global fit to the QMC 
simulation data for the three-leg 5* = 1/2 ladder with 
spatially anisotropic exchange. The best fit for g — 2.1 is 
shown as the heavy solid curve in Fig. 36(a); the spin sus- 
ceptibility contribution is shown as the light solid curve 
and the contributions xo + Cimp/(T' — d) as the dashed 
curve. The parameters of the fit are 



Xo = -0.8(2) X 10- 



cm 



mol Cu 



(35a) 



Ci„ip = 0.00053(26) 



cm'^ K 
mol Cu 



61 = -41(15) K, (35b) 



T- 0.60(4), f 
J Kb 



1814(22) K . 



(35c) 



The relative rms fit deviation (1.3%) is found to be 
nearly independent of the assumed value of g. For an 
allowed (/-value range 2.1(1), the fitted parameter ranges 
become 
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FIG. 36. (a) Magnetic susceptibility % versus tempen 
ture T for the three-leg ladder compound Sr2Cu305 (•)£ 
The heavy solid curve is a fit to the data by theory for iso- 
lated three- leg ladders using Eqs. (jsi]), yielding the intralad- 
der exchange constants J'/J = 0.602 and J/fce = 1814 K. 
The other fit parameters are given in the text. The light 
solid curve is the spin susceptibility x''^'" contribution, and 
the dashed curve the contributions xo + C'imp/(r — 6). 
(b) Expanded plot of the data in (a) at low temperatures. 
A small cusp occurs at ~ 50K as indicated, probably associ- 
ated with the sbprt-range AF ordering seen at ~ 52 K in /^SR 
measurements.Ej 



Xo = -0.8(5) X 10- 



mol Cu 



(36a) 



Ci„,p = 0.0005(4) 



cm^ K 



mol Cu 



= -41(20) K, (36b) 



J^ 
J 



0.60(7) , 



J 

fee 



1810(130) K 



(36c) 



On the other hand, the lowest-T data in Fig. |^(a) do 
not show any direct evidence for the existence of an impu- 
rity Curie- Weiss contribution. In addition, the expanded, 
plot in Fig. |3^(b) shows evidence that the reportecO 
spin-glass transition at ~ 50K affects x(^) and also sug- 
gests that there are pretransitional effects. We therefore 
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refitted the data only above 100 K in Fig. ^(a) assuming 
Cimp = 0. Tlie resulting fitting parameters for the range 
g = 2.1(1) were 



Magnetic Susceptibility of LaCuO 



Xo = -0.4(2) X 10- 



niol Cu 



J 



= 0.66(5) , 



J_ 



1810(150) K 



(37a) 



(37b) 



For the reasons mentioned, these parameters are consid- 
ered to be more reliable than those in Eqs. ( |35| ) and (36). 
The ratio J' / J thus appears to be somewhat larger and 
the value of J a little smaller than the respective values 
in SrCu2 03. 

With regard to the exchange constants obtained in this 
section, it should be kept in mind that we have implicitly 
assumed that the exchange constant along the central 
leg of the three-leg ladder is the same as that along the 
outer two legs. T his is not necessarily the case (see the 
discussion in Sec. VIIl). 



D. LaCu02 5 

x{T) data for two polycrystalline samples of LaCu02.5 
were fitted by our QMC x*(0 simulations. Data for sam- 
ple 1 in the temperato'e range from 4 to 550 K have been 
previously reported. l3 Sample 2 is new; data were ob- 
tained for this sample from 4 to 350 K. The data for both 
samples are shown as filled and open circles in Fig. re- 
spectively. Long-range AF ordering has been found from 
NMR and fiSR measurements at Tn 110-125 kJeI'EI 
as noted in the Introduction. Shown in Fig. ^ are ex- 
panded plots of the measured xC^) data for the two sam- 
ples from 60 to 160 K. Although one could perhaps infer 
the occurrence of an anomaly in each set of data in the 
range between 115 and 130 K, we conclude that there is 
no clearly defined magnetic ordering anomaly in the x{T) 
data for either of our two samples in this T range. 

We fitted the xi^) data for each of the two samples by 
Eqs. ( pT| ) , where x* {t) is our fit to our QMC simulations 
for LaCu02.5-type 3D coupled ladders with no spin gap. 
Because the x*(^) fit function is three-dimensional, we 
could vary the fitting parameters J, J^^/J and J'/ J 
simultaneously to obtain the best fit. We found that with 
any reasonable g « 2, the possibility J'/J « 1 could be 
ruled out by the bad quality of the fits to the data. Good 
fits were obtained for J' / J « 0.5. Setting J'/ J = 0.5 and 
g = 2.1 yielded a fit to the data for each sample with 
parameters 
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FIG. 37. Magnetic susceptibility 
for LaCu02.5 samples 1 (•, Ref. [77 
are theoretical fits for LaCu02.5-t2 
the respective data using Eqs. (pi 
and J'/J — 0.5, yielding the intraladder J and interladder 
J'^^ exchange constants listed for each sample. The other fit 
parameters are given in the text. 
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FIG. 38. Expanded plots of the magnetic susceptibility x 
versus temperature T for LaCu02.5 samples 1 (Ref. and 2 
from Fig. 1^ 
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Sample 1 : 

Xo = -3.0(1) X 10 

Cimn = 0.00187(3) 
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Sample 2 : 

Xo = -3.5(1) X 10 

Ci,„o = 0.00153(1 



-5 



cm 
mol 



cm'^ K 



_J_ 

kB 



mol 



t3D 

1586(14) K, — 



= -5.52(9) K , 
-0.058(1) , (38b) 



as shown by the solid curves in Fig. Allowing the 
parameter J' / J to vary during the fits yielded equivalent 
quality fits with exchange constants J/Zcb = 2690(560) K, 
J'/J = 0.562(6) and J^^/J = -0.038(10) for sample 1 
and J/kB = 1505(39) K, J'/J = 0.484(8) and J^^/J = 
—0.056(2) for sample 2, which are similar to those in 
Eqs. ( |3^ ) assuming J'/J = 0.5. The values of J^^/J 
are close to the value Jq^p/J = -0.036(1) at the QCP 
for FM values and J'/J = 0.5, and are on the 

ordered side of the QCP as expected from the observed 
AF ground state. 

We also carried out fits to the x(T) data in which we 
allowed g to vary along with the exchange constants, but 
due to the strong correlation especially between g and 
J, the estimated standard deviations on the parameters 
were very large and the fitted g and exchange constant 
parameters are therefore uninformative, but arc consis- 
tent within the errors with those given above. Finally, we 
also carried out fits to the data assuming that J'/ J = 0.5 
or 1 and that a spin gap exists in LaCu02.5, but the fits 
for each J'/J yielded J^^ / J values outside the ranges of 
validity of the respective QMC data fit functions, indi- 
cating that the assumption of the existence of a spin gap 
is incorrect, consistent with the fitting results obtained 
above assuming a gapless excitation spectrum. 



E. CaV205 

The xC^) data measured for CaV205 and CaVsOr up 
to 700 K are shown in Fig. we include data for the 
latter compound, which exhibits long-range AF ordering 
below Tn ^ 23 K according to Ref. 231 (we find Tn = 
25 K), because the former compound contains the latter 
as an impurity phase which must be corrected for. M{H) 
isotherms at 5, 100 and 200 K are shown for CaV205 in 
Fig. 
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To within our precision, M cx iJ at 100 and 
200 K, but pronounced negative curvature is apparent at 
5 K for < 1 T. We were able to fit the data at 5 K very 
well by Eqs. ( p3| ) assuming gimp = 1.96 and S'imp = 1/2, 
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FIG. 39. Magnetic susceptibility x versus temperature T 
for CaV205 (sample 1) and CaVsOr. 




FIG. 40. Magnetization M versus applied magnetic field 
H for CaV205 sample 2 at 5, 100 and 200 K. Proportional 
fits to the data at 100 and 200 K are shown. The sofid curve 
through the data at 5 K is a modified-Brillouin function fit 
described in the text. 

as shown by the solid curve in the figure. The fitting 
parameters were 

3 

Xo = 10.47 X 10"^ , (39a) 



fn 



0.0181% 



molV ' 
— -1-4. 17 K (ferromagnetic) 



(39b) 



with a variance of 0.000420 (G cm^/molV)^ 
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FIG. 41. Expanded plot of the magnetic susceptibility x 
versus temperature T for CaV205 sample 1 from Fig. 

In view of the observation of Luke et al. of an anomaly 
in xC^) 1-^"^ of spin-freezing by /iSR at ~ 50 K in 
CaV205ji£] we looked carefully at our x{T) data for this 
compound near this temperature, and indeed found a 
small but clearly defined cusp at 44 K in the measured 
data, as shown in Fig. This anomaly presumably can- 
not arise from the CaVaOy impurity phase in our sample, 
since the AF ordering transition in pure CaVaOy occurs 
at « 25 K. In the following, we limit our theoretical fits to 
the data for CaV205 in the temperature range T > 50K. 

We first fitted the data for CaVaOg from 50 to TOOK 
by 



xiT) = Xo 



imp 



T 



(!-/)(' 



J^CaV30.(y) 

cm^ K 



0.3751 



xit) 



molV/ J'/ZcB 



t = 



1 + ' 



X*{t) , (40a) 



(40b) 



J' 



X = 



J' 



(40c) 



where / is the molar fraction of the sample with re- 
spect to V consisting of the CaVaOy impurity phase, 
-^CaVaOT^ji^ is a fit to the measured susceptibility of this 
impurity phase per mole of V, x*(^) is the reduced spin 
susceptibility of a coupled dimer syst em ac cording to the 
molecular field theory (MFT) in Sec. IIID| , wi th m olecu- 
lar field coupling constant A as defined in Eq. ([40c]) where 



the sum is over all exchange coupling constants Jij from 
a given spin Si to all other spins Sj outside its own dimer. 
The x*''*'™°''(i) of the isolated dimer was given previously 
in Eq. (^ and the exchange constant within a dimer is 
here denoted by J'. 

Our fit of Eqs. ( ^ ) to the data, assuming a fixed g — 
1.96, is shown as the heavy solid curve in Fig. E2I where 
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FIG. 42. Magnetic susceptibility x versus temperature T 
for CaV205 sample 1 (•) and a fit by Eqs. (|4o|), assuming 
g — 1.96, for AF quantum S — 1/2 dimers which are cou- 
pled to each other according to the molecular field approxi- 
mation. The heavy solid curve (not visible over most of the 
T range due to overlap with the data points) is the fit, the 
light solid curve the contribution from the coupled dimers, the 
short-dashed curve the contribution from the CaVaOy impu- 
rity phase and the long-dashed curve the contributions from 
the T-independent susceptibility and impurity Curie term. 
The fit and contribution curves are extrapolated down to 20 K 
and up to 750 K. The fit parameters are given in the text. 



the contributions from the various terms in Eq. (40a) 
are also plotted as indicated in the figure caption. The 
parameters of the fit are 



Xo = 0.78 X 10"^ 
/ 



cm 



5.0% , f- 



molV ' 
= 667K , 



cm K 

Cunp = 0.0011 — 

mol V 



A = 0.31 



(41) 



The fraction / of the sample consisting of CaVaOr im- 
purity phase is close to the value of ~ 4% that we es- 
timate from our x-ray diffraction measurements on the 
same CaV205 sample. The magnitude of J' is surpris- 
ingly large for a d} vanadate. If we assume the same 
types of nearest-neighbor V-V exchange interactions as 
discussed above for the Cu-Cu exchange constan ts in 
SrCu203, then in that notation we have from Eqs. (40c) 
and (ftH) 



2( J + /" -I- J'"] 
A = + ) = 0.31 . 



(42) 



Since we obtained an excellent fit to the data using MFT 
which cannot be significantly improved upon, it is not 
in general possible to establish from the experimental 
x(r) data alone, without further theoretical and/or ex- 
perimental input, which is the V-V dimer bond and to 
which V spin(s) outside a dimer a V atom is most strongly 
coupled. 
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3D Trellis Lattice Fits to Magnetic Susceptibility 
0<J/J'<0.2, -0.2<J7J'<0.2, 0<J'7J'<0.2 
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FIG. 43. (a) Scatter plot of the relative rms fit deviation 
versus twice the sum of the V-V exchange coupling constants, 
excluding the rung coupling constant J', obtai ned f rom fits to 
the magnetic susceptibility of CaV205 by Eq. (31a) assuming 
g = 1.96 and the 3D trellis lattice model, (b) Scatter plot of 
the fitted rung coupling constant J' versus rms fit deviation 
for the fits in (a). 

We assume now that the strongest V-V exchange bond 
in the system, the dimer bond, is across the rungs of the 
two-leg ladders in the structure. We then fitted the ex- 
perimental data by Eqs. (^) using a fixed g = 1.96, but 
where the x*(^) is now given by our fits to our QMC 
simulations for the treUis layer and stacked ladders. We 
carried out 225 fits to the experimental data for the pa- 
rameter ranges < J/ J' < 0.2, -0.2 < J" j J' < 0.2 
and < J'" /J' < 0.2, in increments of 0.05 for each 
parameter. A scatter plot of the rms fit deviation ver- 
sus A in Eq. (^) is shown in Fig. ^(a) for the fits in 
which the fraction / was physical (positive), where the 
minimum in the deviation is seen to occur for A ~ 0.2, in 
approximate agreement with the estimate from MFT in 
Eq. (^) . A scatter plot of the J' values from these fits 
versus rms fit deviation is shown in Fig. ^(b); the eight 
best fits with relative rms deviations below 0.3% give 
J'/fce = 669(3) K, in agreement with Eq. (|4l)). For these 
eight fits, / = 5.3-5.7%, J/ J' and J'" /J' were either 0, 
0.05 or 0.1 and J"/J' was either -0.05, 0, 0.05or0.1, 
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FIG. 44. (a) Magnetic susc eptib ility x versus tempera 



152 and the present work). 



ture T for MgVaOg (•, Ref. 
(b) Expanded plot of the data in (a) at low temperatures. The 
solid curve in (a) is a 300-700 K fit to the data by Eq. 
which is the sum of a constant term and a Curie- Weiss term; 
extrapolations down to 100 K and up to 750 K are also shown. 
In (b), the heavy solid curve is a 1.8-13 K fit to the data 
by Eqs. (0) and (^), where the former expression is the 
low-temperature approximation for the spin susceptibility of 
a two-leg ladder with spin gap A; an extrapolation up to 20 K 
is also shown. The spin susceptibility ^'''""(r) from the fit in 
(b) is plotted as the light solid curve. The fits in (a) and (b) 
both assume g = 1.96. 

subject to the observed constraint that (J -|- J" + 
J'")/ J' — 0.1 for all eight fits. The corresponding ranges 
for xo and Cimp were 0.02 to 0.2 x 10~^ cm'^/mol V and 
0.0011 to 0.0013 cm"^ K/molV, respectively, very similar 
to the values in Eq. (pT|). 



F. MgVaOs 

The x(?^ of MgV205 sample 1 measured by 
Isobe et ali^ in a field = 0.1 T below 300 K and ex- 
tended here up to TOOK is shown in Fig. ^(a); an ex- 
panded plot of the data at low temperatures is given 
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in Fig. ||(b). A broad peak is seen at T™'^'' - 100 K, 
symptomatic of dynamical short-range AF ordering and 
of a dominant AF interaction in the compound. T™'^'^ is 
about a factor of four smaller than that for CaV205 in 
Fig. and we thus expect that the largest AF exchange 
constant in MgV205 is roughly a factor of four smaller 
than in CaV205, i.e. ~ 170 K. Of aU the QMC simu- 
lations we have presented in Sec. ||, the shape of x{T) 
for MgV205 most closely resembles that for the isolated 
ladder with J'/ J ~ 0.2 in Fig. ||. On the other hand, the 
exchange constants in CaV205 found above are evidently 
in the opposite limit J'/ J ^ 1. We begin our analysis 
with the x{T) data at high temperatures T > T™'''' in 
Fig. |4|(a). 

We fitted the x{T) data for MgV205 from 300 to 700 K 
in Fig. 44 [a) by the sum of a constant term and a Curie- 
Weiss term 



xiT) = Xo + 



C 



T-9 



(43) 



where the Curie constant C is given by Eq. ( pig ) as- 
suming g — 1.96. The fit is shown by the solid curve in 
Fig. ^(a), where extrapolations down to 100 K and up 
to 750 K are also shown. From the fit, we obtained the 
parameters 



Xo = 7.144 X 10" 



mol V 



(44a) 



-268. 6K 



-4 6* = 1074 K, (44b) 



where the relationship betwee n and the Jij exchange 
constants was given in Eq. ( p2b| ). The relative rms 
deviation for this fit is 0.20%^ If the data are fitted 
from 200 to 700 K, the fit parameters change slightly to 
Xo = 6.41 X lO^^cm^/molV and 9 = -258K, with a 
larger rms fit deviation of 0.42%. The fact that the ex- 
trapolated fit in Fig. ^(a) describes the data well almost 
down to T'"^'^ indicates that geometric frustratiaa may 
be an important consideration in this compound,Lil con- 
sistent with the sizable AF J'^"^s and J" couplings (see 
Table ^ calculated using, the LDA+U method by Ko- 
rotin et al. for MgV2 05.t£EI On the other hand, the value 
1074 K in Eq. (|44b|) is about a factor of two 



larger than calculated using the exchange constants of 
Korotin et al. in Table |^ We believe that this discrep- 
ancy arises because the fitted data are not at sufficiently 
high temperatures to be in the temperature range where 
the Curie- Weiss law holds accurately, which is typically 
at temperatures T > 10|6'|. We will see below that a 
reasonable fit to all the data from 2K to 700 K can be 
obtained using a model of coupled two-leg ladders. 

Various high-T (T > 200 K) fits of the x{T) data in 
Fig. |4|(a) by Eq. (p3|), including the two discussed above, 
all yielded xo ^i-^ySxlO^^cm'^/molV. On the other 
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FIG. 45. Expanded plots near 90 K of the magnetic suscep- 
tibility X in an applied magnetic field H — 0.1 T versus tem- 
perature T for two samples of MgV2 05. Data from Fig. ^ 
(sample 1) are shown in (a) and those for a different sample 2 
are shown in (b). Both samples show evidence for some type 
of phase transition at ~ 90K, most likely associated with an 
impurity phase. 

at 2K (0 < < 5T) that xo ~ 22 x 10"^ cm^/mol V. 
Since other measurements indicated a finite spin-gap A 
for which x'P'"(7' = 0) = 0, this xo = x{T -> 0) 
was attributed to a large Van Vleck susceptibility. We 
find here that vfT) at low T < 15K can be fitted very 
well by Eqs. (pi|), where the spin susceptibility x*(^) is 
given by the low-i approximation for the 2-leg ladder in 
Eq. (0), using a xo consistent with our range found from 
the high-T fits to x(2^)- Our inferred A is similar to 
the range of values found by Isobe et al. from different 
measurements as discussed in the Introduction. For ex- 
ample, shown in Fig. ^(b) is a 1.8-13 K fit assuming 
Xo = 5 X 10~^ cm'^/mol V and g — 1.96, for which the 
parameters in Eqs. (|l|) and (^Tj) are 



C = 0.00122 



cm'' K 
molV 



= -1.72K 



(45a) 



hand, Isobe et a/.t£3 concluded from analysis of M{H) 
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A = 0.00614- 



;ki/2 



mol V 



= 15.2 K . 

Kb 



(45b) 



There is evidence from expanded plots of x(^) near 
the peak, as shown in Fig. 45 , that some type of phase 
transition occurs near 90 K, most Hkely due to an impu- 
rity phase. This might cause an increase in xo at low T 
with respect to that at high T. 

On the other hand, one could argue that because the 
Xo we derived at high T is much smaller than inferred 
by Isobe et al. from analysis of low-T M{H,T = 2K) 
data, this difference might indicate the absence of a spin 
gap. One might then expect MgV205 to be close to a 
3D QCP; lack of Neel order would then be attributed to 
disorder effects. In that case at low T one would expect 
x(T) = xo+AT^ plus a Curie- Weiss impurity term. For a 
fit to be valid, one expects that the fit parameters should 
not be very sensitive to the temperature range of the fit. 
Therefore, to differentiate the quality and applicability 
of the gapped versus gapless fits, we determined the pa- 
rameters of the two types of fits to the experimental x{T) 
data from 1.8 K to a maximum temperature T™'^^. The 
fit parameters for both types of fits are plotted vs T™'^'^ 
in Figs. ^(a,b) and (c,d), respectively, where the specific 
expressions fitted to the data are shown at the top of the 
two sets of figures, respectively. We see that for the QCP 



fit, the x^/DOF diverges and the fit parameters change 
strongly as T^^^ increases above 5K, whereas x^/DOF 
for the gapped fit remains small and the parameters 
are essentially constant for 5 K < T™'"'' < 13 K. These 
results appear to rule out the gapless 3D QCP sce- 
nario and rather indicate that MgV205 has a spin-gap 
A/fcB = 15.0(6) K. 

We therefore carried out a fit by Eqs. (^) to our 
x{T) data assuming that the spin susceptibility x*(*) is 
given by that for coupled two-leg ladders with spatially 
anisotropic exchange for which there is a spin gap. We 
assumed a MFT coupling between the ladders, so the spin 
susceptibility is given by Eqs. (pO|), where Xoi^^J' /J) is 
our 2D fit to our QMC simulation data for isolated two- 
leg ladders. The fit is shown as the solid curve through 
the data set labeled "Isobe et al." in Fig. n7\ t ogether 
with the derived x*'^'"(T)- Also shown in Fig. [17| are the 
x{T) data up to 900 K obtained for a different sample of 
MgV205 by Onoda et a/.Jl£3 which by comparison with 
the data for the first sample illustrates the rather strong 
variability in x(^) that can occur between different sam- 
ples. Our fit to the data of Onoda et al. and the derived 
-^spm^y-j are also shown in the figure. The parameters of 
the two fits are 



Fit: 3C = 7.1 X 10-5 cm=^/mol V + C/(T - 9) + (A/VT)e-^'''^ 
Fit From 1 .8 K to T'^^'' 
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magnetic susceptibility x versus temperature T for MgV2 05. 



The fit parameters in (a) and (b) are for fits assuming a spin gap and those in (c) and (d) are for fits assuming a gapless magnetic 
excitation spectrum. 
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FIG. 47. (a) Magnetic susceptibility x versus temperature 
T for Mg V2 05 (•) as measured up to 300 K by Isobe et al. 
(Ref. 152 and tbe-present work) and for a different sample 



by Onoda et nl Pi as indicated. Each solid curve is a fit to 
the respective data set by Eqs. (^). A model for the spin 
susceptibility x°^'"(2^) is assumed in which two- leg ladders 
with a spin gap are coupled using the molecular field approx- 
imation, with fit parameters given in the text. The X'^"^{T) 
for each sample was derived by subtracting the respective fit- 
ted constant term and Curie- Weiss impurity term from the 
respective fit, and is shown for the data of Isobe et al. by a 
dotted curve and for the data of Onoda et al. by a dashed 
curve, (b) Expanded plots below 200 K of the data and fits 
in (a). 



Isobe et al: 



Xo = 0.000032(2) 



mol V' 



Cimp = 0.0015(2) 



cm^ K 
mol V 



-2.3(6) K, -^ = 141, 
Kb 



K 



— = 0.333(6) 



A = 3.56(6) , xVDOF = 6.1 x 10"" (cmVmolV) 



Onoda et al: 



Xo = 0.000016(2) 



mol V 



a„,p = 0.005(1) 



cm^ K 
mol V 



J 



-7(4) K, — = 158(2) K 



J 

11 



0.57(2) 



A = 1.62(6) , x /DOF = 7.7 x 10"" (cmVmolV) 



relative rms deviation = 1.65% 



(46b) 



relative rms deviation = 1.31% 



(46a) 



Both fitted values of J/fce and one of the fitted val- 
ues of J'/ J are (fortuitously) close to the respec- 
tive values l|44^^-^nd 0.64 in Table || predicted by 
Korotin et aLll^'t^ from LDA-I-U calculations, lending 
support to the present working hypothesis that MgV205 
has a spi n gap. The MFT coupling constant is given by 
Eq. (20c) to be A = 2J"/J, which from the exchange 
constants of Korotin et al. in Table || predicts A = 0.84; 
this value is significantly smaller than both of our fitted 
A values. Of course, the MFT-based fit is only expected 
to give accurate values of the exchange constants and A 
if |A| ^ 1. In addition, there is no way to include the 
frustrating AF diagonal second-neighbor intraladder ex- 
change coupling J'^"*s, such as in Table |l[ in the present 
modeling framework and this coupling has an unknown 
influence on the values of our fitted parameters. Finally, 
in view of the rather strong variation of the measured 
x(r) for different polycrystaUine samples of MgV205 and 
of possible impurity effects, a definitive evaluation of the 
exchange constants from x{T) data will probably only 
be possible using data for single crystals when such data 
become available. 



VIII. SUMMARY AND DISCUSSION 



We have carried out extensive QMC simulations of 
X*(t) for a large range of exchange parameter combi- 
nations for both isolated and coupled S* = 1/2 two-leg 
Heisenberg ladders, and fitted these and previously pub- 
lished QMC data accurately by interpolating functions. 
Quantum critical points were determined for both 2D AF 
stacked-ladder and 3D FM LaCu02.5-type interladder in- 
teractions between AF two-leg ladders. For each of these 
two interladder coupling configurations, but not for the 
frustrated trellis layer interladder couplings, and for each 
of J'/ J — 0.5 and 1, there is a temperature at which x* {t) 
is independent of the strength of the interladder coupling. 
It would be interesting to know if there is any fundamen- 
tal significance to this result. The dispersion relation 
for the lowest energy one-magnon excitations and for the 
lower boundary of the two-magnon continuum were cal- 
culated in the range 0.5 < J' / J < 1, the regime relevant 
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to cuprate two-leg ladder compounds. LDA+U calcula- 
tions of the exchange constants in SrCu203 were carried 
out. 

Gu, Yu and Shen have derived an expansion for the 
temperature- and magnetic field-dependent free energy of 
the 5 = 1/2 AF two-leg Heisenberg ladder for \ J / Jl\^ 1 
using perturbation theory to third order in J/J'H23 To 
compare our QMC x* (0 data fit in this parameter regime 
with the prediction of this analytic theory, we derived the 
zero-field x* (0 from their expression for the free energy, 
and the result is 
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The first term in Eq. ([4 71), corresponding to J = 0, gives 
the exact susceptibility of the isolated dimer in Eq. (0) , as 
it should. An overview oix*{t) for < J/ J' < 1 is shown 
in Fig. where for the larger J/ J' values the theory is 
not expected to apply. In fact, a pronounced unphysical 
hump in x* is seen to develop at t fa 0.2 for J/ J' > 0.6, 
so this value of J/ J' is an approximate upper limit of 
the J/J' range for which Eq. ( ^7| ) can be useful for all 
temperatures. Of course, since the HTSE of this x*{t) 
is accurate to order as we discussed in Sec. Ill B| , 

Eq. ( ^ ) can be used for any ratio of J/ J' at sufRciently 
high temperatures; however, in this case it is easier and 
simpler just to use the HTSE itself. 

Shown in Fig. ^ is the deviation versus temperature 
of the x*{t) prediction in Eq. for J/ J' = 0.1, 0.2, 
0.3, 0.4 and 0.5 from our accurate global two-dimensional 
fit to our QMC %*(*) data in the range < J/J' < 1. 
This figure shows, for example, that for Eq. ( |47| ) to main- 
tain an accuracy of 1% or better of the maximum value of 
X*(t), the value of J/ J' should not exceed about 0.3. The 
comparison we have done here illustrates that our QMC 
X* (t) data fit functions can be quite useful for easily and 
quantitatively comparing with, and evaluating the accu- 
racy of, other theoretical calculations of the susceptibility 
of spin ladders, in addition to their nominal use for mod- 
eling experimental data as we have done extensively in 
this paper. 

Our QMC x*it) simulations for -0.111 < J'^'^'s/J < 
and 0.4 < J'/ J < 0.65 (see Fig. |) suggest that the spin 
gap does not change significantly for a given J'/ J over 
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FIG. 48. Magnetic susceptibility x versus temperature T 
for the S = 1/2 AF two-leg Heisenberg ladder, obtained by 
us from the free energy calculated for J/ J' from third 

order perturbation theory by Gu, Yu and Shen.Eilj Curves are 
sliown for J/ J' = to 1 in 0.1 increments (top to bottom). 
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FIG. 49. Deviation of the magnetic susceptibility x* versus 
temperature T in Fig. S from our fit to our QMC x* (^) data. 



this range of J'^^^^ / J . We are not aware of any quan- 
titative theoretical calculations of how the spin gap of 
isolated ladders is affected by FM J'^'^s couplings. For 
AF (positive) J'^^'^s interactions (which are geometri- 
cally frustrating), Wang has found from density-matrix 
renormalization-group calculations (for T = 0) that 
the spin gap is very nearly independent of J'^^'^/J for 
< J'^^'^s/j < 0.4 and J'/J = 1, whereas, e.g., for 
J'/ J = 0.5 ttfi-soin gap strongly decreases with increas- 
ing J'^'^s/ jH3'E23 This difference in behavior arises from 
the presence of a phase boundary in the J'^'^s/ J versus 
J'/ J ground state phase diagram between the dimer sin- 
glet phase at small J^^^^ / J and a Haldane phase at large 
jdiag^ J. A similar behavior of the spin gajpzs J'^^'^s/ J for 
J'/J = 1 was obtained by Nakamura et aLtSj for the rail- 
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road trestle lattice, which topologically is a two-leg lad- 
der with one instead of two diagonal couplings per four- 
spin plaquette. Interestingly, Nakamura and Okamoto 
have found that x{T) calculated using QMC simulations 
for the railroad trestle model, the alternating-exchange 
chain model and the anisotropic two-leg ladder model, 
can all accurately fit the experimental x(^) data for 
KCuCls (witiuof course different exchange constants for 
each model) which illustrates that fits to x(^) data 
can establish consistency of a given Hamiltonian for a 
given spin system, but not the uniqueness of that Hamil- 
tonian, as discussed in the Introduction. The natures of 
the ground states and low-lying spin excitations of geo- 
metrically frustrated S ^ 1/2 two.- leg ifeisenberg ladders 
have been extensively studied.E3EMc3 

Turning now to the experimental part of the paper, the 
crystal structure of SrCu203, required as input for our 
LDA-I-U calculations, was reported here together with 
the structure of Sr2Cu305. Experimental x{T) data for 
SrCuzOa, SraCugOs, LaCuOz.s, CaVaOg and MgVaOs 
were fitted by our analytic fits to the QMC simulations 
to obtain estimates of the superexchange interactions be- 
tween the spins- 1/2 in these two- and three- leg ladder 
compounds. A summary of the exchange constants de- 
termined for the cuprate spin ladder materials is given in 
Table 0. 

As shown in Table W% o ur results confirm the prelimi- 
nary conclusion of Ref. p8|for SrCu203 based on analyses 
of x(r) data that J'/ J ~ 0.5 and J/fce 2000 K, assum- 
ing that the spherically-averaged (/-value is in the vicinity 
of 2.1. Due to the insensitivities of the calculated 
to the presence of either a weak FM diagonal intraladder 
coupling J'^'^^s or a FM trellis layer interladder coupling 
J", we could not determine either of these parameters 
from fits to the experimental xC^) data. Setting g = 2.1 



and the exchange constant ratio J'^^'^^/J = —0.1, which 
is the value in Table || predicted by our LDA-I-U calcu- 
lations, we obtained a good fit to the data for J and J' 
values which are within about 5% of the LDA-f U pre- 
dictions. The interladder J" and stacked-ladder J'" cou- 
plings determined from our LDA+U calculations and the 
J'" from fits to our experimental x(T) data are consis- 
tently found to be very small. Our theoretical esimate 
J" /J — 0.002 and our fitted intraladder exchange con- 
stants J and J' place SrCu203 within the spin-gap region 
of the phase diagram in Fig. 5 of Normand et alxB for the 
trellis layer, consistent with the occurrence of an experi- 
mentally observed spin gap. We note, however, that this 
phase diagram in exchange-parameter space assumes that 
the intraladder diagonal coupling J'^'^s — 0, whereas our 
LDA+U calculations indicate a moderately strong FM 
diagonal coupling in SrCu203. 

It is possible to obtain nearly as good a fit to the ex- 
perimental x(T) data for SrCu203 by the isolated ladder 
model (J, J' 7^ 0) for the isotropic ratio J'/ J « 1 (see, 
e.g., Ref. 137 ) as obtained for the strongly anisotropic 
J'/ J « 0.5 that we have inferred, but only if the Cu+^ 
g-factor is strongly reduced from > 2 to a value which 
has been detpMiined, and which we have confirmed, to 
be about 1.4.E^ We are not aware of any case of either a 
Cu"^^ oxide compound, including the Iwered cuprate su- 
perconductors and parent compounds j£d or a Cu+^ defect 
in any oxide system, in which the factor of the = 1/2 
Cu"*"^ ion, whether determined by ESR or inferred indi- 
rectly, is less than 2 (see, e.g.. Table 0, the discussion 
in Sec. |VII A| and Ref. |9l|). Therefore we consider the 
fit for J' /J w 1 and g w 1.4 to be unphysical. As we 
have shown in Fig. |2^, if a physically reasonable 5- value 
is used in the fit, a very poor fit by the isolated two-leg 
ladder model to the experimental x(r) data is obtained 



TABLE VI. Summary of the exchange constants for cuprate spin ladder compounds obtained by fitting the experimental 
magnetic susceptibility data by QMC simulations for the S = 1/2 Heisenberg model for various combinations of the different 
types of exchange interactions. A parameter with an asterisk (*) was fixed to the listed value or range during the fit or 
fits, respectively. The error, or fixed fitting range, in the last digit of a quantity is given in parentheses. In the third set of 
exchange constants for SrCu203, J and J' were constrained by the requirement that the spin gap A/ku = 377 K as found from 
inelastic neutron scattering measurements. Note that for the three-leg ladder compound Sr2Cu305 we have assumed that the 
nearest-neighbor exchange constant along the central leg is the same as along the outer two legs, which is not necessarily the 
case. 

J»°/J 



Compound 
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J/fcB (K) 


J'/J 


jdiag/j 


J" /J 


J'" /J 


SrCu203 


2.1* 

2.10(5)* 
2.1* 


1905(5) 
1970(150) 
1882 


0.488(3) 

0.48(3) 

0.471(1) 










2.1* 


1890(40) 


0.482(13) 


-0.1* 








2.10(5)* 


1950(170) 


0.48(4) 


-0.05(5)* 








2.1* 


2000(60) 


0.465(13) 




-0.1* 






2.10(5)* 


2000(180) 


0.465(40) 




-0.1(1)* 






2.1* 


1894(8) 


0.5* 






0.009(4) 




2.10(5) 


1920(70) 


0.5* 






0.01(1) 


SraCuaOs 


2.1(1)* 


1810(150) 


0.66(5) 








LaCu02.5 


2.1* 
2.1* 


1665(95) 
2400(900) 


0.5* 
0.53(5) 









-0.054(6) 
-0.043(15) 
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FIG. 50. Magnetic spin susceptibility x'^^'" versus temper- 
ature T for SrCu203 sample 1 (o) derived by subtracting 
a constapt term and a Curie-Weiss impurity term from the 
observedEj xiT) data. The solid curve is the x^^'^iT) derived 
from a fit to the data by a model of isolated two-leg ladders 
with the intraladder exchange constants J and J' listed in the 



figure. The filled 
data of Imai et al 



p^bols are the Knight shift K{T) NMR 
for ^'^O in the rungs of the two- leg lad- 
ders in La6Ca8Cu2404i (circles), Sri4Cu2404i (squares) and 
SriiCa3Cu2404i (diamonds), scaled to the x^^^'^{T) data. In 
the case of La6Ca8Cu24 04i, the K{T) data were scaled to 
match the x'''""(2^) data at the highest temperatures; a match 
does not occur over any appreciable temperature range. 

even for the small increase of J'/ J from 0.5 to only 0.7. 

The spin susceptibility x^'^^^iT) derived for SrCu203 
from our fits to the experimental x(T) data using 
Eq. (31a) is nearly the same irrespective of the model 
used for x^^'"(r), because the xo, Cimp and 9 parame- 
ters are nearly the same for the best fits of the various 
models to the data. The x''^'"(T) derived from the data 
for sample 1, obtained using the optimum isolated ladder 
fit in Fig. is shown as the open circles in Fig. lo-, 
gether with the x^'^'^l^) fit (solid curve). Imai et alx23 
have measured the paramagnetic ("Knight") shift K{T) 
for ^^O in the rungs of the CU2O3 two-leg ladders 
in Sri4Cu2404i, SriiCa3Cu2404i and La6Ca8Cu2404i, 
with the field applied perpendicular to the rung axis. The 
relationship between K(T) and x^^™{T) is written as 
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where 



if(T) = i^°''^ -I- if ''P'"(r) , 

^spi„(y) 



^spin^y^ ^ 2F- 



(48a) 



(48b) 



Vleck susceptibility, K^p^'^{T) is the contribution to 
K from hyperfine coupling to the spin susceptibility 
of a Cu spin and should be isotropic apart from the 
small anisotropy in the g factor, the prefactor "2" in 
Eq. (48b) comes from the two Cu neighbors of each 
rung O atom, F is the hyperfine coupling constant of 
an ^^O nucleus to each Cu spin and x*'''™ is in units 
of cm^/molCu. They determined the orbital shifts to 
be K°'^ = -0.02(2)%, -0.02(2)% and 0.04(4)% for 
Sri4Cu2404i, SriiCa3Cu2404i and La6Ca8Cu2404i, re- 
spectively, assuming that x^^"^{T = 0) = 0. We scaled 



their derived K^v'^'iT) data in Fig. 1(a) of Rcf. |12J to 
our experimental x*'''"^(^) data in Fig. IsO, as shown by 
the filled symbols in Fig. If the y^^^ T) is assumed 
to be the same as in SrCu2 03, Eq. (t48b|) yields the hy- 
perfine coupling constants F = 46kOe and 55kOe for 
Sri4Cu2404i and SriiCa3Cu2404i, respectively, which 
are similar tp-the preliminary estimate of F in Ref. 14 
of Imai et alEl The iv:''P'"(T) data for Sri4Cu2404i and 



SriiCa3Cu2404i, in which the doped-hole concentrations 
per Cu atom in the two-leg laddecs are estimated to be 
p ^ 0.06 and 0.12, respectively,EH3 scale very well with 
the x'^'^'^'iT) for (undoped) SrCu2 03, indicating that the 
shapes of x'^^^'^i^) of the two-leg ladders in Sri4Cu2404i 
and SriiCa3Cu2404i are each about the same as in 
SrCu203 (of course, this does not mean that the magni- 
tudes of x^^"^(2^) of all three compounds are the same). 
Remarkably and inexplicably, however. Fig. ^ shows 
that the K^'p^^{T) data for La6Ca8Cu2404i cannot be 
scaled to be in agreement with the x'''""('^) for SrCu203 
over any appreciable temperature range; this ohservation 
was also previously made by Naef and Wang.E23 Of the 



K is the nominally anisotropic and T-independent 
orbital shift arising from the anisotropic orbital Van 



three 14-24-41 compounds, La6Ca8Cu2404i is the one for 
which K{T) would have been expected to scale best with 
X"P'"(r) of SrCu203 because the CU2O3 two-leg ladders 
in both of these two compounds are presumably undoped. 

Our modeling of xC^) for the three-leg ladder com- 
pound Sr2Cu305 by QMC simulations for isolated 
three-leg ladders with spatially anisotropic exchange 
yielded intraladder exchange constants similar to those 
in SrCu203. However, we reiterate that in this fit we 
implicitly assumed that the exchange coupling along the 
inner leg is the same as along the outer two legs of the 
ladder, which is not necessarily the case (see below). 
Thurber et al. have concluded from ^"^Cu NMR measure- 
ments on Sr2Cu305 that a crossover of the instantaneous 
two-spin AF correlation length ^ from a ID behavior 
(f ~ 1/T) to an anisotropic 2D behaiiier (^ ~ q'^-^pJt^^ 
occurs upon coohng below ~ 300K.E13 To investigate 
how this effect might quantitatively affect our fitted ex- 
change constants, we increased the low-T limit of the fit 
from 100 K to BOOK and found no change in the fitted 
exchange constants, to within their respective error bars, 
from those obtained for the low-T fit limit of 100 K. 

The intraladder exchange constants in the two-leg lad- 
der compound LaCu02.5 are found to be J/fce ~ 1700 K 
and J' / J ~ 0.5, which are similar to those we obtained 
for the above strontium cuprate ladders, even though 
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the interladder exchange paths and types are quahta- 
tively different from those in the latter two compounds. 
The interladder couphng is found to be J'^^/J w —0.05, 
which is close to our theoretical value at the QCP, 
j3Dp/j ^ -0.036(1) for J'/J = 0.5, and is on the AF- 
ordered side of the QCP with no spin gap, as expected 
since Lapii02.5 exhibits long-range AF ordering at 110- 
125K.EIH The x{T) data for sample 1, but only up to 
~ 200 K, |5Jiere previously fitted by Troyer, Zhitomirsky 
and UedaED assuming g = 2, J' / J — 1 and a de- 
pendence of x'P'", yielding J/fes = 1340(150) K. Thus 
our results show that when the full data sets for sam- 
ples 1 and 2 are fitted by accurate QMC simulations, the 
fitted exchange constants are quite different from these 
earlier values. Our experimental xC^) data are in agree- 
ment with the previous theoretical conclusion that "due 
to the dominance of quantum fluctuations in this nearly 
critical systern_iio anomaly can be observed at the Neel 
temperature. "EJ 

It is perhaps significant that the average intraladder 
exchange constant {2J + J')/{3kB) ~ 1500K that we 
find in SrCu203, Sr2Cu305 and LaCu02.5 is about the 
same as in the layered cuprate superconductor parent 
compounds ■vmiti spatially isotropic exchange within the 
Cu02 planes.y 

In contrast to the range of intraladder exchange 
constant anisotropy J'/J ~ 0.5-0.7 in Table ^ for 
the cuprate ladder compounds, wc find that CaV205 
is essentially a dimcr compound (or perhaps a spin- 
ladder compound in the strong interchain-coupling limit 
J'/J^l). The AF intradimer exchange constant is 
J' /k-Q = 669(3) K, in agceement with the x(^) analysis of 
Onoda and NishiguchilllJ who obtained J' jk-Q = 660 K. 
These values are both larger than tiie spin gaps 464 K 
and 616 K inferred by Iwase et a/Jlla from ^^V NMR 
paramagnetic shift and nuclear spin-lattice relaxation 
rate versus temperature measurements, respectively, sug- 
gesting the presence of nonnegligible interdimer interac- 
tions. Because these interdimer interactions are so weak, 
we could not determine them unambiguously. We con- 
clude however that they must be included in order to 
obtain the best fit to the x(2^) data: J/ J', J" / J' and 
J'" I J' ^ —0.05 to 0.1, subject to the constraint that 
(J + J" + J'") I J' w 0.1. Our value for the rung cou- 
pling constant J' is about 10% larger thatt-obtained from 
LDA-I-U calculations by Korotin et aZji£j and our esti- 
mated value of the sum {J+J" + J"')/J' is a little smaller 
than the value {J-\- J")/ J' — 0.15 from these calcula- 
tions (see Table ||). However, we did not include a frus- 
trating AF second-neighbor diagonal intraladder interac- 
tion J'^'^s in our QMC simulations of x{T) (due to the 
"negative sign problem" ) , which Korotin et al. find to be 
Jdiag/j' = 0.033. Our exchange constants for CaVaOg 
strongly disagree with those estimated using "empirical 
laws" (involving only the-tiearest-neighbor distance 
and V-O-V bond angle)Eil by Millet et a/., Ell who ob- 
tained J/ J' = 0.80, J" I J' = 0.23 and J'/fce = 730 K, 



not expected to be very accurate. 

Our x(^) fits for CaV205 do not rule out a spin- 
freezing transition at 50K as in£ecred from jjSK and 
x(r) measurements by Luke et a/Jl£3 In fact, close in- 
spection of our x{T) data revealed a small but clear cusp 
at 44 K, qualitatively consistent with their x{T) data. 
In analogy with Sr(Cui_a;Zna;)203 discussed in the In- 
troduction, the spin-gap phase is evidently very delicate 
and the measurements of Luke et al. indicate that the 
otherwise singlet spin liquid ground state in CaV205 can 
be significantly perturbed or destroyed by a relatively 
small concentration of defects. 

The exchange constants in MgV205 were obtained 
by fitting the experimental x{T) data for two polycrys- 
talline samples using our theoretical results for x^^™{T) 
of S — 1/2 anisotropic two- leg Heisenberg ladders which 
are coupled together using the molecular field approx- 
imation. The average intraladder exchange constants 
we obtained for the two samples, J'/J — 0.46(13) and 
J/^B = 151(9) K, are rather close,_t£^ those in Tabic || 



predicted by LDA-I-U calculations.E£3 LDA band struc- 
ture and hopping integral calculations by Korotin et al. 
indicate that the strikingly different exchange interac- 
tions and spin gaps in CaV2 05 and MgV2 05 arise from 
the stronger tilting ofj-the VO5 pyramids in MgV205 
compared to CaV2 05.t£j There is a rather large vari- 
ability in xC^) between various polycrystalline samples 
of MgV2 05. A more definitive analysis of the magnetic 
interactions in this compound will hopefully become pos- 
sible when x{T) data for single crystals become available. 

The similarity between the J'/J (~ 0.5-0.7) ra- 
tios in SrCu203, Sr2Cu305 and LaCu02.5 demonstrates 
that the spatial anisotropy in these nearest-neighbor 
intraladder exchange constants, deduced on the basis of 
bilinear exchange only, is an intrinsic property of cuprate 
two- and three-leg ladders, irrespective of how they are 
(weakly) coupled to each other. Since the nearest- 
neighbor exchange constants in the 2D Cu02 square lat- 
tice are by symmetry necessarily spatially isotropic, this 
result indicates that higher order exchange paths than 
the nearest-neighbor Cu-O-Cu exchange path are present 
which are important to determining the magnetic prop- 
erties. Thus, as discussed in the Introduction, the (ef- 
fective) nearest-neighbor exchange constants in different 
compounds, obtained assuming bilinear exchange only, 
are not simply determined by, e.g., the distance between 
nearest-neighbor transition metal R ions and the R-0- 
R bond angle as has often been assumed (see also the 
discussion of the exchange constants in CaV2 05 above). 
This inference is consistent with theoretical studies of 
the four-spin cyclic exchange interaction around a Cu4 
plaquette in cuprate spin ladders and in high-Tc lay- 
ered cuprate superconductor parent compounds, to be 
discussed below. 

Since the initial work in 1996 indicating that the 
intraladder bilinear exchange psnstants in SrCu203 are 
spatially strongly anisotropic,E2 a great deal of experi- 
but who also include the caveat that these estimates are mental and theoretical research of various kinds in ad- 
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dition to that reported here has been done which fur- 
ther quantifies the spatial anisotropy of the intraladder 
exchange constants in two-leg CU2O3 ladders, which we 
now discuss. 

On the experimental side, Imai et al. have car- 
ried out ^^O NMR investigations of Sri4Cu2404i and 
La6Ca8Cu2404i single crystals and inferred from the ra- 
tio between the spin component of the Knight shift for O 
in the rungs to that in the legs of the CU2O3 two- leg lad- 
ders that the anisotropy in the intraladdec. exchange con- 
stants is J'/ J ~ 0.5 in both compoundsJlHil Then by com- 
paring their results with results which they obtained for 
the square lattice-antiferromagnet Sr2Cu02Cl2, for which 
J is known wcU,Eil they deduced J'/fce = 950(200) K for 
the exchange coupling along the ladder rungs. E£j These 
exchange constants are in good agreement with our re- 
sults for SrCu203.. 



Eccleston et alx23 and Katano et aLllH3 both fitted 
their magnetic inelastic neutron scattering data for the 
two-leg ladders in Sri4Cu2404i and Sr2.5Caii.5Cu2404i 
single crystals, respectively, using an assumed one- 
magnon dispersion relation E{ky) = [E'^'^^'^ sm^{kya) + 
1/2 

A^] , which is quite different from our accurate one- 
magnon dispersion relations in Fig. |2^. These two groups 
respectively obtained A/fce = 377(1) K and 372(35) K, 
iJ'n^'^/fcB = 2245(28) K and 1830(200) K, and E"'^''/A = 
5.95 and 4.9. Using their A and E'^'^^/A values, the 
A/ J versus J'/J in Eq. (§) (Rcf. H) and E'^^^'/A versus 
J'/ J (cf. Fig. ^) obtained from the dispersion relations 
of Barnes and RieraE£j for the two-leg ladder, they esti- 
mated that J'/J = 0.55 and 0.7(2) and J/fce = 1510 K 
and 1040(170) K, respectively. However, Eccleston et al. 
noted that "the low intensity of the signal away from the 
antiferromagnetic zone center means that {E'^^^) is not 
well defined." In addition, we have shown in Sec. |^ 
that two-magnon scattering contributions to the neutron 
scattering function are important near the maximum in 
the one-magnon dispersion relation and must therefore 
be accounted for when extracting the value of E'^^^^ from 
inelastic neutron scattering data. Our fit to our x(r) 
data for SrCu2 03 using Eccleston et a/.s' parameters and 
assuming g = 2.1 yielded an unacceptably poor fit. A 
much better fit was obtained by allowing J and J' to 
vary, but still subject [in Eq. (^] to the constraint that 
A is given by the accurate and reliable neutron scatter- 
ing result A/fce = 377 K, which gave J'/J = 0.47 and 
J/ks = 1880 K. These values are identical within the er- 
rors with our independently determined intraladder ex- 
change constants for SrCu2 03 in Table VI and close to 
our LDA-I-U calcL 



Regnault et al 



tion predictions. 



analyzed their inelastic neutron 
scattering data for the two-leg ladders in a large single 
crystal of Sri4Cu24 04i somewhat differently. They as- 
sumed a dispersion relation parallel to the ladders given 

by E{ky) = [(7rJ/2)2sin2(fcj^a) + A^] and obtained 
A/fcB = 370 K, J'/J - 0.50 and J/kn - 1860 K from fits 
to their data. These exchange constants are essentially 



identical with our values obtained from modeling x(T) 
for SrCu2 03. They also inferred that the inter ladder 
coupling is "extremely weak" , consistent with the very 
small intcrladder exchange constants in Table ^ that we 
infer from our LDA+U calculations and with the very 
small value in Table VI of the intcrladder exchange cou- 
pling perpendicular to the trellis layers determined from 
our fits to our experimental xC^) data. 

Sugai et al. have carried out polarized micro-Raman 
scattering experipieiits on single domains in polycrys- 
talline LaCu02.5£l3 They inferred from the energies of 
the broad two-magnon peaks in the spectra that J/fce — 
1456 K and that the intraladder exchange is nearly 
isotropic, with J'/J = 0.946. Our derived anisotropy is 
much stronger than this. They also estimated from sim- 
ilar measurements on LagCa8Cu2404i and Sri4Cu2404i 
that similarly small anisotropics J'/J — 0.95 and 1 oc- 
cur in the two-leg ladders in these two compounds, re- 
spectively, which are in strong disagreement with the 
anisotropics inferred from the above NMR and neutron 
scattering studies on these compounds, respectively. 

On the theoretical side, large-scale quantum Monte 
Carlo simulations of both the uniform and staggered 
susceptibilities of site-depleted two-leg ladders cou- 
pled in a mean-field approximation were performed 
by Miyazjaki et al. for ladders with spatiall y dstAtrp pic 
exchange.lI3 The magnetic phase diagramtZlEjCy of 
Sr(Cui_a;Zn^)203 was qualitatively reproduced, but the 
optimum doping level for maximum Tn was found to 
be ^ 10-12% and the AF ordering persisted to above 
20% doping, contrary to the experimental values of 4% 
and ~ 10%, respectively. Theoretical analyses of the 
site-diluted two-leg ladder with spatially anisotropic ex- 
change by Laukamp et al. indicated that this long-range 
order can arise at the temperatures and compositions 
found jexperimentally if J'/J = 0.5 but not if J'/J = 1 
or 5.8.0 

Greven and Birgeneau found from Monte Carlo sim- 
ulations that the observed T^{x) in Sr(Cui_a;Zn:r)203 
is consistent with the site-diluted ladders having 
J'/J 0.5, in agreement with our results for SrCu2 03, 
and a constant correlatiorL-lfingth ^/a — 18(2) over the 
experimental doping range .l3 They suggested that the in- 
tcrladder exchange coupling in the direction perpendicu- 
lar to the trellis layers is similar to that between the Cu02 
bilayers in YBa2Cu306.2, i-e. J'"/ J ~ 0.05, whereas we 
find a somewhat smaller value J'" /J — 0.01(1). We 
found theoretically that the QCP for stacked ladders 
with J'/J = 0.5 occurs at an intcrladder coupling given 
by Jq'cp/J = 0.048(2), which increases to « 0.16 for 
J'/J = 1. The fact that SrCu203 is on the gapped spin 
liquid side of the QCP requires that J'" in this compound 
satisfy < J"'/J < Jqcp/J- 

Azzouz, Dumoulin and Benyoussef have carried out so- 
called bond- mean-field theory calculations of the ;^^p'" (T) 
and the NMR nuclear-spin lattice relaxation rate 1/Ti (T) 
for both isolated and coupled two-leg ladders with spa- 
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tially anisotropic exchange.clSl For SrCu203, they ob- 
tained J/fcB = 850 K and J' / J = 0.67 by fating the ^^Cu 
l/Ti(T) measurements of Azuma et al£^ from w 100 
to 200 K by their theoretical prediction (the theory did 
not fit the data at higher temperatures). 

Using a transfer-matrix density-matrix renormaliza- 
tion group method, Naef and Wang recently computed 
the NMR 1/Ti(r) for in the rungs and for ^^^Cu 
in isolated two-leg ladders with both spatially isotropic 
(J'/J = 1, J/fce = 1300 K) and anisotropic (J'/ J = 
0.6, J/k-Q = 2000 K) exchange and compared the re- 
sults with corresponding expeiimental NMR data for 
SrCuaOa and La6Ca8Cu2404i.E£D The ^^{l/Ti){T) cal- 
culation reproduced the experimental ^^(1/Ti)(r) data 
quite well from 210 K to 720 K assuming J'/J — 1, 
whereas the temperature dependence for J'/J = 0.6 dis- 
agreed strongly with that of the experimental data. How- 
ever, poor agreement between the theory and experiment 
for i^(l/ri)(T) was obtained for both values of J'/J, al- 
though the calculation for J'/J = 1 was closer to the 
experimental data than that for J'/J = 0.6. The au- 
thors thus could not reach a firm conclusion about the 
value of J' / J : and suggested that inclusion of additional 
exchange interactions into the Hamiltonian and/or a re- 
evaluation of the hyperfine coupling tensor may be neces- 
sary to bring the theory into agreement with experiment. 

An electronic structure calculation in the local den- 
sity approximation was carried out by Miiller et al. for 
SrCu203.clll The hopping matrix elements estimated by 
fitting the bands were quite anisotropic, with a rung-to- 
leg matrix element ratio of w 0.7, from which J'/J < 1 
was expected. 

Estimates of the intraladder exchange constants J and 
J' and the interladder exchange constant J" within a 
trellis layer of SrCu203 were obtained by de Graaf et al. 
using ab initio quantum chemical calculations fotXlu206, 
CU2O7 and CU4O10 cluster segments of the layer .Efj They 
determined the values J/fce = 1810 K, J'/J = 0.894 and 
J" /J = —0.080. The interladder exchange was thus 
found to be ferromagnetic as expected from the nearly 
90° Cu-O-Cu bond angle between Cu atoms in adjacent 
ladders (see Table IV), with a value in approximate agree- 
ment with the initial estimate J"plJ = —0.1 to —0.2 by 
Rice, Gopalan and Sigrist in 1993.0 The value of J'/J ob- 
tained by de Graaf et al. is much more isotropic than the 
value we obtained from LDA-I-U calculations in Table ||. 
de Graaf et al. suggested that neglect of the frustrat- 
ing trellis-layer interladder coupling J" in the previous 
modeling of experimental data for SrCu203 may be the 
reason that values J'/J ^ 0.5 were obtained rather than 
J'/J « 1. We have demonstrated here, however, that a 
ferromagnetic interladder coupling J" /J = —0.1 has no 
detectable influence on the value of J'/J obtained from 
modeling x(r) data for SrCu203. 

Mizuno, Tohyama and Maekawa studied the micro- 
scopic origip of the superexchange in both ID and 2D 
cupratesElj They found that the hopping matrix ele- 
ments between Cu 3d and O 2p orbitals and between O 



2p orbitals depend strongly on the Madelung potential, 
which is a function of the dimensionality of the spin- 
lattice and details of the crystal structure. Their cal- 
culated exchange constants show qualitatively the same 
variations in magnitude and types of spatial anisotropics 
as observed for the ID chain, two-leg ladder and 2D layer 
cuprates discussed in the present paper. For quantita- 
tive comparisons, however, they pointed out that next- 
nearest-neighbor exchange and the four-spin ring ex- 
change interaction should be included in the calculations. 
In a study of the magnetiG-|Cxcitation spectra of two- 
leg ladders, Brehmer et alc^ concluded that "a moder- 
ate amount of ring exchange reduces the spin gap sub- 
stantially and makes equal bilinear exchange on legs and 
rungs consistent with experimentally observed (magnetic 
inelastic neutron scattering) spectra." Unfortunately, we 
were not able to carry out QMC simulations of x* [t) for 
Hamiltonians containing such cyclic four-spin exchange 
interactions to test this idea, due to a severe "negative 
sign" problem which occurs in the simulations. 

Subsequently, Mizuno, Tohyama and Maekawa exactly 
diagonalized the d-p model Hamiltonian for a CugOiy 
cluster in a two-leg ladder in SrCu203 using open bound- 
ary conditions, and then mapped the eigenenergies onto 
those of the Heisenberg model and thereby determined 
the bilinear J, J' and J'i'^s exchange constaiits and the 
cyclic four-spin exchange interaction J^^y^.tla The bilin- 



ear part of the spin Hamilto]|iiap is the same as given in 
Eq. (0) and the cyclic part if 
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piaqucttcs 

+ [Si ■ Si){S2 ■ 53) - {Si ■ 53) (^2 • Si) 

+ <S'i • ^2 + >S'2 • >S'3 + ^3 • ^4 + ^4 • 5'i 

+ Si-S^ + S2-Si + \ , 



(49) 



where the sum is over all Cu4 piaqucttcs on the two- 
leg ladder, labeled Cui-Cu2-Cu3-Cu4 around a plaquette. 
They calculated J/fce = 2260(60) K, J'/J = 0.77(12), 
jdiag/j ^ 0.015(10) and J'^^"/ J = 0.092(13) for 
SrCu203. Comparison of these exchange constants 
with our experimentally-derived ones indicates that the 
cyclic exchange interaction increases J and reduces the 
anisotropy between J' and J from the (effective) values 
that are obtained assuming only bilinear exchange inter- 
actions. The J and J'/J values are both larger than 
our values in Table || obtained from our LDA-I-U calcu- 
lations. Their J'^^^^ / J is antiferromagnetic rather than 
ferromagnetic as we found from our LDA+U calculations, 
with a magnitude roughly an order of magnitude smaller 
than our (small) value. 

Mizuno, Tohyama and Maekawa also calculated x*(t) 
for the 2x8 ladder by exact diagonalization of their spin 
Hamiltonian.ElS They found that x*(0 very sensitive 
to the presence and value of J'^^'^ for a given set of bilinear 
exchange constants. Using a set of values of the exchange 
constants within their range determined for SrCu203 
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larger than we infer for this compound from our theoret- 
ical fits to the experimental x^) data, both determined 
assuming the presence of only bilinear exchange inter- 
actions. Their proposed cyclic exchangeriiteraction is 
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FIG. 51. Magnetic spin susceptibility x'^'"" versus temper- 
ature T for SrCu2 03 sample 1 (•) derived by subtracting 
a constap^ term and a Curie- Weiss impurity term from the 
observedEJ x(r) data. The solid curve best fitting the data 
is the theoretical x'^^'^C^^) for isolated two- leg ladders with 
the intraladder exchange constants J and J' listed. The 
other solid curve is the first-principles calculation of Mizuno, 
Tohyama and Maekawa for a spin-1/2 2x8 ladder, obtained 
using a spin Hamiltonian derived for SrCu2 03 which also in- 
cludes a next-nearest neighbor intraladder diagonal interac- 
tion J'^'^s g^jj(j ^ foup-spin cycfic exchange interaction J'^^'^ 
with the values hsted.E13 Both calculations of x'^'""(r) assume 
5 = 2.1. 

cited above, they found that their calculation of y^^™{T) 
for g = 2|js in good agreement with the experimental 
x(T) dataE3 for our SrCu203 sample l.cl3 We find even 
better agreement with the experimentally determined 
X'^P'"(r), particularly for the higher temperatures at 
which the calculations of Mizuno, Tohyama and Maekawa 
are expected to be most accurate, using a slightly larger 
g value when computing x*'P'"(T) from their x*(i) cal- 
culation, as shown in Fig. ^ for g — 2.1, where our 
-j^spm j-ji-j isolated anisotropic ladders with param- 

eters in Eq. ( p^ ) is shown for comparison. The good 
agreement with the experimental data of both theoret- 
ical x'^'''"(2^) calculations, derived respectively from dif- 
ferent spin Hamiltonians which cannot be mapped onto 
each other, shows that a good fit of a theoretical x'^^™(2^) 
to experimentally-derived y^^'^^iT) data can demonstrate 
consistency of a given spin Hamiltonian with experimen- 
tal data, but cannot prove the uniqueness of that Hamil- 
tonian, as discussed for (VO)2P207 in the Introduction. 

From the calculations of Mizuno, Tohyama and 
Maekawa of x''^'"(2^) for the two- leg ladder compound 
SrCu203 rwhich fit the experimental data well as just 
discussed,cl3 the largest exchange constant is along the 
ladder legs (chains) with a value J/fce ~ 2300 K, which 
is about 50% larger than in the layered cuprates and also 



not significant in linear chain compoundsE^d Therefore, 
indirect support for the importance of the cyclic exchange 
interaction in the cuprate ladder compounds is that a 
similarly large value J/Zcb ~ 2200 K has been inferred 
from experimental x(2^) data for the linear chain com- 
pound Sr2Cu03 as discussed in the Introduction, which 
is also consistent with their theoretical predictions .E£ll 

It is useful to point out here that the physics of 5 = 1/2 
two-leg ladders with a spin gap is not modified at low 
temperatures k^T ^ min( J, J') by the presence of a 
four-spin cyclic exchange term. As long as there is a 
finite spin gap, the effective field theory is an 0(3) non- 
linear sigma model, characterized completely by two pa- 
rameters, the spin gap (magnon mass gap) A and the 
spin wave velocity c, even if a substantial four-spin cyclic 
exchange term turns out to be present. The low-T mag- 
netic properties are completely determined by these two 
parameters. Fits to experimental x(^) data at low T 
can in principle directly determine the value of A, but 
in practice this is often made difficult by a strong contri- 
bution from a Curie- Weiss impurity term at low T . The 
spin gap and the velocity c can in principle be obtained 
from QMC x(^) simulations and/or exact diagonaliza- 
tion calculations for Heisenberg ladders using calculated 
(e.g. from LDA-I-U) and/or experimentally determined 
exchange constants (where a 4-spin exchange interaction 
could be included). The Heisenberg ladder model studied 
here can thus in any case be viewed as an effective model 
for the low-temperature behavior of gapped spin ladders. 
At higher temperatures, deviations from the predictions 
of the sigma model are implicitly contained in the ef- 
fective exchange constants determined by fitting various 
models to experimental data. This explains why many 
different spin Hamiltonians can fit the same set of exper- 
imental x(2^) data for a given compound containing spin 
ladders with a mass gap. Of course, the values of the 
exchange constants determined from such fits depend on 
the Hamiltonian assumed. 

The physical consequences of an additional four-spin 
cyclic exchange interaction for several other experimen- 
tally observed propiedjies of cuprate spin ladders as well 
as layered cuprateaHlj have also been investigated. For 



two-leg ladders, Matsuda et al. found theoretically that 
the one-magnon dispersion relation E{k) along the lad- 
ders is not very sensi±iife to the presence of the cyclic 
exchange interaction. l£3 Their data for E{k) of the 
CU2O3 two-leg ladders in La6Ca8Cu2404i, as deter- 
mined from inelastic neutron scattering measurements on 
four aligned single crystals at T = 20JL, are consistent 
with the presence of this interaction,E£j but the reso- 
lutions and accuracies of the E{k) data were not high 
enough to be able to discriminate between the validi- 
ties of the bilinear exchange Heisenberg model and one 
with an additional cyclic exchange term. On the other 
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hand, Mizuno, Tohyama and Maekawa calculated that 
even a weak cyclic exchange interaction in the layered 
cuprates strongly influences the dispersioa^ut not the 



intensity, of magnetic excitation spectra.E£ll Sakai and 



Hasegawa have found theoretically that a plateau oc- 
curs for two-leg ladders at T = in the magnetiza- 
tion vs magnetic field at a magnetic moment 1/2 that 
of full Sjaturation if J'^^'^ / J is lai«er than a critical value 
0.05(4)J2£^LE,oger and DelrieuJSM Honda, Kuramoto and 
WatanabeEEJ and Eroles et aLE£3 have discussed the in- 
fluence of the cyclic exchange interaction on the pre- 
dicted magnetic Raman scattering spectrum for insu- 
lating layered cuprates. Lorenzana, Eroles and Sorella 
have concluded that this interaction must be included 
in the Hamiltonian in order to quantitatively describe 
the infrared optical absorption spectra due to phonon- 
assisted multimagnon excitations observed in the lay- 
ered cupratBp^ insulators La2Cu04, Sr2Cu02Cl2 and 
YBa2Cu306-t£j The influence of four-spin exchange on 
the properties of the S = 1/2 2D triangular lattice anti- 
ferromagnet have also been studied. 
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The four-spin cyclic exchange interaction may be im- 
portant to include in theory for and in the quantita- 
tive interpretation of, e.g., dilution experiments on lad- 
der compounds such as the system Sr(Cui_2,Zna;)203, in 
which the three Cu spins around an isolated nonmagnetic 
Zn impurity (acting as a spin vacancy) are members of 
only one Cu4 plaquette, whereas Cu spins in the bulk 
are members of two. Similarly, the spin interactions and 
correlations between Cu spins along the central leg of 
the three-leg ladder compound Sr2Cu305 may be signifi- 
cantly different than along the outer two legs, even apart 
from effects expected from nearest-neighbor interactions 
and the different numbers of Cu nearest-neighbors of Cu 
spins in the central and outer legs, because each spin 
in the central leg is a member of four Cu4 plaquettes, 
whereas each spin in the outer two legs is a member of 
only two. Since the nuclear-spin lattice relaxation rate 
l/Ti is a direct measure of the imaginary part of the 
local spin susceptibility at low frequency and is not yet 
fully understood in the cuprate spin ladder compounds as 
discussed above, it will be interesting to compare experi- 
mental data for 1/Ti of both O and Cu in the cuprate spin 
ladders with corresponding calculations of this quantity 
in which the four-spin interaction is included. 



In conclusion, we have demonstrated that when the 
magnetic spin susceptibilities of cuprate spin ladder com- 
pounds are analyzed in terms of the bilinear Heisenberg 
exchange model, strong anisotropy between the intralad- 
der exchange constants in the legs and rungs is consis- 
tently found, which confirms and extends the conclusion 
of Ref. ^ for SrCu203. This anisotropy strongly violates 
the conventional empirical rules for nearest-neighbor su- 
perexchange interactions in oxides. It is not yet clear 
to us whether the nearest-neighbor exchange is really 
as anisotropic as we deduce, which however is corrob- 
orated by a number of calculations (including our own 
LDA-I-U calculations) and experimental inferences enu- 
merated above, or whether the values we derive are ac- 
tually "effective" values which indirectly incorporate the 
effects of additional terms in the spin Hamiltonian as 
indicated by other calculations and experiments. Analy- 
ses of the temperature-dependent NMR Knight shift and 
1 /Ti in terms of the nearest-neighbor Heisenberg model 
have respectively yielded contradictory results regarding 
the anisotropy, and the analyses of Raman scattering ex- 
periments have indicated a much smaller anisotropy than 
we have inferred. Much work remains to be done to es- 
tablish a spin Hamiltonian and calculational procedures 
which can self-consistently describe the x(T'), NMR, opti- 
cal and inelastic neutron scattering measurements prob- 
ing the magnetism of the cuprate spin ladder materials. 
If additional terms are in fact present in the Hamiltonian 
which are important to determining the magnetic prop- 
erties, the evidence to date suggests that higher-order 
superexchange processes are likely candidates. In the 
case of the magnetic spin susceptibility studied here, the 
presence of a four-spin cyclic exchange interaction can ex- 
ert a strong influence on the strengths of the (effective) 
nearest-neighbor exchange interactions inferred by ana- 
lyzing experimental x(^) data assuming only bilinear ex- 
change interactions, as we have discussed. Such four-spin 
cyclic exchange processes may also strongly contribute 
to various properties of the layered cuprate high-Tc su- 
perconductors and undoped parent compounds as indi- 
cated by recent work on their optical properties. It would 
be very interesting to determine how this interaction in- 
fluences the calculated superconducting correlations in 
doped spin ladder compounds and also in the doped lay- 
ered cuprate high-Tc materials. 
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APPENDIX: QMC x* (t) SIMULATION FITS 

In the following weighted fits to the QMC x*(i) sim- 
ulation data, the quality of a fit to a data set is ex- 
pressed as the statistical per degree of freedom, de- 

fined as x'/DOF ^ (N^ - P)-' E^^.w^ix* " X*'''*)', 
where A^p is the number of data points in the data set, 
P is the number of fitting parameters, the weighting 
function Wi = ^/cff, and Oi is the estimated error for 
the i*'' data point. An additional measure of the qual- 
ity of a fit is the absolute rms deviation (Tims of the fit 
from the data. The fits were carried out on Macintosh 
PowerPC G3 computers (233 and 400 MHz) using the 
software Mathematica 3.0. 



1. Isolated Ladders, Isolated Ladders with 
Ferromagnetic Diagonal Coupling, Stacked Ladders 
and 3D-CoupIed Ladders with J' / J < 1 and Spin 
Gaps 

All of these QMC data were fitted by a single multidi- 
mensional function with 7'|^gj' in Eqs. (^7|), where 
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with J'"'^'^ = J and A;5 EE Ao/J from Eq. (0). The nota- 
tion, e.g., J("/J, means that the respective fit parameter 
coefficient and fit for J'" j J apply only for J' / J = 1. In 
most of the fits, the parameters Di, D2 and D3 were de- 
termined from A^i, A^2, A3 and A^^ by the three HTSE 
constraints in Eqs. (|l8|) and were thus not fitted. For 
ease of implementing our fit functions by the reader, we 
have included the values of the constrained parameters 
in the tables of fitted parameters. 



a. Isolated Ladders with Nearest- Neighbor Couplings 

We first fitted the data for the isolated ladders ( J" = 
Jill ^ Jdiag ^ j3D ^ 0) in Fig. I and additional data not 
shown (see figure caption). Here J™'*'^ = J , Obtaining a 
reliable fit was essential because the fit is used as a foun- 
dation or backbone of the fits to the data for the cases be- 
low incorporating additional exchange interactions. The 
spin gap Agj(J'/J) was set identical to Ao{J' / J) in 
Eq. (R) and thus was not fitted. For J'/ J = for which 
Ao(J7J) = 0, we used the {Nno,D,,o} fit ("Fit 1") 
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Isolated Two-Leg Ladders: 0.1 < J'/J < 1 
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FIG. 52. Semilog plot versus temperature T of the devia- 
tion of the fit from the 378 QMC x* (t) data points for S = 1 /2 
spatially anisotropic Heisenberg ladders with 0.1 < J'/J < 1. 
The error bars are the estimated accuracies of t he QMC data. 
The fit function is given in Eqs. ( |l7[ ) and (Al) with parame- 
ters in Table Ivii. 



parameter set in Vr, 



in Eqs. (|T^ which was obtained 
t < 5 in Ref. 183 for the uni- 
fitted the 30 
for J'/J = 1 

(5) 



,(5) 
(6) 

for the range 0.01 < 
form chain. In addition, we separateljt, 
QMC data points of Frischmuth et al.tZa 
using vj^Q^ in Eqs. (17); we found that using P^g^ ic- 

sulted in a {Nin{J' / J = 1), Din{J'/J = 1)} pa- 
rameter set with several very large parameters (3> 1) 
which would not allow accurate high temperature ex- 
trapolations. We thereby obtained a {Nin{J'/J ~ 
1), Din{J'/J = 1)} nine-parameter set for J'/J — 1, 
which yielded xVDOF = 0.079 and a^^B = 4.76 x 10"^ 
this fit will also be used as one of the end-function fits 
for the range J'/J > 1 in later sections. The fit pa- 
rameters for the general data set 0.1 < J'/J < 0.8 (348 
data points) were constrained by these J'/J = and 
J'/J = 1 fit parameters. The complete parameter set 
{Ninj, Dinj} for this range was then determined and is 



given in Table VII , including the constrained parameters 
Di, D2 and D3. Note that from Eqs. ([l8|), D2 and D3 are 
sixth and ninth order in J'/J, respectively. The good- 
nesses of fit for the 348 data points with 0.1 < J' / J < 0.8 
were xVDOF = 1.07 and a,^s = 0.000159, and for the 
total 378 point data set with 0.1 < J'/J < 1 were 
xVDOF = 1.02 and firms = 0.000154. The deviations 
of the fit from the data are plotted versus temperature 
in Fig. The two-dimensional fit is shown as the set 
of solid curves in Fig. |^, including an example of an 
exchange parameter interpolation x*(t, J'/J = 0.9) for 
which we did not obtain QMC simulation data. The fit 
can also be interpolated in the range < J'/J < 0.1. 
Each of the following fits in this section were done with 
the parameter set {Nno, Ninj, DnOi Dinj} found for 
the isolated ladders fixed, and with all other exchange 
parameters except the one being fitted set to zero. 
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Isolated Ladders with Diagonal Couplings 
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FIG. 53. Semilog plot versus temperature T of the de- 
viation of the fit from the 457 QMC X*(*) data points 
for S = 1/2 spatially anisotropic Heisenberg ladders with 
0.4 < J'/J < 0.65 and with ferromagnetic diagonal couplings 
jdiagyj ^ _o.05, -0.1 and -0.11. The error bars are the 
estimated accuracies of the QMC data. The fit fun ction is 
given in Eqs. (^) and (Al) with parameters in Table VII . 



b. Isolated Ladders with Ferromagnetic Diagonal 
Second-Neighbor Couplings 

For the fit to the 457 QMC data points for the iso- 
lated ladders with 0.4 < J'/J < 0.65 and with fer- 
romagnetic diagonal couplings J'^^^^/ J — —0.05, —0.1 
and —0.11, we coul d on ly use the HTSE constraints on 
Di and D2 in Eqs. (18a) and ( |l8b| ), respectively, because 
for D3 in Eq. (18c) the diagonal co uplin g violates a con- 
dition for the validity of da in Eq. (12c) that no second- 
neighbor couplings occur. The fit yie lded the parameter 
set {A^2nj, D2nj} listcd in Table VII. Because the £'22^ 
parameters are constrained, i.e. determined by other pa- 
rameters according to Eq. (181), the D22 expansion has 
nine terms which are designated in the table by 



D22 = D221 
+ -D223 
£'224 



jdiag 
J 

J' \ / J"^'"^^ 



jdiag . 2 



J 



(7)( . 



(t) 

/ J' \ /J'\^ 

D2..{j)+D22s{j) 



£'225 \ ~J J + -^226 



(^) 



(A2) 



where t he fi rst three parameters are the same as defined 
in Eq. (^). For this fit xVDOF = 0.91 and cr,.ms 
0.000237. The fit for J'/J = 0.4 and 0.6 and J'^^'^^/J = 
—0.05 and —0.1 is shown as the set of solid curves in 
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Fig. 1^ along with the corresponding data and fit for 
jdiag/j ^ rpj^g deviations for all the J'^'^s ^ q data 

are plotted vs temperature in Fig. The fit should not 
be extrapolated into the antiferromagnetic diagonal cou- 
pling regime because poles develop in the fit function at 
low temperatures for J'^^^^/J > 0.05. 



c. Stacked Ladders 

As shown in Fig. |l^, the stacked ladders with J'/ J — 
0.5 have a QCP at J'" /J = 0.048(2). The two regions 
on either side of the QCP require separate fits. We fitted 
the 84 data points for gapped ladders with J'" / J = 0.01, 



0.02, 0.03 and 0.04 and obtained the {N^^j, D^nj, A32} 



parameter set listed in Table |VII| . The value of the fitted 
parameter A32 predicts that the spin gap should decrease 
to zero at J'n np/J = 0.052, close to the value inferred 
from Fig. The deviations of the fit from the data 
are plotted vs t in Fig. [EJ. The goodnesses of fit are 
VDOF = 1.39 and cr.^s = 0.0000464. The relative 



X 



In Table VII, the constrained 



rms deviation is 0.125%. 
parameters Z33„ with n = 1, 2 and 3 are written as power 
series, = J2j -^3nj{J"' / jy ■ The fit is shown as the 
set of solid curves through the data for J'" /J — 0.01-0.04 
in Fig. |. 

We fitted the 70 stacked ladder data points for J' /J = 
1 and J'" I J = 0.05, 0.1 and 0.15, parameters which are 



TABLE VII. Parameters in Eq. (A.1) obtained by fitting quantum Monte Carlo simulations of x*(^) for isolated ladders, 
stacked ladders and LaCu02.5-type 3D-coupled ladders with J' /J < 1 for parameter regimes which have spin gaps. 
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Nln2 


Nl„3 






Dln2 


£•1713 


1 


—0.05383784 


—0.67282213 


0.03896299 


0.01103114 


0.44616216 


— 0.82582213 


0.03896299 


—0.08786886 


2 


0.09740136 


0.12334838 


—0.0253489 


0.00655748 


0.32048245 


— 0.40632550 


0.20252880 


— 0.03801372 
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stacked Ladders with J'/J = 0.5, J"7J = 0.01-0.04 



Deviation of Fit from QIVIC x*(t) Data 



3D-Coupled Ladders with J7J = 0.5 and a Spin Gap 
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FIG. 54. Semilog plot versus temperature T of the de- 
viation of the fit from the 84 QMC X*(*) data points for 
S — 1/2 stacked spatially anisotropic Heisenberg ladders with 
J'/J = 0.5 and AF interladder couplings J'" /J = 0.01, 0.02, 
0.03 and 0.04. The error bars are the estimated accuraci es o f 
the QMC data. The fit fu nctio n is given in Eqs. (|l|) and ( Al ) 
with parameters in Table VII 



Stacked Ladders with J'/J = 1 , J"'/J = 0.05-0.2 
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FIG. 55. Semilog plot versus temperature T of the devia- 
tion of the fit from the 70 QMC x* (i) data points for S' = 1 /2 
stacked spatially isotropic Heisenberg ladders with J'/J = 1 
and AF interladder couplings J'" / J — 0.05, 0.1, and 0.15 for 
which there is a spin gap (•); here arms — 0.000384. The open 
squares are the deviations of a separate fit to the 24 data 
points for J'" /J = 0.2 from those data (cTrms = 0.000256), 
for which there is no spin gap. The error bars are the esti- 
mated accuracies of the QMC data. The fit function for the 
gapped fit is given in Eqs. and (Al) with parameters in 
Table Wul 
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Deviation of Fit from QMC 3c*(t) Data 
J2°/J = -0.02 to -0.035 
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FIG. 56. Semilog plot versus temperature T of the de- 
viation of the fit from the 52 QMC X*(^) data points for 
S — 1/2 3D-coupled spatially anisotropic Heisenberg lad- 
ders with J'/J — 0.5 and ferromagnetic interladder couplings 
J^°/J = -0.02, -0.025, -0.03, and -0.035 for which there 
is a spin gap (•). The error bars are the estimated accuracies 
of the QMC data. Th e fit function for the gapp ed fit is given 
in Eqs. and (Al) with parameters in Table VII . 



deviations of the fit from the data are shown as the filled 
circles in Fig. [s^. The fitted value of A42 predicts that 
Jqqp/J ~ 0.155, at which the spin gap vanishes. 

d. LaCu02.5-type 3D-Coupled Ladders 

For the LaCu02.5-type 3D-coupled ladder data for 
J'/J = 0.5 and = -0.035, -0.03, -0.025 and 

—0.02 (52 data points), which exhibit a spin-gap, we 
fixed the parameter A52 = —145, which yields a QCP 
at J^v/J = -0.0384. We obtained a xVDOF = 1.44 
and (Trnis = 0.0000932 for the {N^nji D^nj, A52} pa- 
rameter set listed in Tabl e |VII| . The fit is shown as the 
set of solid curves in Fig. [14| for these parameter combi- 
nations, and the deviations of the fit from the data are 
shown in Fig. pq. Similarly, for the 84 da±a points for 
J'/J = 1 and J^/J = 0.05, 0.1, and 0.11,EZl which also 
exhibit a spin-gap, we obtained a x^/DOF — 0.92 and 
o'rms = 0.00012 for the {N^nj, Dgnj, A62} parameter 
set listed in Table |VI]| . The fit is shown as the set of 
solid curves for the respective exchange constant combi- 
nations in Fig. 16, and the fit deviations are plotted vs 
temperature in Fig. 



2. Trellis Layers with J'/J < 1 



all in the gapped region with J'" / J < Jqqp/J ~ 0.16, 
and obtained xVDOF = 1.19 and cr„„s = 0.000384 for 



the {N^nj, D^nj, A42} parameter set listed in Table VI] 



The two-dimensional fit is shown as the set of solid curves 
for these exchange constant combinations in Fig. O. The 



The 136 trelfis layer QMC data points with J" = 
-0.2, -0.1, 0.1 and 0.2 for J'/J = 0.5 and 1 in Fig. | 
wexe computed using QMC simulations by Miyahara et 
al.t3 They also computed data for J" = —0.5 and 0.5 



51 



3D-Coupled Ladders with J7J = 1 and a Spin Gap 
15 



10 
5 


E -5 



g 

■> 
Q 



-10 
0.03 



Deviation of Fit from QMC x*(t) Data 
J^°/J = 0.05, 0.1, 0.11 



a = 0.00012 

rms 



0.1 



0.3 
k^T/J 



FIG. 57. Semilog plot versus temperature T of the de- 
viation of the fit from the 84 QMC X*(*) data points for 
S = 1/2 3D-coupled spatially anisotropic Heisenberg ladders 
with J'/ J = 1 and AF interladder couplings J^^ / J = 0.05, 
0.1 and 0.11 for which there is a spin gap (•). The error bars 
are the estimated accuracies of the QMC data. Th e fit func- 
tion for the gapped fit is given in Eqs. ( p^ and (Al) with 
parameters in Table VII, 



which are not shown in the figure. Due to the "nega- 
tive sign problem" arising from the geometric frustration 
in the interactions between spins on adjacent ladders, 
accurate data could not be obtained to low tem- 

peratures. In addition, for the trellis layer the spin gap 
is expected to be nearly independent of the interladder 
coupling. For both of these reasons, we fitted the data 
by expressions having the form of the modified MFT ex- 
pression in Eqs. (p^). For the case of trellis layers with 
J' /J < 1, one has J™^'' — J and we write the function 
f{Jij,t) in Eq. (24b) in the general form 



J ' J 
+ 



-ic33 



„s 3 



+ 



3t^ 



C43 



(A3) 



in which the Cij parameters may be varied and used 
as fitting parameters. For the exact HTSE to 0{l/t^) 
on the right-hand-side of Eq. (|A3|), one would have 
cii = C23 = C33 = C35 = 1 with the remaining param- 
eters being zero. In their fits for —0.5 < J" /J < 0.5 and 
J'/ J = 0.5 and-I in the temperature range 1 < i < 1.5, 
Miyahara et alSB set cn = 1 and then obtained C12 — 
0.3436. Since the correct Weiss temperature in the Curie- 
Weiss law requires cn = 1 and C12 = 0, this parametriza- 
tion does not give the correct Curie- Weiss behavior but 



does yield the correct Curie law in the limit of high tem- 
peratures. 

In our fits to the QMC data in Fig. var ious combi- 



nations of nonzero Cy parameters in Eq. (A3) were tried 



The lowest-order, zero-parameter MFT fit with cn = 1 
and with the remaining Cy parameters being zero yielded 
the fit shown as the set of solid curves in Fig. For this 
"fit" to all the data points, the goodnesses of fit were 
xVDOF = 98.7 and cjr-ms = 0.0020. Excluding the two 
data points with J' /J ~ 0.5, J" /J — ±0.1 at i = 0.1 gave 
a much lower a^-^g = 0.00080. This is the fit function we 
use to fit our experimental x{T) data for SrCu203. We 
also tried various combinations of up to four Cij fitting 
parameters in fitting the data for t > 0.25, with or with- 
out one or more of the constraints above associated with 
the HTSE, but the (Tims of the fit could not be signif- 
icantly improved compared to that of the above MFT 
"fit" with zero fitting parameters. In the attempts with 
one or more nonzero 04^, the fits diverged significantly 
from the trend of the data at low temperatures because 
of the lack of enough data points at low temperatures to 
constrain these parameters. 



3. Stacked Ladders with No Spin Gap 

Our 118 QMC data points for the stacked ladders with 
J'/J = 0.5 and J"7J = 0.05, 0.1, 0.15 and 0.2 are on 
the side of the QCP at Jqcp/J = 0.048(2) with no spin 
gap and were fitted by the expression with J™^'^ = 
in Eqs. (|l^, with 



and AS, 



/ J'" \ I J'" \ 2 

N-a - iV„o + (— ) + ^"2 (— ) , (n = 1 - 5) 

/ J'" X / V" X 2 

A. = A.o + All (— ) + A.2 (— ) . (n 1 - 6) (A4) 

The Z3i, D2 and D3 series coefficients were not fitted 
but were determined from the A^i, iV2 and A'3 series ac- 
cording to the three HTSE constraints in Eqs. (p^). The 
low-T expansion is of th,e_coEi:ect form for a 2D quantum 



critical point, x* oc t\ 



184-186 



on the side of the QCP 



with no spin gap, we find that x* continues to be linear 
in t at low t, as is al so pre dicted. The parameters of the 
fit are given in Table VIII, for which x^/DOF = 1.09 and 
CTrms = 0.000431 were obtained. The fit is shown as the 
set of four solid curves for J'"/ J = 0.05, 0.1, 0.15 and 0.2 
in Fig. H Note that the fits for all four J'"/ J values cross 
at the temperature t « 0.16; i.e., at this temperature the 
X* is nearly independent of J'" / J in the gapless regime. 
The deviations of the fit from the data are plotted vs t 
in Fig. 

Our 24 QMC data points for the stacked ladders with 
J'/J ^ 1 and J'" /J = 0.2 are also on the side of the QCP 
at Jqqp/ J ~ 0.16 with no spin gap and were therefore 



also fitted by the eight-parameter expression in Eq. ( A4) 



with J'"^'^ = J, V)q} and A*, = in Eqs. (0), with 



52 



stacked Ladders with J'/J = 0.5, J"'/J = 0.05-0 
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FIG. 58. Semilog plot versus temperature T of the devi- 
ation of the fit from the 118 QMC X*(i) data points for 
5=1/2 stacked spatially anisotropic Heisenberg ladders with 
J'/J = 0.5 and AF interladder couplings J'" /J = 0.05, 0.1, 
0.15 and 0.2 for which there is no spin gap. The error bars are 
the estimated accuracies of the QMC data. The fit function 
is given in Eqs. (|l7|) and (A4) with parameters in Table VIII 



parameters shown in Table |^ The three Z?i, D2 and 
parameters were not fitted but were determined from 
the Ni , N2 and N3 fitting parameters as constrained by 
the HTSE. The fit is shown as the solid curve in Fig. |lj 
and the deviations of the fit from the data as the open 
squares in Fig. |5^. The goodnesses of fit are x^/DOF = 
0.60 and cr„ns = 0.000256. 



One expects that x* oc at the 3D quantum critical 
point, and for further increases in IJ^EVJ] one expects 
X* = x*(0) + bt'^ with a finite x* (0)1230 Thus we fitted 
all of the above QMC data using Ag^ 



in Eqs. (Il7|), with 



and using vj^^}^ 



'J' 

Nn = NnO + iV„i ( y 



' J'\2 



iV„, — 



3D 



JJ\ J 



"-(7)(V)- 



D„ = D 



Dm 



J' 



Dr, 



J3D 2 
J J 

= 1 - 3, 5) 



(A5) 



n3 



J 



3D 



J 



Dm 



J3D 2 



J 



The three HTSE constraints in Eqs. ( |l8[ ) were enforced, 
so the series for £)„ with n = 1, 2 and 3 are determined 
from A^i, A^2 and/or iVa series, and Eqs. @ with A* = 
yield the £)„ series terms in Eq. (A5) plus the respective 
additional terms 



4. LaCu02.5-Type 3D Interladder Couplings 
with J'/J = 0.5 0.6, 0.7, 0.8, 0.9 
and 1.0 and with No Spin Gap 

The LaCu02.5-Type 3D interladder couplings —0.1 < 
J3D/J < _o.035 and presumably 0.05 < J^^/J < 0.2 
wuth J'/J = 0.5, and 0.12 < J" /J < 0.2 with J'/J = 
1,EZI are on the side of the QCP with no spin gap and a 
finite x*(^ = 0). We estimate that our QMC data with 
J3D/J = 0.05 and J'/J = 0.6 and with J^^/J = 0.1, 
0.15 and 0.2 and J'/J = 0.6, 0.7, 0.8 and 0.9 are all in 
the gapless regime. These data total 564 data points. 



TABLE IX. Parameters in Eq. (A4) obtained by fitting 
quantum Monte Carlo simulations of x* (t) for stacked ladders 
with J'/J — 1 which have no spin gap. 



-0.6740395 
0.3915161 
0.007222304 
0.01429191 
0.005028154 




0.1759605 
0.2035826 
0.2061725 
0.1497920 
-0.01677125 
0.05307760 



TABLE VIII. Parameters in Eq. (A4) obtained by fitting quantum Monte Carlo simulations of x*(i) for stacked ladders with 



J'/J 


= 0.5 which have 


no spin gap. 










n 




iVnl 






Dm 


Dn2 


1 


-0.07395387 


-1.1046032 


-7.7000424 


0.5510461 


-0.6046032 


-7.7000424 


2 


0.17151127 


-1.16738784 


5.24878372 


0.4065401 


-1.8947418 


0.1339557 












D23: -3.8500212 




3 


-0.00486493 


0.02453618 


-0.53626512 


0.1700718 


-0.9299952 


-0.0235946 












Dzz: 2.4315744 


Dsi- -1.9250106 


4 


0.001860855 


-0.019272190 


0.082869123 


0.05484099 


-0.53559069 


0.82169392 


5 


-0.0003207081 


0.0074168695 


-0.0175909734 


0.001093210 


-0.051478243 


0.223840382 


6 








-0.000527212 


0.025391112 


-0.063286649 
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3D-Coupled Ladders with No Spin Gap 
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FIG. 59. Semilog plot versus temperature T of the de- 
viation of the fit from the 564 QMC X*(^) data points for 
S = 1/2 3D-coupled spatially anisotropic Heisenberg ladders 
with J'/J = 0.5 and J^°/J = -0.1, -0.08, -0.06, -0.04 
and 0.05, J'/J = 1 and J™/J = 0.12, and J'/J = 0.5, 0.6, 
0.7, 0.8 and 0.9 and J^°/J = 0.1, 0.15 and 0.2 (•), for which 
there is no spin gap. The error bars are the estimated accu- 
racies of the QMC data. The fit function for the gapped fit 
is given in Eqs. (|l7|) and (|A5|) with parameters in Table 
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This fit function satisfies the low-temperature Umit pre- 
diction x(T) = A + BT^. The 32 parameters of the fit 
in Eqs. (A_5) and the series in Eqs. (A6) for the con- 
strained parameters are given in Table for which 
xVDOF = 0.95 and a,ms = 0.000214 were obtained. 
The three-dimensional fit is shown as the sets of solid 
curves in Fig. [l8| an d for the gapless parameter regimes 
in Figs. |lj and |l6[ The deviations of the fit from the 
data are plotted vs temperature in Fig. M. 



5. Isolated Ladders, Trellis Layers and Stacked 
Ladders with J'/J > 1 



Since these spin lattices for all the exchange constant 
combinations with J'/J > 1 that we simulated have a 
spin gap, we were able to fit all the QMC data with the 
same multidimensional function, using T^^gj* in Eqs. (|l^), 
with 

Nn = iV„o + A^lnl (-^) + ^1"2 (-^) ^ + A^l„3 (^) ^ 
/ J \ 4 / J" \ / J" \ 2 

I J'" \ I J'" \ 2 

+ ^-(-)+^-2(-) 

+ ^-3(;^)(^),(n^l-6) 

/ J \ /J\2 I J \ 



+ i?l„4(;^)%i^2„t(^)+i?2„2(5)' 
/ J'" X / J'" X 2 

+ ^-3(;^)(^),(n=l-6) 



— - 1 + An(- 



J' 



J\2 



J\3 
J' 



with J™^'^ = J'. AH of the fits in this section utilized the 
three respective HTSE constraints on Di, D2 and D3 
in Eqs. ([l8|). The seven A^„o and Z?„o fitting coefficients 
were first determined to high accuracy by fitting to the 
exact expression for the 5=1/2 dimer given in Eq. (H), 
for which J/J' — O and A/ J' — 1. We were able to fit a 
498-point double-precision representation of this function 
for 0.02 < t < 4.99 to a variance of 1.2 x lO'^^. The fit 
parameters were further constrained by requiring that 
the {Nn, Dn, Afitl value s for J/ J' = 1 be identical with 
those found in Sec. A 1 a| above for J' / J = 1 from the fit 
to the QMC data for J'/J < 1. 

The 29 Ninj, Dinj and Aij coefficient fitting param- 
eters for the isolated ladders were determined by fitting 
to the 119 QMC simulation data points for J/ J' = 0.1, 
0.2, 0.3, 0.4, 0.5, 0.7 and 0.9 and the results are given 
in Table yielding a x^fDOF = 1.24 and a,.n,s = 
4.19 X 10^ for the fit to those data. The fit is shown 
as the set of solid curves in Fig. ^, including two curves 
showing exchange parameter example interpolations for 
X*{t, J/ J' = 0.6, 0.8) for which we did not obtain QMC 
simulation data. The fit deviations from all the fitted 
data are plotted vs temperature in Fig. |6^. In the fol- 
lowing fits, the {iV„o, Dno, Ninj, Dinj, Aij} set of pa- 
rameters for the isolated ladders was held fixed. 
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Isolated Ladders with 0.1 < J/J' < 0.9 
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FIG. 60. Semilog plot versus temperature T of the de- 
viation of the fit from the 119 QMC data points 
for S = 1/2 spatially anisotropic Heisenberg ladders with 
0.1 < J/ J' < 0.9. The error bars are the estimated accu- 
racie s of the QMC data. The fit function is given in Eqs. ( |l^ ) 
and (A7) with parameters in Table XI. 



The 18-parameter three-dimensional fit to the 162 
QMC treUis layer x*(0 data points with intraladder 
couplings J/ J' = 0.1 and 0.2, each with interladder 
couplings J" /J = -0.2, -0.1, 0.1 and 0.2£a yielded 
xVDOF = 1.17 and a,ms = 0.000302. As in the fit to the 
trellis layer QMC data for J'/ J < 1 above, the data could 
be fitted well assuming that the spin gap is independent 
of the interladder coupling. The set of {N2nj^ D2nj} pa- 
rameters obtained, including the constrained -D21, D22 
and D23 series, is given in Table pcg_ and the fit is plotted 
as the set of solid curves in FigsT^ and ||, respectively. 
The fit deviations are plotted vs temperature in Fig. |6^. 
The series for D2n with n = 2 an d 3 have the following 
terms in addition to those in Eq. (A7) 



D2niidd = D2n3[^-jjJ + D2ni\^-p J +D2n5y-jjJ 
+ D2ns[jj) +D2n9[j;) 
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1' . 



D2n\^{^-^ + D2n\i{-p^ 
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'J"\2 



/ J \ 6 



D2n20y--pj + D2n2iy--p 
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D2n23 ( j7 j + D2n24y-JJ 



J\2 



D 



2n25 



f- 



/ J" \ 4 

D2n26 y~jr) 



(A8) 



The 27-parameter three-dimensional fit to the 119 
stacked ladder data points for J' /J = 0, 0.1 and 0.2 
and J'" /J = 0.1 and 0.2 yielded xVDOF = 0.81 and 
frms = 0.000208, and is plotted as the set of solid curves 
in Fig. O. The fit deviations are plotted vs temperature 
in Fig. The set of {N^njj D^nj} parameters are given 



TABLE X. Parameters in Eq. (A5) obtained by fitting quantum Monte Carlo simulations of x*(^) LaCu02.5-typc 3D 
coupled 5 = 1/2 Heisenberg ladders with no spin gap. 



n 


N„o 


iVnl 


N„2 


N„3 


Nni 


iVnS 


1 


-0.4102143 


0.6183598 





0.1499184 


-16.56455 


3.798390 


2 


0.2411046 


-0.2644301 





0.2756975 


-0.03651210 


1.150476 


3 


-0.0318917 


0.09842195 


-0.05803464 


-0.01381445 


-0.5306434 


0.4387261 


5 


0.0008320664 


-0.002259251 


0.001012209 





0.03673690 


-0.005869753 


n 




Dm 


D„2 


D„3 


D„4 


D„5 


1 


0.08978571 


0.8683598 
Die: 0.5 





0.1499184 


-16.56455 


3.798390 


2 


0.2859975 


-0.05780379 


0.2795899 


0.6598366 


-8.068789 


-1.016509 






D26- -0.2051071 


D27: 0.03747960 


Dag: -8.282277 


D29: 1.899195 


D210: 0.9495975 


3 


0.06944034 


0.1810730 


-0.1754189 


0.02929887 


-4.792592 


2.246872 






Dae: 0.1205523 


D37: 0.06892439 


Das: 0.06507728 


D39: 0.6127174 


D310: -1.782950 






D311: 0.1189616 


D312: 0.01873980 


D313: -4.141139 


D3i4:0.4747988 


D315: 0.9495975 


4 


0.01957895 


0.03025742 


-0.007991254 


0.1063007 


-1.631987 


0.3758439 


6 


-0.0002842768 


0.001377244 


-0.001162417 





0.09843710 


0.0408774 
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TABLE XI. Parameters in Eq. (A7) obtained by fitting quantum Monte Carlo simulations of x*(*) for isolated 5=1/2 
Heisenberg ladders, trellis layers and stacked ladders with J' / J > 1, all of which have spin gaps. 
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FIG. 61. Semilog plot versus temperature T of the devia- 
tion of the fit from the 162 QMC x* (*) data points for 5 = 1 /2 
spatially anisotropic Heisenberg two-leg ladder trellis layers 
with intraladder couplings J /J' = 0.1 and 0.2, and interlad- 
der intralayer couplings J" /J' = -0.2, -0.1, 0.1 and 0.2 (•), 
for which there is a spin gap. The error bars are the esti- 
mated accuracies of the QMC data. The fit function for the 
QMC data is given in Eqs. and (A7) with parameters in 
Table Ixil 
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FIG. 62. Semilog plot versus temperature T of the de- 
viation of the 29-parameter three-dimensional fit from the 
119 QMC data points for stacked 5=1/2 spatially 

anisotropic Heisenberg two-leg ladders with intraladder cou- 
plings J/ J' — 0, 0.1 and 0.2, and interladder couplings 
J'" / J' = 0.1 and 0.2 (•), for which there is a spin gap. The 
error bars are the estimated accuracies of the QMC data. The 
fit function for the QMC data is given in Eqs. (0) and 
with parameters in Table Kfl. 



in Table ^CJ The constrained series for with n = 2 
and 3 h ave the following terms in addition to those in 
Eq. (0) 



£1 
J' 



3n8 



J_\3 
J' 



~J' 



/ J \ 31 / J'" \ 3 

/ J x 10 

/ J'" \\ /J\5 /J\6 

+ (-) ^-"W 
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i„19(j7) 



1' 



D3n2Q 
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1' 



J_\4 
J' 



J 



D 



J\5 



3n21 



-D3n23 1 j7 j + ^3ri24 y-JJ 



J\2 



D 



3n25 



f- 



f J'" \ 4 
^3n26 \~jr) 



(A9) 



6. 5 = 1/2 Three- and Five-Leg Ladders 

In these n-leg ladders, the spin gap A = 0. The spins 
are not all equivalent, since the outer leg spins have three 
nearest neighbors whereas the inner leg spins have four. 
Therefore, only the Curie- Weiss terms to 0{\/t^) in the 
HTSE in Eqs. (|ri|) and (PX can be utilized to constrain 
the fits. Frischmuth et auH^ have carried out QMC sim- 
ulations of X*{t) for ladders with n = 1-6 and spatially 
isotropic exchange and also for n = 3 with anisotropic 
exchange. 

For the three-leg [n = 3) ladders with anisotropic, 
exchange, we fitted the 126 QMC x*{t) data pointsEZE 
for J'/J = 0.4, 0.5, 0.6 0.7 and 1 by Eqs. (0) using 
jmax = A = 0, the constrained Di given by Eq. ( |l8a| ), 

Adi given by the average value (2/3)(3-|-2 J'/ J), and v'^^}^ 
in Eq. (p7|). The error bars for J' / J — 1 (which were not 
available to us) were estimated from those of the other 
data sets. The -/V„ and Dn parameters are written as 
power series in J'/J according to 

= Nno + iVi„ (y ) + N2n (y ) ' , (n = 1 - 4) 

D„ = A.o + Din[j) + D2n (y ) ' . (n = 1 - 5) (AlO) 
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FIG. 63. Semilog plot versus temperature T of the devi- 
ation of the fit from the 126 QMC x*(*) data points for 
5=1/2 spatially anisotropic three-leg Heisenberg ladders 
with J'/J = 0.5, 0.5, 0.6, 0.7 and 1 (•), for which there is 
no spin gap. The error bars are the estimated accuracies of 
the QMC data. The f it function for the QMC data is given 
in Eqs. ( p^j ) and (AlC) with parameters in Table XII 



The {Nn, Dn} series coefficients are given in Table XI]. 
The two-dimensional fit is shown as the set of solid curves 
through the data in Fig. ^ where extrapolations to t = 2 
and interpolation curves for J'/J = 0.8 and 0.9 are also 
shown. The deviations of the fit from all the data are 
plotted vs temperature in Fig. The qualities of the 
fit are xVDOF = 1.08, (Trms = 0.000111, and the relative 
rms deviation is 0.227%. 

We also obtained^n unweighted fit to the x*{t) data of 
Frischmuth et alJtl3 for the isotropic five-leg ladder (for 

which the error bars were not available to us) using 'P^^^ 



in Eq. (y), A = 0, Di given by Eq. ( |l8a| ) and Mi su 

by the average coordination number 4 — (2/n) — 
with parameters 

iVi = -0.2732853 , N2 = 0.09333487 , 
N3 = 0.006660300 , Di = 0.6267147 , 
D2 = 0.3077097 , L»3 = 0.04438012 , 



Da = 0.07488932 



DOF 



5.6 X 10" 



(All) 
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